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INTRODUCTION 

The  objective  of  this  research  is  to  exploit  signal  processing  methods  to  improve  ultrasound  resolution  in  the 
detection  and  classification  of  small  breast  lesions.  Ultrasound  is  routinely  used  in  conjunction  with  conventional 
x-ray  mammography  for  the  evaluation  of  variety  of  breast  masses,  and  it  is  particularly  useful  for  the 
differentiation  of  solid  from  cystic  breast  lesions.  However,  ultrasound  resolution  is  limited  and  the  major 
obstacle  for  the  achievement  of  high  resolution  is  acoustic  wavefront  distortion  which  is  developed  as  coherent 
waves  propagate  through  the  inhomogeneous  breast  [1-6].  The  primary  distortion  effect  is  the  scattering  of 
coherent  energy.  As  a  result,  the  target  strength  is  reduced,  the  image  lobe  is  broadened  [1-2],  the  background 
level  is  increased  and  therefore  the  image  contrast  is  reduced  [7].  The  secondary  distortion  effect  is  the  multipath 
interference  which  produces  ghost  image  artifacts  or  false  targets  in  addition  to  true  targets  in  the  image  [3]. 
Adaptive  compensation  of  wavefront  distortion  is  an  active  area  of  current  research  [8].  Two  basic  phase 
deaberration  algorithms  have  been  developed  in  [9-11].  These  algorithms  adjust  time  delays  of  received 
waveforms  at  different  elements  of  a  phased  array  to  compensate  delay  errors  due  to  an  inhomogeneous  medium 
and  consequently  reduce  the  scattered  energy.  These  phase  deaberration  algorithms  are  very  useful  to  partially 
remove  phase  distortion  and  build  up  the  strength  of  the  coherent  field.  Wavefront  amplitude  distortion  remains, 
as  does  residual  phase  distortion,  and  produces  a  significantly  high  background  level  in  the  image.  Measurements 
of  44  in  vivo  breasts  obtained  from  the  Hospital  of  the  University  of  Pennsylvania  have  shown  that  this 
background  level  in  the  image  caused  by  the  amplitude  distortion  is  17  dB  higher  than  that  obtained  from  a 
homogeneous  medium  [5]. 


The  simple  models  of  Fig.  1  illustrate  the  effects  of  wavefront  distortion  upon  a  point  source  image  or  its 
angular  spectrum  (intensity  distribution),  and  the  effects  of  phase  deaberration.  Without  medium  induced 
distortion,  the  diffraction-limited  coherent  spectrum  is  shown  in  Fig.  la.  With  medium  induced  scattering,  the 
overall  angular  spectrum  is  the  sum  of  coherent  and  scattered  spectra  [12]  as  shown  in  Fig.  lb.  Because  of 
scattering,  the  spectrum  is  broadened  and  reduced  in  strength.  The  background  level  is  increased  and,  therefore, 
image  contrast  is  reduced.  A  phase  deaberration  process  removes  some  of  phase  distortion  and  therefore 
strengthens  the  coherent  spectrum  as  shown  in  Fig.  Ic;  yet  a  significantly  high  incoherent  background  level 
remains  in  the  spectrum  due  to  residual  phase  errors  and  wavefront  amplitude  distortion.  Any  phase-deaberration 
procedure  that  provides  a  better  phase  error  estimate  would  reduce  this  background  level  [13].  Algorithms  that 
compensate  wavefront  amplitude  distortion  can  further  reduce  this  background  level  [13]  .  In  multipath 
interference  (Fig.  Id)  the  scattered  energy  removed  by  the  phase  deaberration  process  strengthens  the  primary 
spectrum  (sharp  peak)  by  folding  the  removed  scattered  energy  into  the  primary  field.  The  result  is  improvement 
in  the  ratio  of  primary  to  interference  lobe  level  (Fig.  le). 


diffraction 

limited 

spectrum 


e 


Coherent  spectrum 

> 

Scattering  spectrum 
due  to  amplitude 
and  phase  distortions 


Coherent  spectrum 


Scattering  spectrum 
due  to  amplitude 
and  residual  phase 
distortions 


(a) 


(b) 


Multipath 
interference 
spectrum 

Scattering  spectrum  Ac  _ 
due  to  amplitude  gj, 
and  phase  distortions 


Coherent  spectrum 


Scattering  spectrum 
due  to  amplitude 
and  residual  phase 
distortions 


(d)  (e) 

Fig.  1.  Statistical  models  of  angular  source  intensity  distribution,  (a)  diffraction  limited  spectrum,  (b)  Scattering  induced 
distortion.  The  spectrum  is  the  sum  of  coherent  spectrum  and  scattering  spectrum,  (c)  After  phase  compensation  of  (b).  (d) 
Multipath  interference  and  scattering  induced  distortion,  (e)  After  phase  compensation  of  (d). 
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Full  compensation  of  wavefront  distortion  presents  a  more  difficult  problem  because  there  is  no  general  way 
to  reduce  wavefront  amplitude  distortion.  Liu  and  Waag  [14]  have  improved  time-delay  compensation  by 
backpropagating  the  received  waveforms  to  an  optimal  distance  and  then  performing  time-delay  compensation  at 
this  distance.  Backpropagation  partially  compensates  waveform  distortion  produced  by  wave  propagation  from 
the  backpropagation  distance  to  the  receiving  aperture.  We  have  developed  a  novel  spectrum  amplitude 
compression  technique,  we  call  toward  inverse  filtering  (TIF),  during  the  first  two  years  of  this  project  for  this 
purpose  [7,13]  (see  Appendices  A  and  B).  Spectrum  amplitude  refers  to  the  complex  modulus  at  each  discrete 
frequency  within  the  signal  band.  TIF  ,when  used  in  combination  with  phase  deaberration  and  backpropagation, 
provides  a  general  means  for  reducing  wavefront  amplitude  distortion  and  therefore  reducing  the  incoherent 
background  level  in  the  image. 

High  resolution,  high  contrast  breast  imaging  implies  seeing  a  black  cyst  60  dB  weaker  than  a  brightly 
speckled  target.  This  implies  that  30  dB  contrast  has  to  be  achieved  in  one-way  transmission  and  reception. 
Time-delay  algorithms  have  been  applied  to  in  vivo  breast  data  measured  from  one-dimensional  phased  arrays 
[2,5]  and  the  improvement  in  image  contrast  is  not  significant.  The  ideal  adaptive  compensation  must  be  a  two 
dimensional  (2-D)  problem  of  capturing  the  entire  2-D  wavefront  and  constructing  a  2-D  compensation  vector  for 
the  array  [15-16].  In  the  last  two  years,  we  found  that  by  applying  time-delay  correction  to  in  vitro  one-way  2-D 
breast  data,  the  algorithms  suppressed  scattered  energy  and  folded  it  into  the  target  image,  restoring  the  mainlobe 
diffraction  shape  to  the  -20  dB  region.  [6,17].  We  also  found  that  an  additional  4-5  dB  improvement  can  be 
obtained  by  applying  Liu  and  Waag's  algorithm  [13].  Additionally,  with  TIF,  we  are  able  to  improve  all  the 
correction  algorithms  and  restore  the  mainlobe  diffraction  shape  down  to  -30  dB  [13]  which  is  the  necessary  one¬ 
way  condition  to  image  small  breast  lesions  and  lesion  architecture. 

There  are  two  remaining  tasks  that  need  to  be  completed  before  we  can  apply  wavefront  compensation 
algorithms  to  clinical  patients.  First,  we  have  demonstrated  30  dB  suppression  of  the  distortion  sidelobes  in  one¬ 
way  propagation  experiments,  and  we  need  to  further  demonstrate  that  TIF  is  useful  in  two-way  pulse-echo 
imaging.  Second,  2-D  compensation  requires  2-D  transducer  arrays  which  are  not  commercially  available. 
Therefore,  we  need  to  develop  a  2-D  array  to  implement  wavefront  compensation  algorithms. 

In  our  original  proposal,  we  proposed  to  design  a  2-D  concave  transducer  which  would  reduce  breast- 
curvature  caused  refractive-distortion  effects  upon  images  in  addition  to  being  a  tool  for  implementing  2-D 
compensation  algorithms.  However,  concave  shape  introduces  complications  in  array  design  and  fabrication. 
Since  breast-curvature  caused  distortion  effects  can  be  largely  reduced  by  using  a  planar  transdueer  compressing 
against  the  examined  breast,  we  have  focused  our  research  on  design  and  fabrication  of  such  a  2-D  planar  array. 

In  this  report,  we  present  our  preliminary  experimental  results  of  wavefront  compensation,  using  two-way 
pulse-echo  phantom  data.  We  also  give  our  progress  report  on  the  development  of  a  2-D  array.  Deabrration 
results  are  reported  in  Part  A  and  array  development  is  in  Part  B. 

In  addition  to  ultrasonic  wavefront  deaberration  and  2-D  array  fabrication  work,  we  have  done  new  work 
on  coupling  ultrasound  and  near  infrared  diffused  light  imaging  methods.  The  objective  is  to  increase  detection 
sensitivity  and  specificity  of  breast  lesions.  Experimental  results  with  phantoms  are  very  encouraging,  and  clearly 
indicate  that  this  novel  method  could  be  potentially  very  important.  This  work  was  not  proposed  in  the  original 
proposal,  but  it  is  directly  related  to  Pi's  breast  cancer  research.  Therefore,  it  is  reported  in  Part  C  and  the  preprint 
submitted  for  publication  is  in  Appendix  C. 

PART  A.  WAVEFRONT  COMPENSATION  IN  PULSE-ECHO  USING  TIF 
I.  EXPERIMENTS  AND  DATA  PROCESSING 

The  two-way  pulse-echo  data  were  down  loaded  from  the  world-wide  web  at  http://bul.eecs.umich.edu 
and  the  detailed  description  of  the  data  can  be  found  in  [Al].  Briefly,  all  data  sets  were  acquired  with  a  1-D,  64- 
element,  15.4-mm  array  operating  at  3.33  MHz  with  a  fractional  bandwidth  of  about  35%.  Each  complete  data 
set  consisting  of  4096  RF  A-scans  sampled  at  a  rate  of  17.66  MHz  with  10  Bit  analog-to-digital  conversion. 
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‘Each  rfecord  represents  a  single  element  acting  as  a  transmitter  and  an  independent  element  acting  as  a 
receiver.  Three  data  sets  were  obtained  from  a  standard  graphite-gel  AIUM  resolution  phantom  with  three 
aberration  plates,  made  of  RTV  silicone  rubber,  placed  between  the  transducer  array  and  the  phantom. 
The  first  plate  introduces  only  a  very  small  defocusing,  and  hence,  is  used  as  the  control  plate  or  reference 
plate.  The  second  plate  causes  phase  distortions  less  than  n,  while  the  third  plate  produces  phase 
distortions  within  2  ;r.  In  our  preliminary  experiments,  we  have  selected  first  and  third  data  sets  and  refer 
to  these  data  sets  as  Plate  0  and  Plate  2. 


B-scan  images  without  any  correction  were  formed  using  conventional  shift-and-add  beamforming 
process  upon  these  two  data  sets.  The  simple  time-delay  phase  deaberration  algorithm  [A2]  was  used  to 
correct  phase  errors  of  Plate  2.  TIE  was  implemented  after  time  domain  phase  correction  and  the  steps 
were  shown  in  the  following  diagram  (right).  First,  Fast  Fourier  Transforms  (FFT)  of  temporal 
waveforms  at  each  element  position  of  aperture  were  taken  within  a  selected  range  window.  Let  complex 

70  (n) 

signal  at  element  position  n  of  discrete  frequency  bin  number  m  be  represented  as  >l„,(n)e  ”*  .  TIF 

Nth  root  operation  was  performed  upon  amplitude  A^{n)  and  the  resulting  complex  signal  was 


Only  square  root  was  used  in  our  experiments,  beeause  of  signal-to-noise  ratio 
consideration.  Inverse  FFTs  of  (n)e'^'^'"^”^  at  each  element  position  were  taken  to  obtain 


processed  RF  signals.  B-scan  images  were  formed  using  shift-and-add  upon  processed  RF  signals. 
All  images  were  formed  with  fixed  focusing  on  transmission  and  reception. 


Fig.al.  Diagram  of  image  formation  with  time-delay  correction  (left)  and  TIF  correction  (right). 
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II.  PREUMINARY  RESULTS 

Fig,a2  is  the  B-.scun  image  of  Plate  0  and  is  u.sed  a.s  the  reference  image.  Several  bright  targets 
(wires)  and  two  anechoic  regions  ("cysts")  are  shown  in  the  figure.  Wire  targets  allow  u.s  to  evaluate  2- 
way  system  point-spread  function  and  cysts  allow  us  to  evaluate  system  contrast  resolution.  Since  cysts 
and  most  tumors  appear  anechoic,  deaberralion  procedures  should  make  them  more  "black"  in  order  to 
identify  them.  A  portion  of  B-.scan  image.s  of  Plate  2  is  shown  in  Fig.ei3.  Because  of  the  limitation  of 
computing  power,  only  the  portion  which  contains  the  wire  targets  imd  the  lower  ancchoic  region  in  Fig.ft2 
is  computed  and  shown  in  Fig.a3.  The  images  shown  in  Fig.3a  are  results  of  no  correction  (part  a),  with 
phase  correction  (par  b)  and  TTF  correction  (part  c).  Comparing  parts  a  and  b.  we  can  see  that  the  system 
point-spread  function  is  improved  by  phase  deaberration.  The  sidelobes  of  the  close  by  wire  target  as 
pointed  by  upper  airows  in  both  images  arc  improved.  However,  the  contra.st  improvement  of  anechoic 
region  is  marginal.  With  TIF  correction  (part  c),  the  anechoic  region  is  beginning  to  show  up.  There 
results  arc  first,  and  clearly  indicate  that  TIF  offers  better  improvements  on  low  contrast  lesions. 

III.  FUTURE  WORK 

We  will  furiher  quantify  the  improvement  brought  about  by  TIF  by  measuring  the  contrast  ratio  of 
energy  inside  and  outside  ancconic  regions.  In  addition,  theoretical  work  will  be  undertaken  to  study  the 
optimal  compression  parameter  of  TIF.  Finely,  the  amplitude  distribution  as  a  function  of  the  TiF 
compression  parameter  will  he  .studied.  Statistics  of  amplitude  distribution  are  u.seful  for  tissue 
characterization. 
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Fig.a2.  B-Scan  image  of  Plate  0.  Wires  are 
as  anechoic  regions. 
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Fig.aSa.  A  portion  of  B-scan  image  of  Plate  2.  Upper  arrow  points  to  the  sidelobes  of  the  close  by  wire  target 
and  lower  arrow  points  to  the  anechoic  region. 

Plata  2,  with  phasa  corraction 

50 

100 
150 

200 

250 

300 

350 

400 

450 
500 

60  100  150  200  250  300  350  400  450  500 

Fig.aSb.  A  portion  of  B-scan  image  of  Plate  2  after  phase  deaberration.  Upper  arrow  points  to  the  sidelobes  of 
the  close  by  wire  target  and  lower  arrow  points  to  the  anechoic  region.  The  sidelobe  level  has  been  improved  but 
the  improvement  of  contrast  of  anechoic  region  is  marginal. 
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Plate  2,  with  TIF 
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Fig.aSb.  A  portion  of  B-scan  image  of  Plate  2  after  TIF.  The  arrow  points  to  the  anechoic  region.  The  contrast  of 
anechoic  region  is  improved  and  the  cyst  beginning  to  show  up. 
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PART  B.  DESIGN,  FABRICATION  AND  TESTING  OF  A  TWO-DIMENSIONAL  ULTRASONIC 
ARRAY 

I.  INTRODUCTION 

The  goal  is  the  development  of  two-dimensional  (2-D)  piezoelectric  array  transducers  operating  at  3.75 
-10  MHz.  The  techniques  are  developed  and  refined  in  stages.  A  2-D  array  of  2.5  cm  x  2.5  cm  size  operating  at 
1  MHz  is  first  built.  The  transducer  head  contains  the  transducer  layer  with  ground  metalization  on  the  front  side 
and  individual  electrodes  on  the  back  side.  Connections  to  the  back  electrodes  of  all  operational  piezo-elements  are 
provided  through  wires  soldered  to  a  PC  board.  A  quarter  wavelength  impedance-matching  front  layer  made  of 
silver-filled  epoxy  is  glued  to  the  ground  metalization  side  of  the  transducer  layer.  Design,  fabrication  and 
preliminary  testing  of  the  fabricated  array  are  given  below. 

II.  DESIGN  OF  THE  ARRAY 

A.  Transducer  Materials  and  Interconnects 

A  1-3  type  of  piezoelectric  composite  consisting  of  PZT  (lead  zirconate  titanate)  ceramic  rods  in  a  polymer 
matrix  serves  as  the  material  from  which  two  dimensional  addressable  transducer  arrays  are  constructed  [Bl]. 
Composites  (combinations  of  two  distinct  solid  phases  with  different  connectivity  patterns)  can  greatly  enhance 
material  utilization  and  properties  by  the  synergy  generated,  the  properties  of  the  combination  being  better  than  the 
individual  properties.  Composites  provide  a  much  wider  range  of  material  properties  than  those  of  the  standard 
ceramic  and  polymer  transducers  and  are  being  actively  studied  in  the  medical  ultrasonic  field  PB2]. 

In  our  design  the  individual  rods  with  end  faces  individually  electroded  serve  as  transducer  elements. 
Ideal  characteristics  are  strong  electromechanical  coupling  for  efficient  energy  conversion,  low  acoustic  impedance 
for  good  transmission  into  tissue,  and  a  reasonably  large  dielectric  constant  for  efficient  electrical  impedance 
matching  [B3].  Further  advantages  of  strictly  periodic  composite  materials  are  control  of  lateral  transducer 
resonances  and  reduction  of  crosstalk  between  individual  electrodes  in  transducer  arrays. 

By  utilizing  the  optimized  1-3  composite  transducer  design,  high  transmit  efficiency,  high  receive 
sensitivity  and  broad  bandwidth  are  achieved.  Proper  choice  of  transducer  material  and  component  dimensions 
optimizes  electromechanical  coupling  and  minimizes  mechanical  cross-talk  and  spurious  low  frequency  lateral 
modes  [B2],  [B3-B6].  We  note  that  the  1-3  composite  exhibits  -30%  higher  thickness  mode  electromechanical 
coupling  constant,  k3t,  than  bulk  ceramic  and  that  individual  pixel  k3t  in  the  1-3  composite  approaches  the  free 
k33  of  an  "ideal"  long  slender  unconstrained  piezo  element  [B2].  The  monolithic  bulk  sample  exhibits  self- 
pinning  and  hence  has  k3t  <  k33. 

B.  Design  of  Front  and  Back  Layers 

Fig.bl  is  a  sectional  view  of  the  array  configuration  along  with  front  and  back  impedance-matching  layers. 
The  transducer’s  front  surface  is  interfaced  with  a  uniform  XIA  impedance-matching  front  layer.  A  patterned 
composite  backing  layer  of  polymer  and  microfabricated  ceramic  pillars  [B7]  can  provide  the  piezoelements  with 
the  appropriate  matching  mechanical  impedance,  thus  allowing  its  broadband  operation  with  minimal  signal 
degradation.  The  backing  layer  is  patterned  to  match  the  active  elements  of  the  transducer.  The  material  in  the 
backing  layer,  which  is  between  the  positions  of  the  operative  piezo-elements  is  made  of  a  compliant  material 
(such  as  polyurethane)  so  as  to  maintain  the  mechanical  independence  of  the  individual  ceramic  elements. 

The  X I A  plate,  covering  the  front  surface  of  the  metallized  transducer,  can  be  fabricated  from  a  composite 
material.  Commercial  impedance  matching  materials  such  as  carbon  and  glass  filled  phenolic  or  silver-filled 
epoxy  are  available  and  can  be  machined  to  the  appropriate  thickness  and  glued  in  place.  The  patterned  back  layer 
can  be  constructed  from  a  matching  inactive  array  obtained  from  the  same  1-3  composite  supplier.  The  back  layer 
microfabrication  process  will  then  consist  of  inserting  conducting  paths  alongside  the  matching  pillars  or 
columns.  They  provide  connection  from  the  metallized  electrodes  to  the  commercial  multiplexer  boards. 
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matching  PZT  Impedance 

Fig.bl.  Sectional  view  of  array  configuration  plus  front  and  back  impedance  matching  layers. 


m.  FABRICATION  OF  THE  ARRAY 

Our  fabrication  of  the  array  with  a  front  layer  and  a  back  layer  were  done  in  a  sequence  of  prototypes.  Our 
first  one  was  the  array  with  back  electrodes  and  ground  plane  metalization.  The  next  step  took  this  same  first 
prototype  and  glued  a  quarter  wavelength  matching  front  layer  made  of  silver-filled  epoxy  to  the  gold  ground 
metalization.  The  final  step  will  add  a  matching  back  layer  with  patterned  conductive  vias  behind  the  transducer 
layer.  Electronic  and  acoustic  testing  has  been  performed  at  each  stage.  A  notable  change  in  the  SNR  should  be 
evident. 

A.  PZT  Composite  Materials  and  Metallization 

We  have  developed  the  procedures  for  making  the  final  transducer  array  package  by  a  series  of  steps.  The 
arrays  utilize  PZT  piezoelectric  ceramic -polymer  1-3  composite  materials  obtained  from  Materials  Systems, 
Inc.(MSI)  [B 1].  Our  design  utilizes  the  MSI  periodic  arrays  of  ceramic  rods  within  a  polymer  matrix  to  construct 
a  2-D  array  with  element  size  0.75  X  and  element  spacing  0.5  X  (see  Fig.b2).  The  array  materials  were  obtained 
as  2.5  cm  x  2.5  cm  modules.  Using  DC  Gold  Sputtering  metalization,  carried  out  in  our  Dental  School  facilities 
[University  of  Pennsylvania],  electrodes  were  applied  to  the  entire  front  surface  to  serve  as  the  ground  plane.  DC 
Sputtering,  which  occurs  at  room  temperature,  was  used  because  the  polyurethane  (which  holds  together  the  piezo 
rod  array)  swells  at  temperatures  above  100  degrees  Celsius.  This  was  deduced  from  our  previous  preliminary 
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'  array  with  only  12  transducers  which  resulted  with  an  open  circuit  on  the  ground  plane.  The  thin  electrode  of 
about  300  angstroms  of  gold  was  thickened  through  Gold  Electroplating  at  65  degree  Celsius,  so  as  to  have  at 
least  1000  angstroms  of  Gold.  The  non-grounded  plane  of  the  array  had  individual  contacts  on  selected  array 
elements.  These  contacts  were  manually  painted  silver  epoxy  pads  which  matched  the  silver  epoxy  that  was  used 
to  bond  the  individual  wires  to  the  pads. 


Fig.  b2.  2.5  cm  x  2.5  cm,  IMHz  transducer  module. 

The  module  shown  in  the  picture  is  enlarged  by  2.  The  module  consists  of  306  elements. 

B.  Bonding  and  Interconnections 

A  7  cm  by  7  cm  PC  board  with  40  microns  of  copper  metalization  on  both  sides  was  used  to  provide 
mechanical  support  to  the  individual  connections  of  selected  array  elements.  From  the  PC  board  soldered  wires 
were  connected  to  the  array  of  piezo  rods  (see  Fig.b3).  Since  the  piezo-composite  fabrication  process  of  MSI 
yielded  a  quasi-periodic  element  arrangement,  the  photolithographic  process  masks  for  the  etching  of  the  copper 
metalization  of  the  PC  board  were  obtained  from  digitizing  actual  piezo-composite  array  photos  [B9],  a  technique 
we  have  used  in  the  past.  This  allowed  for  ohmic  patterns  which  correspond  to  the  piezo  rods.  The  copper 
pattern  on  the  front  side  of  the  PC  Board  consists  of  .5  mm  squares  (which  locations  correspond  to  the  piezo 
rods)  with  lines  to  pads  (which  correspond  to  a  54  pin  socket).  The  design  of  the  ohmic  lines,  not  including  the 
contacts  for  the  pin  socket,  can  be  seen  on  Fig.  b3.  The  side  view  of  the  copper  pattern  and  its  location  on  one 
side  of  the  PC  board  can  be  seen  on  Fig.  b4.  The  squares  have  0.25  mm  holes  which  are  filled  with  a  copper  wire 
(0.13mm  in  diameter)  which  stick  out  to  a  length  which  would  contact  the  piezo  rods.  The  only  copper  pads  that 
do  not  have  lines  are  four  center  elements  which  are  to  be  connected  to  individual  copper  wires.  On  the  backside 
of  the  board  (see  Fig.  b4,  #5)  there  are:  a  copper  pattern,  spacers,  and  locating  pins.  The  copper  pattern,  which 
consists  of  filled  in  circles  (which  was  possible  with  the  digitized  picture  of  the  piezo  rods  and  photolithographic 
processes),  served  as  additional  support  for  the  wires  of  center  receiving  elements  and  for  wires  that  needed  solder 
repair,  more  accurate  positioning  of  the  locating  pins,  and  visual  confirmation  that  the  wires  would  accurately 
connect  to  a  piezo  rod.  Spacers  were  soldered  wires  which  provide  evenly  spaced  pressure  points  needed  when 
connecting  the  board  to  the  array.  Locating  pins  are  soldered  copper  wires  (0.8mm  in  diameter)  into  1.1  mm 
holes  on  the  PC  board.  The  array  has  six  matching  holes  corresponding  to  the  farthest  elements  from  the  center 
which  are  not  in  use.  (See  Fig.  b5  for  the  cross  section  view  of  array  assembly.) 

C.  Electronic  Connections 

There  were  54  receiving  elements  which  were  connected  to  a  pin  socket.  This  allows  for  the  manual 
changing  of  each  receiving  element  by  inserting  a  wire  that  is  soldered  to  a  coaxial  cable.  Future  work  includes  a 
multiplexer  with  built  in  preamplifier.  Four  transmitting  elements  were  chosen  from  the  center.  These  had 
individual  wires  which  plug  into  a  contact  to  a  coaxial  cable. 

IV.  TRANSDUCER  ARRAY  TESTING 
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'  A.  Elefctrical  Testing 

Selected  elements  of  the  array  were  characterized  throughout  the  design  process  to  provide  data  for 
optimization  of  interfacing  between  the  active  piezoelectric  elements  and  the  first  stage  of  electronics.  These 
measured  characteristics  are  important  as  the  overall  performance  of  the  imaging  system  is  crucially  dependent  on 
the  sensitivity  of  the  front-end  design  and  its  SNR. 


Fig.bS.  The  array  of  squares  connected  to  lines  are  the  receiving  elements.  The  four  center  squares  without  lines 
are  transmitting  elements.  Six  of  the  surrounding  squares  without  lines  are  locations  of  the  position  pins. 


Sideview  of  Transducer  Connections 


Fig.b4.  The  ohmic  contacts  for  the  individual  elements  are  illustrated  by  the  above  picture  and  the  below 
description.  1.  four  wires  to  transmitting  transducers.  2.  54  pin  socket.  3.  copper  pattern  on  the  PC  board.  4. 
PC  board  with  wires,  spacers  and  locating  pins.  5.  backside  of  PC  board.  6.  array  holes  for  6  locating  pins.  7. 
ground  metalization  side  of  array. 


Optimization  of  SNR  involves  knowledge  of  the  piezoelectric  element  impedance  and  information  on  the 
possible  changes  in  the  impedance  due  to  ohmic  contacts  of  the  piezoelectric  structure  using  different  auxiliary 
materials.  Similar  measurements  should  also  be  made  for  the  final  array  prototype  with  selected  individual 
elements  appropriately  connected. 

On  an  Impedance  Meter  each  of  the  57  elements  were  tested  for  impedance  from  0.1  to  19.6  kHz  in 
increments  of  0.5  kHz.  52  (3  transmitting  and  49  receiving)  elements  demonstrated  uniform  impedance  ranging 
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'  from  9.77  to  0.341  Mohms  (see  Fig.b6).  One  of  the  57  elements  served  as  a  control  for  bad  contacts,  because  it 
did  not  have  a  wire  connecting  from  the  PC  board  to  the  array.  Four  transducers  had  high  impedance  similar  to 
the  control  (e.g.  at  lOkHz  ranging  from  2.2  to  ISMohms).  Thus  it  can  be  deduced  that  the  difference  in 
impedance  is  due  to  bad  contacts  on  the  wire.  Of  the  intended  57  total  elements,  52  proved  to  have  good  contacts. 
The  average  capacitance  was  23.1  pF. 


Cross  Section  of  Array  Assembly 


Fig.b5.  The  array  assembly  consists  of  a  PC  board  with  copper  wires  soldered  to  copper  contacts,  soldered 
spacers  and  soldered  locating  pins.  On  the  top  side  of  the  PC  board  are  copper  contacts  and  lines  that  lead  to  a  54 
pin  socket.  Locating  pins  fit  into  holes  in  the  array.  Copper  wires  are  bonded  to  the  rod  metalization  with  silver 
epoxy.  Spacers  prevent  bent  wires  and  short  circuits  and  bond  to  the  array  with  silver  epoxy  for  mechanical 
support. 


Fig.bb.  The  impedance  (1/ljwCl)  of  52  out  of  57  elements  were  uniform  and  of  the  order  of  Mohms  (<0.4Mohm 
for  f>20kHz).  The  standard  deviation  is  drawn  by  error  bars. 
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B.  Acoustic  Testing 

Fig.b7  is  a  block  diagram  of  our  water  tank  testing  system.  The  first  tests  were  made  in  transmission, 
where  each  selected  transducer  element  acts  as  a  transmitter  and  a  needle  hydrophone  (0.6  mm  in  diameter)  as  a 
receiver.  An  ultrasound  pulser  (Reflectoscope)  generates  RF  pulses  that  drive  the  selected  element  one  at  a  time. 
The  received  RF  signals  are  amplified  by  the  low  noise  preamp  and  amp  inside  the  pulser.  A  digital  scope  (Tek 
2230)  samples  the  16  bit  RF  data  at  500  MHz  which  is  well  above  the  Nyquist  frequency.  The  digital  data  are 
transferred  to  PC  through  a  GPIB  board. 


water  2-D  GPIB 


Fig.bV.  Block  diagram  of  water  tank  testing  system. 


Fig.bS  is  a  typical  impulse  response  of  a  single  element  without  front  and  back  layers.  The  SNR  is  24  dB 
which  is  very  good.  SNR  should  be  further  improved  when  the  front  and  back  layers  are  added  to  the  transducer 
layer.  As  of  the  writing  of  this  report,  a  quarter  wavelength  front  layer  has  been  glued  to  the  transducer  array  but 
the  contact  is  not  good  for.  making  measurements.  We  need  to  reglue  the  front  layer  with  better  contact. 


Fig.bS.  Impulse  response  of  a  single  element. 


V.  SUMMARY 

From  this  preliminary  prototype  we  deduced  that  the  patterned  PC  board  with  soldered  wires  worked  as  a 
good  mechanical  support  of  ohmic  lines  and  provided  accurate  contact.  The  few  that  were  detected  to  have  bad 
contact  can  easily  be  prepared  by  replacing  the  wire  and  resoldering.  The  PC  board  can  thus  be  reused  with 
different  arrays  of  the  same  make  by  removing  the  present  array.  Removing  the  present  array  will  allow  us  to 
glue  a  back  layer  between  the  PC  board  and  the  array.  This  removal  is  possible  because  the  array  is  held  together 
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■  with  sliver  epoxy.  When  force  is  applied,  the  epoxy  will  separate  from  the  array.  The  next  step  involves  using 
this  same  prototype  at  first,  however,  a  back  layer  will  be  added  to  the  transducer  layer  . 

VI.  FUTURE  WORK 

A.  ELECTRONICS  OF  THE  FABRICATED  ARRAY 

Currently,  transducer  testing  is  performed  by  manually  selecting  transmitting  and  receiving  array 
elements.  A  high  voltage  multiplexer  which  allows  electronic  selection  of  any  array  elements  has  been  purchased 
from  Matrix  Cooperation  and  the  expected  delivery  date  is  middle  October  this  year.  With  this  multiplexer,  we 
will  be  able  to  perform  acoustic  testing  with  the  fabricated  array  described  below. 

B.  ACOUSTIC  TESTING 

The  fabricated  array  with  back  layer  will  be  fully  tested  acoustically.  Testing  parameters  are  the  following. 

Impulse  re.sponse:  Similar  to  the  test  discussed  early. 

Acoustic  cross-talk:  This  measurement  will  be  made  by  operating  a  selected  individual  element  of  the 
array  in  transmit  mode  and  simultaneously  monitoring  the  output  signals  at  the  terminals  of  the  adjacent  elements 
of  the  array  operating  in  the  receive  mode. 

Pulse-echo  sensitivity.  Pulse-echo  sensitivity  will  be  measured  by  the  standard  method  of  bouncing  the 
launched  acoustic  pulse  off  a  perfect  reflector  and  monitoring  the  transmitted  and  received  signals. 

Single  element  directivity  pattern:  The  directivity  pattern  of  a  single  element  will  be  measured  in  pulsed 
transmit  mode.  A  needle  hydrophone  probe  will  be  placed  in  the  plane  of  the  far  field  of  the  selected  element  and 
the  lateral  beam  profile  of  the  element  will  be  measured. 

Radiation  pattern  of  entire  array:  The  2-way  radiation  pattern  of  the  fabricated  array  will  be  measured  in 
pulse-echo  mode  using  synthetic  aperture  technique.  The  center  transmitting  elements  will  provide  a  broad 
transmitting  beam.  As  a  result,  the  measured  2-way  radiation  pattern  is  the  1-way  radiation  pattern  of  the 
receiving  aperture.  Main  lobe  and  sidelobe  characteristics  can  be  evaluated  from  the  measured  radiation  pattern. 

C.  ARRAY  DEVELOPMENT 

After  the  completion  of  acoustic  testing  of  the  IMHZ  array,  we  will  have  complete  knowledge  about  array 
fabrication  procedures  and  pitfalls.  Our  next  milestone  is  the  development  of  a  3.75  MHz,  0.75X.  array  with  an 
element  size  of  0.3  mm.  This  frequency  is  suitable  for  breast  imaging. 

PARTC.  COUPLING  ULTRASOUND  AND  DIFFUSED  LIGHT  IMAGING  METHODS 

I.  INTRODUCTION 

Standard  ultrasound  is  excellent  in  differentiating  cysts  from  solid  lesions  and  is  routinely  used  in 
conjunction  with  x-ray  mammography.  Ultrasound  can  resolve  objects  a  few  mm  in  size  in  a  homogeneous 
medium,  and  its  resolution  in  the  inhomogeneous  breast  is  steadily  improving  as  described  early.  However, 
ultrasound  sensitivity  and  specificity  are  low  which  may  be  due  to  the  fact  that  even  when  tissues  are 
pathologically  different,  their  ultrasound  properties  are  often  quite  similar  [Cl].  Near  infrared  (NIR)  diffusive 
light  is  expected  to  have  high  sensitivity  and  specificity  because  it  affords  four  novel  biochemistry  parameters 
(blood  volume,  blood  02  saturation,  mitochondrial  light  scattering  and  levels  of  contrast  agents)  which  differ  in 
tumors  and  normal  breast  tissues.  Three  of  the  four  parameters  are  easily  obtainable  noninvasively.  The  fourth 
requires  i.v.  delivery  of  optical  contrast  agents  which  may  achieve  a  significant  increase  of  specificity  and  afford  a 
positive  benefit/risk,  particularly  if  the  contrast  agents  indicate  not  only  the  extent  of  the  tumor  vascular  bed  and 
extravasation.  Further,  NIR  may  eventually  indicate  more  fundamental  biological  characteristics  of  tumors  such 
as  mitochondrial  membrane  potential,  intracellular  calcium  levels,  and  plasma  membrane  potentials,  as  a  potential 
result  of  researches  in  this  rapidly  developing  field  [C2].  However,  NIR  diffused  light  suffers  poor  spatial 
resolution  because  of  the  diffuse  nature  of  NIR  light  in  tissue.  Currently,  typical  reconstructions  can  distinguish 
simple  structures  of  approximately  1  cm  in  size,  sharp  edges  are  typically  blurred  by  a  few  millimeters  [C2].  The 
novel  combination  of  ultrasound  and  NIR  light  methods  reported  here  has  the  potential  to  improve  detection 
sensitivity  and  specificity  of  either  single  method.  Combination  took  advantage  of  ultrasound  imaging  capability 
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’  and  MR  high  sensitivity/specificity  potentials  to  increase  detection  sensitivity  and  specificity  of  small  targets. 
MR  data  provided  a  Region  of  Interst  (ROI)  for  human  observers  to  search  for  possible  targets  and  lesions  in 
ultrasound  images.  As  a  result,  detection  sensitivity  and  specificity  has  been  improved  significantly  in  phantom 
studies.  Further  clinical  investigations  are  needed.  If  successful,  this  could  potentially  result  in  a  low  cost 
screening  device  which  could  assist  x-ray  mammography  for  early  detection  of  breast  cancers  and  reduction  of  the 
number  of  normal  breast  biopsies,  and  could  have  a  highly  beneficial  impact  on  the  early  treatment  of  breast 
tumor. 

A.  OPTICAL  SPECTROSCOPY  AND  IMAGING 

Over  the  last  two  decades,  the  scope  and  importance  of  research  being  conducted  on  optical  methods  for 
spectroscopy  and  imaging  has  increased  dramatically  [C2].  Spectroscopy  is  useful  for  the  measurement  of  time- 
dependent  variations  in  the  absorption  and  scattering  of  large  tissue  volumes.  For  example,  brain  oxymetry 
(hemoglobin  spectroscopy)  can  reveal  internal  bleeding  caused  by  head  injury.  Imaging  is  critical  when 
localization  of  heterogeneities  inside  tissues  is  involved,  such  as,  an  early  detection  of  breast  or  brain  tumors. 
Optical  spectroscopy  and  imaging  of  body  structure  and  function  is  made  possible  by  a  spectrum  window  that 
exists  within  tissues  in  the  700  -  900  nm  region  (near  infrared),  in  which  photon  transport  is  dominated  by 
scattering  rather  than  absorption  (Fig.cl).  Thus,  to  a  very  good  approximation,  NIR  photons  diffuse  through 
relatively  thick  tissues,  such  as  10-20  cm  of  a  large  human  breast. 

Compared  with  other  imaging  methods,  the  optical  method  has  numerous  advantages.  NIR  light  is 
expected  to  have  high  sensitivity  and  specificity  because  certain  characteristics  differ  in  tumors  and  normal  breast 
tissue:  1)  The  amount  of  blood  needed  to  serve  the  tumor's  metabolic  needs  is  increased  over  that  of  normal 
background  tissue;  2)  Hemoglobin  desaturation  in  tumors  is  increased  due  to  the  high  oxygen  demand  of 
cancers;  3)  Light  scattering  properties  of  the  tumor  are  expected  to  be  enhanced  due  to  the  increased 
mitochondria  population  with  respect  to  the  background  of  normal  cells  [C3];  4)  The  angiogenesis  phenomenon 
[C4]  allows  tumor  identification  by  observing  the  delivery  of  contrast  agents,  such  as  the  well-characterized  NIR 
absorber  indocyanine  green  (ICG),  through  the  permeable  vascular  beds  to  the  extravascular  space  surrounding 
the  tumor.  The  intense  absorption  band  of  ICG  at  800  nm  affords  a  high  contrast  of  the  tumor  extravascular 
space  as  compared  with  the  background  capillary  network  of  the  normal  tissue.  In  addition,  ICG  has  an 
interesting  redshifted  fluorescence  property  [C5].  It  absorbs  radiation  at  the  source  wavelength  of  780  nm  and 
reradiates  photons  at  830  nm.  In  this  process,  the  tumor  extravascular  space  is  converted  into  a  source  of 
secondary  diffuse  waves  with  reduced  background.  Finally,  the  contrast  agents  could  be  linked  to  tumor-specific 
antibodies  or  to  specific  sites  in  the  tumor  such  as  their  mitochondria  membranes.  In  addition,  the  NIR  device 
requires  relatively  low  manufacturing  and  operating  costs,  NIR  light  is  non-invasive  to  the  human  body  and  the 
device  is  compact  and  easily  portable. 
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Fig.cl.  In  vitro  absorption  spectra  of  hemoglobin  and  water,  showing  a  spectral  window  in  tissue  in  the  NIR.  The 
window  is  brought  about  by  first  a  decrease  in  blood  (oxy-  and  deoxyhemoglobin)  absorption  and  then  an  increase  in  water 
absorption  with  increasing  wavelength  (from  [C2]). 
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■  However,  there  are  intrinsic  difficulties  in  obtaining  images  with  high  spatial  resolution  because  of  the 
diffuse  nature  of  light  propagation  in  human  tissue.  In  essentially  all  imaging  algorithms  are  based  on  computer 
tomographic  (CT)  reconstruction  techniques  [C8-C9,C38-C44].  Tomographic  approaches  require  one  to  invert 
the  diffusion  equation  or  to  use  iterative  techniques  and  the  problem  is  more  complicated  than  the  x-ray  CT 
approach  because  photons  do  not  travel  along  straight  lines.  Currently,  typical  reconstructions  can  distinguish 
simple  structures  of  approximately  1  cm  in  size;  sharp  edges  are  typically  blurred  by  a  few  millimeters  [C2]. 

B.  OPTICAL  PHASED  ARRAY: 

A  different  approach  to  localization  for  NIR  detection  is  given  in  [CIO].  In  Ref.  [CIO],  a  simple  optical 
phased  array  system  (OP AS)  for  object  positioning  has  been  constructed  that  measures  the  interference  of 
intensity-niodulated  NIR  diffuse  light  signals  (in-phase,  0^  and  out-of-phase,  180°)  at  detector  positions 
equidistant  from  two  spatially  separated  sources.  In  the  absence  of  heterogeneities  (Fig.c2a)  the  signals  will 
cancel,  and  the  null  plane  or  line  is  always  at  the  bisecting  plane  or  line  and  the  deviation  between  null  and 
bisecting  lines  is  zero.  If  a  hidden  heterogeneity  is  present  in  the  medium  (Fig.c2b),  the  balance  of  the  sources  on 
the  bisecting  line  is  perturbed  and  the  null  line  will  deviate.  This  deviation  is  used  as  an  indicator  in  the  scanner  to 
detect  an  inhomogeneous  object  with  respect  to  a  homogeneous  background.  The  line-of-sight  of  the 
heterogeneity  is  obtained  from  the  newly  formed  null  line  through  either  mechanical  or  electronic  scanning.  This 
type  of  "pointing"  device  affords  highly  sensitive  localization  of  a  very  small  (~1  mm)  single  objects  to  an  angular 
accuracy  of  a  few  degrees  [C22].  In  vivo  experiments  have  demonstrated  clear-cut  amplitude  nulls  and  sharp 
phase  transitions  [C10,C1 1].  Simple  imaging,  referred  to  as  optical  phased  array  imaging  or  null  imaging  in  the 
following  text,  can  be  made  by  displaying  spatial  derivatives  of  the  phase  and/or  amplitude  curves  obtained  from 
scanning  in  one  and  two  dimensions  [C12].  A  larger  derivative  may  indicate  the  existence  of  a  larger  mass  while 
a  smaller  derivative  may  indicate  a  smaller  mass.  However,  exact  object  shape  and  size  information  and  optical 
properties  (absorption,  scattering  coefficients,  etc.)  cannot  be  obtained  from  the  null  imaging. 

INTENSITY  MODULATED  LIGHT 


C.  PULSE-ECHO  ULTRASOUND  IMAGING 

The  concept  of  pulse-echo  ultrasound  imaging  is  well  known  [C13].  A  short  burst  of  acoustic  pulses  is 
generated  by  a  piezoelectrical  transducer  array,  which  is  capable  of  converting  electrical  energy  to  acoustic  energy 
and  vice  versa,  and  is  directed  into  the  body  along  any  desired  scan  direction;  echoes  are  generated  and  travel  back 
to  the  transducer  for  detection.  These  echoes  vary  in  intensity  according  to  the  type  and  structure  of  tissue.  A  B- 
scan  image  is  formed  by  sweeping  a  narrow  acoustic  beam  through  a  plane  and  producing  display  of  the  received 
echoes,  such  that  there  is  a  correspondence  between  the  display  scan  line  and  the  direction  of  acoustic  propagation 
in  the  tissue.  Ultrasound  has  many  advantages;  its  manufacturing  and  operating  costs  are  low;  it  is  non-invasive 
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’  to  the  'human  body  and  it  provides  real-time  images.  The  use  of  ultrasound  to  examine  the  breast  was  first 
described  in  [14].  Now,  ultrasound  is  frequently  used  in  conjunction  with  x-ray  mammography  to  differentiate 
cysts  from  solid  lesions.  However,  ultrasound  applications  with  respect  to  early  detection  of  small  solid  lesions 
were  disappointing  [Cl 5].  Ultrasound  sensitivity  in  the  detection  of  palpable  solid  masses  less  than  2  cm  in  size 
was  71%,  which  was  about  10%  less  than  the  sensitivity  of  physical  examination  [C16].  The  sensitivity  in  the 
detection  of  small  invasive  lobular  carcinomas  (<1  cm)  was  only  25%  [C17].  Ultrasound  specificity  was  low  due 
to  the  considerable  overlap  of  acoustic  characteristics  between  benign  and  malignant  lesions  [C18-C19].  More 
recently,  the  imaging  capability  of  ultrasound  equipment  has  been  improved.  For  example,  the  results  reported  by 
Stavros  et  al.  were  very  encouraging  [C20].  Although  no  lesion  size  was  documented  in  [C20],  of  125  malignant 
lesions,  123  were  correctly  classified  as  intermediate  or  malignant  (98.4%  sensitivity). 

II.  EXPERIMENTS  AND  RESULTS 

A  novel  combination  of  the  two  suboptimal  systems,  namely  ultrasound  and  OP  AS,  reported  in  [C21] 
and  the  Appendix  C  has  shown  a  superior  performance  in  detection  sensitivity  and  specificity  compared  with 
ultrasound  alone  in  phantom  studies.  The  combination  took  advantage  of  ultrasound  imaging  and  OPAS 
localization  capabilities,  and  directed  human  observers  to  look  for  possible  targets  and  lesions  within  a  ROI  in  the 
ultrasound  image.  The  ROI  was  determined  from  highly  sensitive  optical  line-of-sight  information.  As  a  result, 
the  detection  sensitivity  and  specificity  of  a  few  mm  object  buried  in  a  scattering  phantom  has  been  improved 
significantly.  Experimental  set-ups,  examples  and  statistics  are  described  in  the  following  paragraphs. 

In  our  phantom  studies  [C21],  we  incorporated  the  simple  NIR  positioner  into  an  ultrasound  imaging 
system  (Fig.c3).  An  ultrasound  phased  array  probe  is  coaxially  mounted  with  the  two  laser  diode  sources, 
referred  to  as  the  coaxial  probe  in  the  figure,  and  the  coaxial  probe  and  the  optical  detector  assembly  is 
mechanically  translated  together  on  a  rigid  frame  by  a  stepping  motor.  This  simple  combined  system  was  tested 
using  models  consisting  of  a  single  emulated  tumor  of  a  few  mm  size  buried  in  a  scattering  phantom.  The 
homogeneous  scattering  phantom  was  made  of  an  acrylamide  gel  (7.5  x  12  x  14  cm^),  and  scattering  particles  at 
optical  and  ultrasound  wavelength  scales  were  added  to  the  acrylamide  solution  before  the  gel  was  formed.  In 
each  test,  the  combined  system  was  mechanically  translated  to  scan  for  the  small  object.  In  a  typical  experiment, 
1  mm  step  size  was  used.  At  each  scanning  position,  optical  amplitude  and  phase  data  and  an  ultrasound  sector 
scan  image  were  acquired.  The  object  line-of-sight  information  was  obtained  from  the  optical  data  and  was  coded 
into  ultrasound  images  after  scanning.  A  total  of  30  targets  with  various  acoustic  and  optical  absorption  and 
scattering  properties  were  studied,  and  ultrasound  images  without  and  with  line-of-sight  information  were  graded 
on  a  five-point  scale  (see  Table  I)  by  two  experienced  radiologists  in  the  Breast  Imaging  Section  of  the  Radiology 
Department.  A  typical  example  of  optical  scanning  and  ultrasound  image  evaluation  is  demonstrated  as  follows. 
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Fig.c3.  (a)  Coaxial  ultrasound  and  optical  probe  and  optical  detector  assembly,  (b)  front  view  of  the  coaxial  probe. 


Fig.c4  shows  optical  amplitude  and  phase  responses  when  concentric  black-plastic  tubes  of  1.9  and  1.7 
mm  diameters  filled  with  0.5%  intralipid  (optical  scattering  medium)  were  inserted  into  the  scattering  phantom  at 
approximately  6  cm  depth.  The  testing  object  is  an  optical  absorber  which  also  generates  an  inhomogeneous 
ultrasound  echotexture  (see  Fig.c5)  as  some  tumors  may  do.  The  sharp  phase  transition  and  amplitude  null 
correspond  to  the  object  line-of-sight.  As  a  result  of  mechanical  translation,  a  sequence  of  ultrasound  images  was 
obtained.  Since  weak  targets  in  a  sequence  of  images  were  more  readily  detected  than  in  a  single  image,  triple 
images  taken  in  1  mm  steps  at  the  neighborhood  of  the  optical  phase  transition  and  amplitude  null  were  selected  in 
each  test.  The  criteria  of  triple  image  selection  are  described  in  detail  in  Appendix  A.  Since  the  optical  assembly 
and  ultrasound  probe  are  coaxially  mounted,  the  object  in  the  triple  images  always  appears  along  the  ultrasound 
scan  axis.  To  avoid  giving  any  prior  knowledge  to  human  observers  about  the  object  location,  we  chose  different 
image  windows  in  different  tests  to  randomize  the  object  location.  In  each  test,  the  image  windows  used  for  triple 
images  were  the  same  and  the  object  position  shift  was  caused  by  the  1  mm  lateral  translation  of  the  imaging 
system.  In  this  example,  three  images  taken  at  phase  transition  and  amplitude  null  positions  (16,  17,  18  mm)  are 
shown  in  Fig.c5.  The  image  of  the  object,  indicated  by  an  arrow,  is  highly  obscured  by  speckle  noise.  Two 
observers  evaluated  the  triple  image  using  the  form  given  in  Table  I.  They  were  told  that  they  could  mentally 
average  these  images  to  filter  out  artifacts  and  noise.  One  observer  picked  up  the  target  and  marked  category  4  on 
the  evaluation  form.  When  the  same  images  coded  with  optical  line-of-sight  information,  shown  by  dashed  lines 
in  Fig.c6,  were  presented  to  the  observer,  his  detection  confidence  increased  from  category  4  to  5.  Another 
observer  completely  missed  the  target  in  ultrasound-alone  images  and  marked  category  1  on  the  form.  When  the 
images  with  object  line-of  sight  information  were  presented,  the  second  observer  picked  up  the  target  and  selected 
category  3  on  the  form.  This  example  shows  the  improvement  in  detection  sensitivity  when  the  combined 
method  is  used. 


Table  I.  Evaluation  form 
Please  check  for  possible  lesions  or  targets. 

Target  #1 

1 .  No  lesion  or  target  present 

2.  Probably  no  lesion  or  target  present 

3.  Lesion  or  target  may  be  present 

4.  Lesion  or  target  probably  present 

5.  Lesion  or  target  definitely  present 

Target  #2 

1.  No  lesion  or  target  present 

2.  Probably  no  lesion  or  target  present 

3.  Lesion  or  target  may  be  present 

4.  Lesion  or  target  probably  present 

5.  Lesion  or  target  definitely  present 

Target  #3... 

Target  #4...  ^  ^  i  •  i  ^ 

Lesions,  or  targets  may  be  hyperechogenic  or  hypoechogenic.  If  you  identify  any  lesions  please  circle  and 

number  them  on  images.  _  . 
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(a)  (b 

Fig.c4.  Optical  response  when  concentric  black  plastic  tubes  of  1.9  and  1.7  mm  diameters  filled  with  0.5%  intralipid 
were  inserted  into  the  scattering  phantom,  (a)  Amplitude  response  and  (b)  phase  response.  The  amplitude  null  and 
phase  transition  are  indications  of  the  object  line-of-sight.  _ 
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Fig.cS.  The  triple  images  taken  at  lateral  positions  16, 17, 18  mm.  Target  is  indicated  by  arrow. 
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Fig.c6.  The  same  triple  images  as  shown  in  Fig.c5  but  optical  line-of-sight  information  was  coded  in  the  images  as 
indicated  by  dashed  lines.  Target  is  indicated  by  arrow. 

Figs  c7-cl9  show  another  example.  Fig.c7  shows  optical  amplitude  and  phase  responses  when  a  4.5  mm 
white  Teflon  bar  was  inserted  into  a  scattering  phantom  at  approximately  6  cm  depth  and  was  imaged  using  the 
combined  system.  The  testing  object  is  an  optical  scatterer  which  generates  a  slightly  hyperechoic  echotexture. 
The  triple  images  taken  at  the  neighborhood  of  optical  phase  transition  and  amplitude  null  at  1mm  steps  (20,  21, 
22  mm)  are  shown  in  Fig.cS  with  the  object  indicated  by  the  upper  arrow  in  the  images.  The  image  of  the  object 
is  highly  obscured  by  speckle  noise.  One  observer  missed  the  target  entirely  and  marked  category  1  on  the 
evaluation  form.  When  the  images  with  object  line-of-sight  (Fig.c9)  were  presented,  the  observer  identified  the 
target  and  circled  category  3  on  the  form.  Another  observer  missed  the  target  completely  and  selected  an  artifact, 
pointed  by  the  lower  arrow,  and  circled  category  4  on  the  form.  Once  the  images  with  object  line-of-sight  were 
presented,  the  observer  identified  the  target  and  circled  category  3  on  the  form.  This  example  shows  the 
combined  method  can  significantly  improve  detection  sensitivity  and  specificity  of  ultrasound. 


Fig.cV.  Optical  response  when  a  4.8  mm  white  Teflon  rod  was  inserted  into  the  scattering  phantom,  (a)  Amplitude 
response  and  (b)  phase  response.  The  amplitude  null  and  phase  transition  are  indications  of  the  object  line-of-sight. 
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Fig.c9.  The  same  triple  images  as  shown  in  Fig,c8  but  optical  line-of-sight  information  was  coded  in  the  images  as 
indicated  by  dashed  lines.  Target  and  artifact  are  indicated  by  arrows. 

We  have  studied  a  total  of  30  positive  trials  and  21  negative  trials.  In  the  positive  trials,  the  testing  objects 
used  were  optical  absorbers  and  scatterers  which  emulated  to  a  large  extent  tumor  optical  properties  and  generated 
various  acoustic  echo  textures.  The  sizes  of  the  testing  objects  range  from  1.6  mm  to  8.2  mm  with  an  average  of 
3.64  mm  (std.  1.97  mm).  In  the  negative  trials,  images  of  the  scattering  phantom  obtained  at  three  consecutive 
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'lateral  positions  in  1  mm  steps  were  presented.  Receiver-Operator  Characteristic  (ROC)  curves,  commonly  used 
to  depict  figures  of  merit  for  clinical  methods  of  detecting  lesions  or  differentiating  among  diseases  [C22,C23], 
are  measured  from  the  readings  of  the  two  radiologists.  Average  empirical  ROC  curves  of  ultrasound  (solid)  and 
ultrasound-light  (dashed)  of  two  independent  readings  from  the  two  observers  are  shown  in  Fig.clO.  The 
superior  performance  of  ultrasound-light  is  shown  by  the  higher  ROC  curves.  The  area  under  the  ROC  curve 
provides  a  measure  of  discriminatory  efficacy:  an  area  of  100%  denotes  an  ideal  system  and  50%  denotes  an 
ineffective  system.  The  calculated  areas  under  ROC  curves  of  ultrasound  and  ultrasound-light  of  observer  #1  are 
53.9%  and  79.2%,  and  of  observer  #2  are  49.5%  and  85.8%,  respectively.  The  mean  values  are  51.7%  for 
ultrasound  and  82.5%  for  ultrasound-light.  The  difference  in  area  between  the  two  methods  is  statistically 
significant  (p<0.05,  see  Appendix  of  [C21]). 


Fig.clO.  Average  ROC  curves  of  two  readings  from  two  observers.  The  ROC  curves  of  ultrasound  method  are  plotted 
in  solid  lines  and  ROC  curves  of  ultrasound-light  are  plotted  in  dashed  lines. 

III.  CONCLUSION 

A  novel  imaging  method  which  combines  ultrasound  and  near  infrared  diffuse  light  has  been  introduced 
in  this  paper.  The  combined  method  utilizes  object  line-of-sight  information  provided  by  a  highly  sensitive  optical 
positioner  to  confine  human  observers  along  the  object  direction  to  search  for  possible  lesions  or  targets  in 
ultrasound  images.  As  a  result,  the  detection  sensitivity  and  specificity  of  a  small  target  of  a  few  mm  size  have 
been  improved  significantly  over  the  conventional  ultrasound.  The  average  areas  under  ROC  curves  from  2 
radiologists  are  much  higher  for  the  combined  method  (82.5%)  than  for  ultrasound  alone  (51.7%).  Thus  the 
combined  method  affords  a  significant  improvement  in  detection  sensitivity  and  specificity  compared  with 
conventional  ultrasound.  Further  clinical  investigations  are  required  to  evaluate  this  new  method  in  breast  cancer 
detection.  Currently,  we  are  applying  funds  to  explore  the  utility  of  this  new,  exciting  method  through  clinical 
trials. 
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ABSTRACT:  There  are  two  types  of  wavefront  distortion  inside  the 
female  breast:  incoherent  scattering  and  coherent  Interference.  Adap¬ 
tive  wavefront  compensation  algorithms  developed  so  far  are  useful 
for  correction  of  incoherent  distortion  caused  by  scattering.  The  per¬ 
formances  of  these  algorithms  differ  according  to  the  extent  of  wave- 
front  amplitude  distortion.  It  is  shown  analytically  and  experimentally 
in  this  article  that  the  matched  filtering  approach,  while  optimal  for 
detection,  is  not  so  for  imaging  when  the  wavefront  amplitude  is 
distorted  in  addition  to  phasefront  distortion.  Matched  filtering  algo¬ 
rithms  increase  wavefront  amplitude  variance  and  therefore  decrease 
image  contrast.  The  inverse  filtering  approach,  while  ideal  for  fidelity, 
is  not  stable  when  the  signal -to-noise  ratio  is  low.  An  approach  toward 
inverse  filtering,  amplitude  compression  operation  in  addition  to 
phase  deaberration.  Is  introduced  in  this  article.  Analysis  and  experi¬ 
ments  show  that  its  performance  is  superior  to  matched  filtering  algo¬ 
rithms  and  to  time-delay  type  correction  algorithms  without  introduc¬ 
ing  stability  problems.  ©  1996  John  Wiley  &  Sons,  Inc.  Int  J  Imaging  Syst 
Technol,  8,  322-335,  1997 

Key  words:  ultrasound;  wavefront  compensation;  inverse  filtering; 
wavefront  model 


I.  INTRODUCTION 

In  principle,  high  ultrasound  resolution  can  be  achieved  with 
very  large  ultrasound  transducers.  However,  large  transducers  are 
more  susceptible  to  wavefront  distortion  caused  by  inhomoge¬ 
neous  tissue,  because  the  transducers  cover  large  volumes  of 
tissue  and  the  beams  interact  with  more  spatial  variations  in  tissue 
composition.  The  distortion  inside  the  female  breast  falls  into 
two  categories:  incoherent  scattering  reduces  the  target  strength, 
broadens  the  image  lobe  [1,2],  and  raises  the  background  level, 
and  therefore  lowers  the  image  contrast  [3,4].  Refraction  creates 
coherent  multipath  interference  that  produces  false  targets  or  im¬ 
age  artifacts  in  addition  to  true  targets  in  the  image  [4,10].  The 
interference  problem  is  more  severe  when  the  aperture  is  large. 

Numerous  efforts  have  been  made  to  correct  wavefront  distor- 
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tion  produced  by  breast,  liver,  and  abdominal  wall  [9,13-17]. 
Basic  deaberration  algorithms  [5-7]  are  recognized  to  be  suitable 
only  for  weak  scattering  that  can  be  modeled  as  a  thin  random- 
phase  screen  located  in  the  plane  of  the  receiving  aperture  [8- 
9,18].  Extensions  of  the  basic  algorithms  are  made  to  correct 
amplitude  distortion  caused  by  strong  scattering  or  distributed 
scattering  [8,9] .  A  model-based  approach  that  uses  a  priori  infor¬ 
mation  of  the  speeds  of  the  rectus  muscle  layers  inside  the  abdom¬ 
inal  wall  has  been  developed  to  correct  double  image  artifacts 
caused  by  refraction  [17].  The  proposed  technique  is  not  model 
based  and  therefore  is  freed  from  the  necessity  of  having  a  priori 
information  about  the  source  of  the  distortion. 

The  motivation  of  this  study  stems  from  the  following.  The 
algorithms  [5-9]  that  require  no  a  priori  knowledge  about  the 
propagating  medium  function  differently  upon  incoherent  scatter¬ 
ing  and  coherent  interference;  therefore,  it  is  important  to  under¬ 
stand  the  effects  of  the  two  distortion  processes  upon  the  image 
and  then  apply  suitable  correction  mechanisms  to  the  algorithms. 
The  scattered  energy  increases  with  the  propagation  depth  in  the 
breast;  therefore,  the  performance  of  descattering  algorithms  [5- 
9]  is  expected  to  deteriorate  with  depth.  Modeling  and  quantifi¬ 
cation  of  depth-dependent  distortion  and  correction  are  extremely 
valuable  for  the  design  of  next-generation  high-quality  breast 
scanners. 

In  this  article  we  analyze,  through  modeling  and  experiment, 
image  distortion  caused  by  scattering  and  coherent  interference, 
and  the  extent  to  which  image  quality  improvements  are  affected 
by  deaberration.  Wavefront  compensation  algorithms  based  on 
matched  filtering  are  first  recognized  to  be  unsuitable  for  imaging 
in  distributed  scattering  because  the  algorithms  increase  wave- 
front  amplitude  distortion  and  therefore  decrease  image  contrast. 
A  compression  operation  on  wavefront  amplitude,  akin  to  inverse 
filtering,  which  can  significantly  reduce  wavefront  amplitude 
variance  is  introduced.  We  demonstrate  through  experiments  that 
a  combination  of  the  amplitude  compression  and  either  phase 
deaberration  at  the  aperture  or  backpropagation  plus  phase  deab¬ 
erration  at  some  optimal  distance  from  the  aperture  can  sig¬ 
nificantly  improve  image  contrast  and  restore  the  mainlobe 
diffraction  shape  to  a  -30-dB  level.  Linear  regression  analysis 
of  side-lobe  energy  versus  propagation  depth  from  available  data 
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Figure  1 .  Approximate  behavior  of  the  coherent  intensity  4  and  the 
incoherent  intensity  // . 


Figure  2.  Angular  source  intensity  distribution.  Statistical  model,  (a) 
Before  phase  deaberration,  showing  reduced  coherent  spectrum  and 
of  incoherent  spectrum,  (b)  After  phase  deaberration. 


indicates  that  good  focusing  with  high-contrast  resolution 
throughout  the  100-mm  propagation  depth  by  deaberration  is 
possible. 

II.  SCATTERING  MODELS 

The  wave  propagation  problem  is  formulated  similarly  to  the 
line-of-sight  propagation  problem  through  a  tenuous  distribution 
of  particles  [19].  A  narrowband  wave  is  used  to  simplify  the 
description  and  derivation.  The  analytic  results  can  be  taken  as 
rules  of  thumb  applicable  to  the  wideband  case  to  characterize 
contrast  resolution  due  to  wavefront  amplitude  and  phase  distor¬ 
tions. 

An  ultrasound  wave  (either  spherical  or  planar)  propagates  in 
an  inhomogeneous  medium  (Fig.  1).  It  consists  of  coherent  and 
incoherent  waves.  The  incoherent  wave  is  also  called  a  scattered 
wave.  As  the  wave  propagates,  the  coherent  intensity  is  attenu¬ 
ated  due  to  absorption  and  scatter,  and  the  incoherent  intensity 
li  is  scattered  broadly. 

A.  Weak  Scattering.  The  field  £  at  a  receiving  aperture  is 
predominantly  coherent  and  the  magnitude  of  the  incoherent  field 
is  much  smaller  than  that  of  the  coherent  field.  E  is  the  sum  of 
the  coherent  and  incoherent  fields. 

E{x)  =  E,{x)  +  Ei{x)  (1) 

EXx)  >  £/(Jc),  {EAx\)Ei{x2))  =  0 

and  {Ei{x))^0  (2) 

where  (  )  is  the  ensemble  average  and  x  is  the  spatial  variable 

at  the  receiving  aperture. 

The  spatial  correlation  or  mutual  coherence  function  of  the 
field  at  the  aperture  is  given  by  F  =  {E{xl)E{x2)) .  Since 
{Ei{x))  =  0  and  {Ec{xl)Ei{x2))  =  0,  the  spatial  correlation 
function  satisfies  F  =  F^  +  F/ .  If  the  receiver  is  located  in  the 
far  field  or  within  the  focal  region  of  the  source,  the  angular 
source  intensity  distribution  or  intensity  image  is  the  Fourier 
transform  of  the  spatial  correlation  function  and  equals  the  sum¬ 
mation  of  the  coherent  and  statistically  smoothed  symmetrical 
incoherent  spectra  [Fig.  2(a)].  For  each  realization,  however, 
the  incoherent  spectrum  can  have  valleys  and  hills.  The  Fourier 
relationship  between  the  incoherent  spatial  correlation  function 
F,  and  incoherent  spectrum  is  commonly  known  as  the  Van- 
Cittert  Zemike  theory  [20]. 

The  wavefront  distortion  is  primarily  in  the  phasefront  and 
can  be  modeled  approximately  as  due  to  a  thin,  random-phase 
screen  located  in  the  plane  of  the  receiving  aperture  and  defined 


by  a  spatial  correlation  distance.  Such  distortion  is  correctable 
by  phase-deaberration  or  time-delay  compensation  algorithms 
[5-7].  The  algorithms  operate  upon  the  incoherent  field.  They 
fold  the  scattered  energy  into  the  coherent  field  and  thereby  re¬ 
store  the  full  diffraction-limited  image  [Fig.  2(b)]. 

The  most  relevant  wavefront  measurement  in  a  weak  scatter¬ 
ing  field  is  the  spatial  correlation  function  of  the  phasefront  or 
time-delay  profile  as  reported  by  several  groups  [1,16,21,22]. 
From  this  information,  the  performance  of  the  phase-deaberration 
process  can  be  predicted. 

B.  Strong  Scattering  or  Distributed  Scattering.  The  coher¬ 
ent  component  is  attenuated  and  the  incoherent  component  is 
scattered  broadly.  Wavefront  amplitude  is  distorted  in  addition 
to  phase-front  distortion.  A  more  suitable  representation  of  the 
field  is  the  Rytov  model.  A  distributed  random  distortion  medium 
can  be  represented  as  a  cascade  of  thin,  independent  distortion 
layers  [Fig.  3(a)]  or  phase  screens  [19].  At  the  ith  layer,  the 
field  can  be  written  as  Ei  =  A/exp(7<^/).  Phase-front  distortion 
is  developed  owing  to  time  shift  through  each  layer.  Wavefront 
amplitude  distortion  is  developed  as  the  phase-aberrated  wave- 
front  propagates  through  the  medium  to  successive  layers  [Fig. 
3(b)] .  Let  Eq  =  Aoexp(7<^o)  be  the  reference  field.  The  distorted 
field  is  then  a  product  of  contributions  from  each  layer 

E(x)  =  EoExE2'  •  'Ei  -  '  'Em  (3) 


The  amplitude  of  E{x)!Eq{x)  denoted  as  At  is  At  = 
A1A2'  -  '  Ai'  •  '  An  and  the  phase  +  (/>2  +  *  *  *  <^)/  *  *  •  + 

(I)n.  Assume  that  the  layers  are  statistically  independent  and  have 
common  statistics;  therefore,  Ai  =  A  6Ai  and 
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Figure  3.  The  field  is  the  product  of  contributions  from  each  layer. 
Wavefront  amplitude  fluctuation  increases  as  wave  propagates. 
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Figure  4.  Statistical  models  of  angular  source  intensity  distribution, 
(a)  Before  phase  correction,  (b)  After  phase  correction. 


A,  =  (A)-n(i  +  f^) 

1  +  X  +  higher-order  products  j  (4) 

Since  A  is  close  to  unity  and  M,  /A,  i  =  1,  2,  •  •  -  are 
independent  variables,  the  amplitude  variance  of  A{  normalized 
to  (A)^  is  «  NasA/Ay  where  cr|^/^is  the  normalized  amplitude 
variance  of  a  single  screen,  grows  linearly  with  the  total 
number  of  phase  screens  N  and  therefore  with  the  propagation 
distance  D.  The  phase  variance  =  N(j\,  where  cr^  is  the 
variance  of  a  single  phase  screen;  therefore,  <7%^  grows  linearly 
with  N  and  D . 

The  overall  spectrum  (sum  of  coherent  and  incoherent  spec¬ 
tra)  is  broadened  [Fig.  4(a)]  and  reduced  in  strength.  The  back¬ 
ground  level  is  increased,  and  therefore  image  contrast  is  reduced. 

The  scattered  energy  increases  with  the  propagation  depth. 
The  energy  ratio  (ER),  which  is  the  energy  outside  the  main 
image  lobe  to  the  energy  inside  the  main  image  lobe,  is  a  useful 
tool  to  quantitatively  evaluate  the  depth-dependence  distortion 
upon  image  contrast.  The  relationship  between  ER  and  propaga¬ 
tion  depth  D  (see  Appendix  for  derivation)  is 

ER  =  +  aj^)  =  CN{al  +  gIaia) 


=  C  ^  (4  +  (5) 

a 

where  C  is  a  function  of  the  source  illumination  profile  and  d  is 
the  thickness  of  the  screen  and  is  approximately  the  correlation 
distance  of  the  medium.  ER  linearly  grows  with  the  propagation 
distance.  This  relationship  is  analogous  to  the  average  side-lobe 
floor  in  the  radiation  pattern  of  an  array  with  errors  [23]  and  to 
the  average  side-lobe  floor  in  an  arbitrary  source  distribution 
imaged  through  a  distorting  medium  [18].  Equation  (5)  can  be 
applied  to  wideband  as  an  approximation. 

So  long  as  certain  correlation  properties  exist  in  the  wavefront, 
phase  compensation  algorithms  are  still  useful  in  strong  scattering 
to  partially  remove  phase  distortion  and  build  up  the  strength  of 
the  coherent  field  [Fig.  4(b)].  Wavefront  amplitude  distortion 
remains,  as  does  residual  phase  distortion,  and  produces  a  sig¬ 
nificantly  high  incoherent  background  level  in  the  spectrum.  Any 
phase-deaberration  procedure  that  provides  a  better  phase  error 
estimate  than  that  of  the  basic  algorithms  [5-7]  can  reduce  this 
background  level.  Wavefront-deaberration  algorithms  that  take 
amplitude  distortion  into  account  can  further  reduce  this  back¬ 


ground  level  if  the  algorithms  reduce  wavefront  amplitude  vari¬ 
ance. 

Consider  the  distortion  vector  f  =  ElEo  ~  A^expOc^,).  Figure 
5  shows  such  a  distortion  vector  f  for  some  instantaneous  sample 
of  the  radiation  field.  The  optimum  compensation  weight  vector 
or  transformation  Wopt  =  A7^qx^{- is  the  rightmost  one 
which  carries  the  distortion  component  of  the  complex  sample 
to  the  intersection  of  the  x-axis  and  the  unit  circle,  for  then 
both  amplitude  and  phase  are  exactly  corrected.  The  leftmost 
transformation  w  =  A(Cxp{- j(j)c)  is  matched  filtering  (MF) 
which  squares  the  amplitude  and  conjugates  the  phase.  MF  is 
theoretically  optimum  for  maximizing  signal-to-noise  ratio 
(SNR)  on  a  white,  Gaussian  channel.  It  is  an  optimum  solution 
for  detection  but  not  when  fidelity  is  an  important  criterion,  as  in 
imaging.  But  it  does  satisfy  the  requirement  for  phase  correction. 
Examples  of  MF  algorithms  are  dominant  scatterer  algorithm 
(DSA)  [24]  and  time-reversal  mirror  (TRM)  [8].  We  should 
expect  that  MF  algorithms  are  useful  when  wavefront  distortion 
resides  principally  in  phase,  with  amplitude  relatively  unaffected. 
Otherwise  MF  algorithms  would  be  less  effective  than  if  ampli¬ 
tude  were  ignored  and  phase  only  were  corrected.  Later  we  show 
that  by  squaring  the  amplitude,  an  MF  algorithm  increases  the 
wavefront  amplitude  variance  and  thereby  increases  the  energy 
in  the  side-lobe  region. 

Consider  what  happens  when  a  phantom  consisting  of  a  ran¬ 
dom  phase  screen  is  placed  at  a  receiving  transducer.  The  element 
signal  samples  are  correct  in  amplitude  and  distorted  only  in 
phase  by  a  random  additive  component.  Phase  aberration  correc¬ 
tion  (PAC),  consisting  of  phase  conjugation  (PC),  is  then  the 
optimum  compensator.  The  signal  amplitude  is  unchanged.  This 
operator  is  represented  by  the  second  (from  the  left)  transforma¬ 
tion.  Time-delay  compensation  (TDC)  [5-7]  is  in  the  same  class. 

The  rightmost  transformation  which  corrects  both  amplitude 
and  phase  is  the  inverse  filter  (IF).  It  is  theoretically  ideal  for 
fidelity  but  has  drawbacks.  Because  the  IF  adjusts  the  channel 
gain  to  be  the  reciprocal  of  the  signal  strength,  at  points  in  the 
receiving  aperture  where  signal  strength  is  weak  the  enhanced 
channel  gain  raises  the  noise  to  the  point  where  SNR  can  be 
impaired.  This  is  particularly  troublesome  when  there  is  coherent 
refractive  interference  in  the  receiving  array.  A  second  potential 


DEABERRATION  TRANSFORMATIONS 


MF  =  matched  filter  (SNR) 

DSA  =  dominant  scatterer  algorithm 
TRM  =  time  reversal  mirrior 

PC  =  phase  conjugation 
TDC  -time-delay  com^sation 
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following  phase  correction 


Figure  5.  Deaberration  transformations,  f  is  distortion  vector  for 
instantaneous  sample  of  the  radiation  field.  Transformations  to  the 
real  axis  all  correct  phase.  Left  (MF)  maximizes  SNR  but  increases 
amplitude  distortion.  Right  (IF)  maximizes  imaging  fidelity  but  in¬ 
creases  noise.  SR  is  a  good  compromise. 
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problem,  which  we  have  avoided  by  not  working  to  close  to  the 
IF  point,  is  that  an  IF  is  unstable  when  the  distorting  medium  has 
zeroes  in  the  complex  plane.  Thus,  the  rightmost  transformation 
carries  the  correction  too  far. 

The  complex  weights  of  MF,  PC,  and  IF  are,  respectively,  A^, 
1,  and  each  multiplied  by  exp(-y<^;).  The  signal  after 
weighting  is  A? ,  A, ,  and  1 .  The  first  and  last  are  far  from  optimum 
for  the  reasons  given.  The  third  from  the  left  transformation  is  a 
compromise.  The  correction  vector  is  Ar^^^exp(-7<^f).  The  sig¬ 
nal  after  weighting  is  A^^  (called  later  the  square  rooter).  An¬ 
other  transformation  to  the  right  of  the  square  rooter  is  fourth 
rooter  (signal  amplitude  is  A^"^  after  weighting).  In  general,  the 
correction  vector  is  Ar^^~‘^^^exp(-y<^,)  and  signal  amplitude 
is  after  transformation  of  the  Mth  rooter.  In  the  next  section, 
we  demonstrate  experimentally  that  low-order  rooters  (square 
rooter  and  fourth  rooter),  in  addition  to  phase  deaberration  at  the 
aperture,  can  significantly  improve  contrast  resolution  without 
encountering  the  stability  problem. 

The  DSA  developed  in  narrowband  radar  [24]  is  designed  for 
adaptively  focusing  a  large  distorted  receiving  array  upon  a  point 
target  and  is  a  matched  filter  process  in  space.  DSA  creates  a 
compensating  weight  vector  w  from  the  complex  conjugate  of 
the  received  echo,  i.e.,  w  =  A,exp(”7(^;).  The  DSA-corrected 
receiving  beam  is  the  Fourier  transform  of  wAt&xp{j(j)f)  =  A?. 
Consequently,  DSA  uses  wavefront  amplitude  to  provide  addi¬ 
tional  SNR  gain  than  phase  conjugation  only,  for  which  w  = 
exp(  —  y</)^) .  However,  when  the  wavefront  amplitude  is  distorted, 
the  amplitude  weight  further  increases  the  side-lobe  level  in  the 
beam.  The  measure  in  the  increase  is  the  increased  wavefront 
amplitude  variance.  Since  A?  =  A^A^*  •  •  Af  •  •  •  A%\  therefore. 
Equation  (4)  is  changed  to 

A?=(A)-fi(i+f^y 

=  +  X  2  -A  +  higher-order  products^  (6) 

The  variance  of  square  amplitude  ^  4^0-^/^,  which  is  four 
times  larger  than  without  amplitude  weight.  Assuming  phase  er¬ 
ror  is  compensated,  ER  will  rise  from  C(Dld)a8AiA  to 
C {D/d)4alA/A,  which  is  about  a  6-dB  decrease  in  image  contrast. 

The  TRM  developed  by  Fink  [8]  provides  optimal  focusing 
of  a  wideband  transducer  upon  a  target  through  a  scattering  me¬ 
dium  and  is  a  matched  filter  process  in  both  space  and  time 
(waveform).  The  operation  reverses  the  received  waveforms,  re¬ 
transmits  them  through  the  medium,  and  then  focuses  the  array 
on  the  target.  The  transmitting  beam  pattern  is  the  convolution 
of  the  received  waveforms  and  their  time-reversed  version.  In  the 
frequency  domain,  TRM  is  a  phase  conjugation  procedure  (with 
amplitude  weight)  at  each  frequency  within  the  signal  band.  It 
is  anticipated  that  TRM  focusing  produces  poorer  illumination 
in  the  neighborhood  of  the  target  than  does  phase  conjugation 
focusing,  for  the  reason  discussed  in  the  previous  paragraph. 

Square  root  amplitude  weighting  changes  Equation  (4)  to 


=  (A)^^^^l  -h  X  +  higher-order  products^  (7) 


input  amplitude 


Figure  6.  Plots  of  the  input  and  output  amplitude  transforms. 


The  amplitude  variance  «  Nct^aiaI^^  which  is  four  times 
less  than  with  unit  amplitude  weight.  Assuming  phase  error 
is  compensated.  ER  will  reduce  from  C{Dld')a\AiA  to 
C(D/4£/)(Jm/4,  which  implies  about  a  6-dB  increase  in  image 
contrast.  Mth  rooter  will  reduce  the  amplitude  variance  by  2  M 
if  the  SNR  is  still  adequate.  This  nonlinear  transformation  upon 
wavefront  amplitude  is  a  compression  operation  (Fig.  6).  The 
detailed  evaluation  of  the  effects  of  this  compression  operation 
upon  image  can  be  performed  by  studying  the  amplitude  distribu¬ 
tions  of  images  before  and  after  the  operation.  This  is  a  subject 
of  a  later  study. 

Backpropagating  the  received  wavefront  to  an  optimal  back- 
propagation  distance  and  performing  phase  deaberration  at  this 
optimal  distance  (BPT)  [9]  provides  first-order  correction  of 
wavefront  amplitude  distortion  due  to  propagation  from  the  back- 
propagation  distance  to  the  receiving  aperture  and  therefore  per¬ 
forms  better  than  phase  deaberration  at  the  aperture.  In  the  next 
section,  we  demonstrate  experimentally  that  the  BPT  transforma¬ 
tion  is  somewhere  to  the  right  of  TDC  and  PC  and  left  of  PC 
and  rooter.  We  also  demonstrate  that  square  rooter  in  addition  to 
BPT  can  significantly  improve  the  performance  of  BPT. 

The  relevant  measurement  in  a  strong  scattering  region  is 
wavefront  amplitude  variance  [15]  (narrowband)  or  wavefront 
energy  fluctuation  [16,21]  (wideband)  in  addition  to  arrival  time 
variance.  Since  wideband  waveforms  can  be  Fourier  decomposed 
into  a  set  of  narrowband  complex  wavefronts  within  a  signal 
band,  waveform  distortion  therefore  implies  complex  wavefront 
amplitude  distortion,  and  vice  versa.  In  the  following  text,  wave¬ 
form  distortion  and  wavefront  amplitude  distortion  are  used  indis- 
tinguishably. 

C,  Correction  upon  Scattering  Samples.  To  develop  an  un¬ 
derstanding  of  deaberration,  our  studies  have  concentrated  on 
one-way  propagation  measurements  and  corrections  of  point 
source  fields.  The  two-dimensional  (2D)  measurement  system 
and  procedure  were  described  elsewhere  [11,16].  Briefly,  a  breast 
specimen  was  placed  between  hemispherical  source  and  receiving 
array.  The  2D  array  consisted  of  a  92-mm  ID  linear  array  trans¬ 
lated  46  mm  perpendicular  to  its  axis  to  form  a  synthetic  2D 
array  92  X  46  mm.  Element  pitch  in  the  receiving  transducer  was 


Vol.  8,  322-335  (1997)  325 


(a)  -6dB  (water) 


(c)  -20  dB  (brsOOS) 


0.72  mm  and  a  reflecting  mask  reduced  the  receiving  elevation 
to  1.44  mm.  Frequency  was  3.7  MHz.  Bandwidth  was  2  MHz. 
Waveforms  were  measured  at  each  element,  from  which  2D 
wavefronts  were  reconstructed  as  functions  of  time.  3D  image 
data  were  calculated  by  ( 1 )  Fourier  decomposition  of  the  tempo¬ 
ral  waveforms  at  each  (jc ,  y)  position  in  the  aperture,  (2)  calcula¬ 
tion  of  a  complex  CW  2D  image  at  each  Fourier  frequency  in 
the  focal  plane  by  using  the  angular  spectrum  technique  [18], 
and  (3)  summation  of  the  2D  images  to  form  the  3D  transient 
image.  The  final  2D  image  used  in  this  report  was  obtained  by 
detecting  the  peak  pressure  value  at  each  {xf,  yf)  position  in 
the  image  plane  within  the  transient  period.  Contour  maps  of  2D 
point  source  images  at  different  thresholds  were  used  to  quantita¬ 
tively  evaluate  image  quality.  The  implementations  of  different 
correction  methods  and  correction  results  were  reported  in  [11]. 
The  total  measurements  made  were  16  with  breast  sample  thick¬ 
ness  varied  from  1  to  4  cm.  The  following  is  a  typical  example 
of  correction  upon  a  scattering  sample. 

Figure  7(a)  and  7(b)  shows  —  6-dB  contour  plots  of  images 
obtained  from  the  water  path  and  the  3.5-cm  breast  tissue  path 
(brs005).  Contour  spacing  is  1  dB.  Abscissa  and  ordinate  are 
azimuth  and  elevation  in  millimeters  in  the  image  plane.  Figure 


(b)  -6  dB  (brsOOS) 


Figure  7.  2D  contour  maps  of  point  source  images, 
(a)  Water  data  showing  diffraction  pattern  of  system. 
Outer  contour  is  -6  dB  level,  (b)  Measured  through 
the  3.5-cm  tissue  (brsOOS),  showing  asymmetric  scat¬ 
tering  effect.  Outer  contour  is  -6  dB  level,  (c)  Image 
brs005  at  -20  dB  contour  level.  The  outer  contour 
shows  symmetrical  scattering  pattern.  No  distinct 
structure  or  multipath  is  observed.  Note  the  change  of 
scale  from  (a)  and  (b). 


7(a)  shows  the  system  diffraction  pattern  at  -6  dB.  The  azi¬ 
muthal  beamwidth  is  about  1.5  mm  (at  a  distance  of  180  mm), 
meaning  that  the  point  (or  lateral)  resolution  is  1.5/180  or  8.3 
mrad.  Because  the  size  of  the  receiving  array  in  elevation  is  half 
the  size  in  the  array  direction,  the  width  of  the  image  in  elevation 
is  twice  as  large  as  in  azimuth.  When  tissue  is  present  (b),  the 
image  lobe  is  broadened.  In  (b)  the  beamwidth  has  grown  from 
1.5  to  3  mm  in  azimuth  and  3.1  to  9,3  mm  in  elevation,  and  the 
point  resolution  has  worsened  approximately  2: 1  in  azimuth  and 
3:1  in  elevation.  A  symmetrical  scattering  pattern  appears  when 
the  threshold  is  reduced.  Figure  7(c)  shows  a  — 20-dB  contour 
plot  of  the  sample  image  brs005.  Contour  spacing  is  2  dB.  Energy 
is  spread  out  over  a  large  area.  No  distinct  structure  or  multipath 
lobe  can  be  identified  from  the  image. 

Figure  8(a-g)  shows  the  correction  result  of  the  sample  im¬ 
age  brs005  by  using  different  methods.  Part  (h)  is  the  -30-dB 
level  contour  of  a  water  path  image.  The  contour  spacings  in 
Figure  8  are  the  outer  contour  levels  divided  by  10.  Figure  8(a- 
c)  shows  TDC  at  -6-,  -20-,  and  -30-dB  levels  [see  Fig.  7(b) 
and  7(c)  for  comparison].  Image  quality  improved  significantly 
and  the  system  diffraction  pattern  is  restored  up  to  “20-dB  level. 
The  symmetrical  scattering  pattern  outlined  by  outer  contours  at 
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(b) 


TDC  at  -  20  dB 


(d) 


TDC  at  -30  dB 


PC  at  -30  dB 


(f) 


PC  and  square  rooter 
at  -30  dB 


(h) 


Figure  8.  Distortion-corrected  contour  maps  of  the  sample  image  brs005  [see  Fig.  7(b)  and  4(c)  for  comparison],  (a)  TDC  at  -6  dB.  (b) 
TDC  at  -20  dB.  (c)  TDC  at  -30  dB.  (d)  PC  at  -30  dB.  (e)  BPT  at  -30  dB.  (f )  PC  and  square  rooter  at  -30  dB.  (g)  Wideband  DSA  at  -30 
dB.  (h)  Water  at  -30  dB.  Note  the  scale  change  of  (b-h)  from  (a). 


the  —  30-dB  level  is  caused  by  residual  phase  errors  after  phase 
deaberration  and  incoherent  amplitude  distortion.  Figure  8(d)  is 
the  correction  result  of  PC  at  the  -30-dB  level.  Because  of  the 
complete  removal  of  phase  errors  except  for  a  27r  jump,  the 


scattering  pattern  is  further  improved  as  compared  to  Figure  8(c). 
Figure  8(e)  is  the  result  of  BPT,  and  the  scattering  pattern  is 
much  better  than  that  of  PC  at  the  aperture.  The  optimum  back- 
propagation  distance  of  this  tissue  sample  is  30  mm.  It  is  interest- 
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Figure  9.  The  field  is  the  original  field  plus  refracted  fields. 


ing  to  note  that  the  beamwidth  is  narrower  than  that  obtained 
through  other  methods.  The  reason  is  that  the  actual  aperture  is 
moved  30  mm  closer  to  the  source  and  the  beamwidth  is  (/^-optimal 
backpropagation  distance)  \/L  instead  of  R(\/L),  where  R  is  the 
distance  between  source  and  the  receiving  aperture,  and  L  is  the 
aperture  size.  Figure  8(f)  is  the  correction  result  of  PC  and  square 
rooter  on  wavefront  amplitude  at  -30  dB.  The  pattern  is  very  close 
to  the  system  diffraction  pattern  [see  Fig.  8(h)  for  comparison]. 
Figure  8(g)  is  the  result  of  wideband  DSA  at  the  -30-dB  level. 
The  correction  result  is  worse  than  those  of  other  methods. 

Correction  statistics  of  a  simple  metric,  the  lowest  level  at 
which  main-lobe  diffraction  shape  is  maintained  when  different 
methods  are  used,  are  given  in  [1 1] .  The  average  levels  of  using 
wideband  DSA,  PC,  and  PC  and  rooter  are  14.8  (standard  devia¬ 
tion  2.34),  19.5  (2.61),  and  25.4  dB  (2.55  dB),  respectively. 
These  results  closely  agree  with  the  analysis  in  Section  II  of  a 
6-dB  decrease  and  increase  in  image  contrast  from  squaring  the 
amplitude  (DSA)  to  amplitude  intact  (PC)  and  amplitude  intact 
to  square  root  of  the  amplitude  (PC  and  rooter).  The  average 
levels  of  using  TDC  and  BPT  are  19.3  (2.43)  and  24.0  dB  ( 1.99) . 

ill.  SCATTERING  AND  REFRACTION  MODEL 
A.  Strong  Scattering  and  Refraction  Model.  In  a  refractive 
medium,  coherent  rays  can  be  bent  and  split  after  they  pass 
through  tissue  beds  with  different  sound  speeds  [10].  The  inter¬ 
ference  of  refracted  coherent  ray  bundles  with  incident  rays  pro¬ 
duces  image  artifacts  which  appear  as  false  targets.  A  simple 
representation  of  such  a  field  is 

E{x) 

=  EaEiE2-  +  £,2(x)^f'AW^  (8) 

where /i(x)  and/2(x)  define  directions  of  secondary  fields  caused 
by  refraction  (Fig.  9). 

The  angular  source  intensity  distribution  is  the  superposition 
of  the  intensity  image  lobe,  coherent  interference  lobes,  and  inco¬ 


herent  spectrum  [Fig.  10(a)].  The  strengths  and  the  numbers  of 
interference  lobes  depend  on  orientations,  curvatures,  sizes,  and 
numbers  of  refractive  bodies  in  the  insonified  medium,  and  may 
not  be  directly  related  to  propagation  depth.  Although  the  chance 
of  incident  rays  passing  through  more  refractive  bodies  increases 
with  the  propagation  depth,  the  strengths  of  interference  lobes 
resulting  from  multiple  refraction  may  not  be  significant  com¬ 
pared  with  those  resulting  from  single  refraction.  One  quantitative 
measure  of  coherent  interference  phenomenon  is  the  ratio  of  a 
large  interference  lobe  or  artifact  to  image  lobe.  For  nine  samples 
that  were  determined  to  have  significant  levels  of  coherent  inter¬ 
ference,  the  ratio  varied  randomly  from  sample  to  sample  and 
showed  no  increase  with  thickness.  The  average  strength  (peak 
value)  is  10  dB  (standard  deviation  4  dB )  [1 1] .  Another  quantita¬ 
tive  measure  is  the  total  number  of  large  refractive  lobes  in  each 
image.  The  average  is  2  and  standard  deviation  is  0.53.  The  lobe 
number  is  independent  of  propagation  depth.  The  refractive  lobes 
are  tightly  clustered  around  the  main  lobe  with  an  average  radius 
of  5.6  mm  (standard  deviation  1.8),  which  is  about  2°.  This 
is  because  the  speed  variations  across  glandular  tissue  and  fat 
boundaries  are  small  (5-10%). 

The  correction  mechanism  in  a  strong  scattering  and  refraction 
medium  is  interesting.  If  the  primary  field  is  stronger  than  second¬ 
ary  refraction  fields,  the  scattered  energy  removed  by  the  deaber¬ 
ration  process  strengthens  the  primary  field  by  folding  the  re¬ 
moved  scattered  energy  into  the  primary  field.  As  a  result  the 
ratio  of  image  lobe  to  interference  lobe  is  improved  [Fig.  10(b)] . 
If,  however,  the  primary  field  is  weaker  than  the  refracted  field, 
the  scattered  energy  removed  by  the  deaberration  process  may 
strengthen  the  refracted  field.  The  corrected  result  could  then  be 
worse  than  without  correction.  However,  we  have  seen  no  evi¬ 
dence  of  this  so  far  in  our  experiments  reported  in  [11]. 

Time-reversal  mirror  or  wideband  DSA  on  reception  is  a 
matched  filter  process  in  both  space  and  time,  resulting  in  maxim¬ 
ization  of  the  signal  relative  to  the  background  white  noise.  Corre¬ 
lation-based  time-delay  compensation  algorithms,  on  the  other 
hand,  perform  a  different  function;  they  adjust  time  delay  to 
maximize  image  sharpness  or  image  contrast.  In  a  refractive  me¬ 
dium,  the  problem  is  more  complicated.  Neither  TRM  or  wide¬ 
band  DSA  nor  time-delay  algorithms  have  a  way  of  knowing 
which  signal  is  real  and  which  is  an  artifact  caused  by  coherent 
multipath.  Because  TRM  or  wideband  DSA  is  a  spatial  matched 
filter,  its  weight  vector  forms  a  lobular  radiation  pattern  that 
matches  the  sum  of  both  the  real  and  false  signals.  Hence  all 
targets,  whether  real  or  false,  are  enhanced.  In  our  experiments 
reported  in  [11],  with  nine  samples  that  have  well-defined 
multipath  artifacts,  we  found  an  average  5  dB  ( standard  deviation 
2.5 )  increase  in  peak  artifact  to  image  lobe  ratio  by  using  wide¬ 
band  DSA  rather  than  by  using  time-delay  compensation. 
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Figure  10.  Before  (a)  and  after  (b)  phase  compensation. 
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BPT  backpropagates  waveforms  to  the  optimum  distance  be¬ 
fore  applying  time-delay  compensation.  It  removes  wavefront 
amplitude  distortion  caused  by  propagation  and  reduces  residual 
phase  distortion  by  improving  the  performance  of  time-delay  cor¬ 
rection,  but  will  not  remove  large  lobes  caused  by  refraction.  A 
simple  thought  experiment  is  convincing:  Backpropagating  two 
plane  waves  produces  two  plane  waves  at  whatever  distance  the 
waves  are  backpropagated. 

Thus,  it  is  important  to  understand  the  deaberration  operations 
and  to  design  appropriate  algorithms  that  offer  optimal  perfor¬ 
mance  in  the  corresponding  distorting  medium.  Aperture  mea¬ 
surement  alone  cannot  distinguish  between  wavefront  distortions 
caused  by  refraction  or  by  strong  scattering.  Distinction  can  be 
made  in  the  image  domain  where  refraction  appears  as  large, 
coherent  interference  lobes  while  strong  scattering  spreads  into 
an  incoherent  background.  While  the  refracted  energy  may  or 
may  not  be  significantly  larger  than  the  energy  in  the  incident 
wave,  it  nevertheless  produces  image  artifacts  which  will  destroy 
the  quality  of  a  high  resolution  image. 

B.  Correction  upon  Samples  with  Well-Defined  Multipath. 

In  most  breast  samples  (nine  of  16),  scattering  and  coherent 
interference  are  both  present.  Figure  11(a)  and  11(b)  is  -10- 
dB  contour  plots  of  images  obtained  from  the  water  path  and  the 
4-cm  breast  tissue  path  (brs006).  Contour  spacing  is  2  dB.  In 
tissue  image  (b),  three  lobes  instead  of  one  central  image  lobe 
appear.  Inner  contours  of  the  image  lobe  are  close  to  the  system 
diffraction  pattern  [Fig.  11(a)],  while  the  outer  contours  are 
highly  irregular  owing  to  scattering.  The  point  resolution  is 
worsened  by  approximately  2:1  in  azimuth  and  3:1  in  elevation. 
Two  refractive  lobes  appear  at  (—4,  —2.5)  and  (3.5,  1.5),  and 
their  strengths  are  —6  and  -8  dB.  A  closer  examination  is  ob¬ 
tained  from  45°  cuts  (Fig.  12)  of  the  Figure  11(a)  and  11(b) 
images.  In  Figure  12(a),  the  peak  sidelobe  is  -36  dB,  which  is 
the  best  we  can  obtain  from  the  measurement  system.  In  the 
tissue  image,  two  large  side  lobes  (—6  and  -8  dB)  appear  on 
either  side  of  the  main  lobe.  The  two  nulls  between  the  lobes  are 
deep  (left:  —13  dB  below  peak,  right:  -15  dB),  which  is  a 
typical  coherent  interference  phenomenon  [10]. 

A  symmetrical  scattering  pattern  generally  appears  when  the 
threshold  is  reduced,  as  shown  in  Figure  11(c)  with  a  —  16-dB 
threshold  contour  and  contour  spacings  of  approximately  1.5  dB. 
The  outer  contours  show  a  roughly  symmetrical  pattern  typical  of 
scattering,  while  the  inner  contours  show  an  asymmetric  lobular 
interference  pattern. 

Figure  13  shows  correction  results  of  the  sample  image  brs006 
after  applying  different  methods.  The  contour  spacings  are  the 
outer  contour  levels  divided  by  10.  Figure  13(a)  is  the  -16-dB 
contour  map  of  applying  TDC  [see  Fig.  11(c)  for  comparison]. 
Image  quality  improved  significantly  because  of  the  minimization 
of  the  phase-front  distortion  caused  by  scattering.  The  area  within 
the  -16-dB  inner  contour  is  reduced  by  a  factor  of  11.  The 
mainlobe  pattern  is  close  to  the  system  diffraction  pattern,  al¬ 
though  the  area  is  still  1 .6  times  larger  than  the  area  in  water  at 
the  same  level.  Two  large  lobes  remain  as  expected  and  appear 
as  two  additional  sources.  Figure  13(h)  is  the  correction  result 
at  a  -16-dB  level  of  applying  wideband  DSA.  It  is  evident  that 
the  procedure  enhances  interference  in  this  example;  the  artifacts 
are  5  dB  higher  than  in  Figure  13(a).  Figure  13(b)  and  13(c) 
shows  -25-  and  -30-dB  contour  plots  of  TDC.  In  Figure  13(b) 
the  outer  contours  show  an  asymmetric  interference  pattern  with 


scattered  energy  distributed  around  the  two  interference  lobes. 
In  Figure  13(c)  the  outer  contours  show  the  more  symmetrical 
scattering  pattern  caused  by  incoherent  wavefront  amplitude  dis¬ 
tortion  and  residual  phase  distortion.  Figure  13(d)  is  the  result  of 
PC  at  -30  dB.  The  approximately  symmetrical  scattering  pattern 
outlined  by  outer  contours  in  Figure  13(c)  is  improved  because 
of  better  phase  deaberration.  The  energy  is  more  concentrated 
around  the  image  lobe  and  interference  lobes.  Part  (e)  is  the 
result  of  BPT,  and  the  scattering  pattern  is  further  improved  as 
compared  with  the  result  of  PC  at  the  aperture.  The  backpropaga- 
tion  distance  of  this  sample  is  50  mm.  Part  (f )  is  the  result  of 
PC  and  square  rooter  at  -30  dB.  The  pattern  is  very  close  to  the 
system  diffraction  pattern.  PC  and  fourth  rooter  shown  in  Figure 
13(g)  restores  the  diffraction-limited  image  to  a  -30-dB  level. 

It  is  interesting  to  note  that  the  strengths  of  the  two  interference 
lobes  are  also  reduced  by  8  dB  after  TDC  and  PC.  They  are 
further  reduced  by  5,  7,  and  15  dB  after  BPT,  PC  and  square 
rooter,  and  PC  and  fourth  rooter,  respectively.  The  improvement 
comes  from  the  scattered  energy  which  is  originally  distributed 
around  the  main  lobe  and  each  of  the  refracted  multipath  lobes. 
This  energy  is  coherently  added  back,  by  the  deaberration  pro¬ 
cess,  to  the  dominant  lobe,  which  in  this  case  is  the  main  lobe. 
The  overall  result  is  an  improvement  of  the  image  lobe  to  the 
interference  lobe  ratio. 

Determination  of  whether  interference  was  coherent  or  scat¬ 
tered  was  based  on  whether  large  lobes  could  be  removed  by 
applying  phase-deaberration  algorithms.  Another  indicator  is 
whether  a  lobe  location  moves  (scattering)  or  remains  relatively 
the  same  (refraction)  after  the  phase  deaberration  process.  Both 
means  were  used  to  identify  whether  these  lobes  were  produced 
by  refraction  or  scattering.  Large  lobes  in  nine  samples  are  be¬ 
lieved  to  be  caused  by  the  interference  process.  In  these  samples, 
we  have  found  an  average  of  a  10-dB  (standard  deviation  4) 
improvement  in  the  range  of  5-16  dB  after  TDC,  owing  to  the 
coherent  strengthening  of  the  mainlobe. 

IV.  EXPERIMENTS 

The  significance  of  appropriate  amplitude  weighting  discussed  in 
Section  II  is  experimentally  evaluated  using  ER  introduced  in 
Equation  (5).  When  phase  error  is  removed,  ER  is  proportional 
to  a  A,.  The  detailed  measurement  procedures  of  ER  and  are 
reported  in  [11].  Figure  14  shows  regression  plots  of  ER,  mea¬ 
sured  after  reducing  the  phase  error  using  PC,  wideband  DSA, 
and  PC  and  square  rooter,  versus  a\.  The  three  methods  compen¬ 
sate  phase  errors.  ER  and  a\  measured  in  water  using  the  three 
correction  procedures  are  also  incorporated  into  the  three  data 
sets  to  represent  the  system  performance.  The  linear  relationships 
are  statistically  significant.  The  average  ERs  are  103%  (standard 
deviation  33.9%),  25%  (6.9%)  and  6.8%  (1.6%)  when  wideband 
DSA,  PC,  and  PC  and  rooter  are  used,  respectively.  An  approxi¬ 
mate  factor  of  4  reduction  in  ER  from  amplitude  squaring  (wide¬ 
band  DSA)  to  amplitude  intact  (PC)  and  amplitude  intact  to 
square  root  of  amplitude  (PC  and  rooter)  is  obtained.  The  results 
agree  with  the  analysis  in  Section  II  and  indicate  that  intelligent 
use  of  amplitude  weighting  is  crucial  for  improving  image  con¬ 
trast. 

Scattered  energy  increases  with  the  propagation  depth  as 
predicted  by  Equation  (5)  and  performances  of  descattering 
algorithms  deteriorate  with  propagation  depth.  In  [11],  regres¬ 
sion  plots  of  ER  versus  approximate  thickness  of  16  breast 
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(c)  -16  dB  (brs006) 

Figure  11.  2D  contour  maps  of  point  source  images,  (a)  water  data  showing  diffraction  pattern  of  system.  Outer  contour  is  “10  dB  level, 
(b)  measured  through  the  4-cm  tissue  (brs006),  showing  highly  asymmetric  interference  pattern.  Outer  contour  is  -10  dB  level,  (c)  Image 
brs006  at  -16  dB  contour  level.  The  outer  contour  shows  more  symmetrical  scattering  pattern.  Note  the  change  of  scale  from  (a)  and  (b). 


samples  were  given  by  using  seven  descattering  procedures. 
The  regression  plot  without  correction  was  also  given.  The 
linear  relationships  were  statistically  significant.  Figure  15 
shows  partial  results  taken  from  [11].  The  five  curves  in  the 
figure  are  obtained  by  using  wideband  DSA  (circle),  TDC 


(solid  square),  PC  (triangle),  PC  and  rooter  (solid  triangle), 
and  no  correction  ( solid  circle ) .  ER  reduction  from  amplitude 
weighting  (wideband  DSA)  to  square  root  of  amplitude 
weighting  (PC  and  rooter)  is  clearly  shown  by  the  regression 
plots.  The  regression  lines  are  extrapolated  to  100-mm  propa- 


(a)  (b) 

Figure  12,  (a)  45^"  cuts  of  2D  Images:  (a)  water,  (b)  brs006. 
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Figure  13.  Distortion  corrected  contour  maps  of  the  sample  image  brs006  [see  Fig.  11(b)  and  11(c)  for  comparison],  (a)  TDC  at  -16  dB  level, 
(b)  TDC  at  -25  dB  level  and  (c)  TDC  at  -30  dB  level,  (d)  PC  at  -30  dB  level,  (e)  BPT  at  -30  dB.  (f)  PC  and  square  rooter  at  -30  dB.  (g)  PC 
and  fourth  rooter  at  -30  dB.  (h)  Wideband  DSA  at  -16  dB.  Artifacts  are  5  dB  higher  than  time-delayed  correction  shown  in  Figure  13(a). 


gation  depth,  which  is  considered  the  average  maximum  pene-  between  ER  and  image  contrast,  which  is  defined  as  the  lowest 

tration  depth  required  for  breast.  level  of  the  main-lobe  diffraction  shape.  Image  contrast  increases 

Before  we  discuss  the  predicted  performance  of  the  descattering  monotonically  with  ER.  Table  I  lists  the  measured  average  ER 

algorithms  at  100-mm  depth,  we  need  to  find  the  relationship  (first  row)  and  the  measured  average  image  contrast  (second  row) 


Vol.  8,  322-335  (1997)  331 


i 


Figure  14.  Linear  regression  plots  of  ER,  after  wideband  DSA  (upper),  phase  conjugation  (middle)  and  phase  conjugation  and  square  rooter 
(bottom),  versus  wavefront  amplitude  variance. 


when  different  correction  methods  as  well  as  no  correction  are 
used  [11].  The  relationship  between  ER  and  image  contrast  can 
be  predicted  by  using  ER  and  image  contrast  obtained  from  PC 
(or  any  correction  method)  as  baseline  (bold  numbers  in  Table 
I).  For  example,  the  measured  image  contrast  with  PC  is  19.5  dB. 
The  ER  ratio  of  PC  alone  to  PC  and  rooter  is  25.1  1/6.8  =  3.7, 
corresponding  to  an  increase  of  5.7  dB.  Thus,  the  predicted  image 


contrast  is  19.5  +  5.7  =  25.2  dB,  which  closely  agrees  with  the 
measured  value  25.4  dB.  The  third  row  in  Table  I  lists  the  predicted 
image  contrast  based  on  the  measured  ER.  The  predicted  values 
are  in  reasonably  good  agreement  with  the  measured  values  except 
for  the  pair  without  correction.  The  reason  is  due  to  ER  measure¬ 
ment  error  in  images  without  correction.  For  some  severely  dis¬ 
torted  images,  it  is  often  difficult  to  identify  the  mainlobe  region. 


propagation  depth  (mm) 

Figure  15.  Linear  regression  plots  of  ER  versus  propagation  depth  obtained  from  without  correction  (solid  circle),  wideband  DSA  (circle), 
TDC  (solid  square),  PC  (triangle),  and  PC  and  rooter  (solid  triangle).  The  correlation  coefficients  are  given  by  R.  One  outlier  was  removed 
from  (a). 
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Table  I.  Relation  of  ER  and  image  contrast  (predicted  and  measured). 


Without 

Corr. 

TDC 

PC 

Wideband 

DSA 

BPT 

PC  & 
Square 
Rooter 

BPT& 

Square 

Rooter 

PC  & 
Fourth 

Rooter 

ER 

214.7% 

44.2% 

25.1% 

103.5% 

12.8% 

6.8% 

3.7% 

2.1% 

Image  contrast 
Measured 

4.2  (dB) 

19 

19.5 

14.8 

24.0 

25.4 

29.4 

31.1 

Predicted 

10.2  (dB) 

17 

19.5 

13.3 

22.4 

25.2 

27.8 

30.3 

At  100-mm  depth,  ER  after  PC  reaches  85%  which  corre¬ 
sponds  to  14.0  dB  image  contrast.  Therefore,  PC,  the  best  linear 
phase  correction  procedure  at  the  aperture,  is  likely  to  result  in 
only  14  dB  image  contrast  at  this  depth.  ER  after  wideband 
DSA  reaches  324.3%,  which  corresponds  to  only  ~  8.4  dB  image 
contrast.  ER  after  PC  and  rooter  reaches  22%,  which  provides  ~ 
20.1  dB  image  contrast.  ER  after  TDC  reaches  148%,  which 
offers  only  ~  11.8  dB  image  contrast. 

The  data  used  to  test  the  analysis  developed  in  Section  II  are 
sparse;  nevertheless,  they  provide  information  for  designers  of 
next-generation  high-quality  breast  scanners. 

V.  DISCUSSION 

A.  Application  to  Diffuse  Scattering  Medium.  The  com¬ 
pression  operation  upon  wavefront  amplitude,  in  addition  to  phase 
deaberration  or  BPT,  introduced  here  is  applicable  to  a  diffuse 
scattering  medium.  Experiments  with  TDC  and  rooter  and  TDC 
alone  show  that  ER  can  be  reduced  from  44.2%  to  28.6%,  the 
ratio  of  which  indicates  1.3  dB  improvement  in  image  contrast. 
However,  experiments  with  PC  and  rooter  and  PC  alone  show 
that  ER  can  be  reduced  from  25.1%  to  6.8%,  which  is  about  6 
dB  improvement  in  image  contrast.  The  large  difference  is  due 
to  the  fact  that  PC  corrects  phase  in  every  signal  spectrum,  while 
TDC  makes  a  single  correction  which  applies  to  all  frequencies. 
The  consequence  is  that  TDC  leaves  significant  phase  residual. 
Therefore,  to  successfully  improve  image  contrast,  good  phase 
deaberration  is  the  first  crucial  step  to  reduce  ER,  followed  by 
amplitude  compression  such  as  rooter. 

B.  Application  to  Coherent  Interference  upon  Pulsed- 
Echo  Imaging.  High-resolution,  high-dynamic  range  breast 
imaging  implies  seeing  a  black  cyst  60  dB  weaker  than  a  bright 
speckled  target.  In  this  article,  we  show  that  30  dB  image  contrast 
on  reception  is  obtainable  by  taking  low-order  rooters  provided 
that  phase  error  is  completed  compensated.  In  pulsed  echo,  re¬ 
fracted  energy  from  the  targets  insonified  by  the  subbeam  will 
arrive  at  angles  of  1-3°  from  the  target  direction  [11]  and  there¬ 
fore  in  the  side-lobe  region  of  array.  The  side-lobe  level  of  the 

receiving  beam  will  be - 30  dB.  The  average  level  of  refracted 

multipath  is  found  to  be  -10  dB  relative  to  the  target  image 
[11],  and  therefore  the  average  multipath  signal  level  entering 
the  system  is  ~  ( 10  +  30)  dB  =  —40  dB  below  the  echo  strength 
of  the  illuminated  target.  Cancelation  algorithms  [12]  must  sup¬ 
press  these  multipath  signals  to  achieve  60  dB  contrast  resolution. 
Hence,  about  20  dB  further  cancelation  is  required. 

VI.  CONCLUSIONS 

We  conclude,  through  analysis  and  experiment,  that  ( 1 )  matched 
filtering  compensation  algorithms  while  optimal  for  detection  are 


not  suitable  for  imaging  in  distributed  scattering,  because  the 
algorithms  increase  wavefront  amplitude  variance  and  therefore 
decrease  image  contrast;  and  (2)  low-order  rooters  upon  wave- 
front  amplitude,  akin  to  inverse  filtering,  in  addition  to  phase 
deaberration  can  significantly  improve  image  contrast  without 
encounter  stability  problem.  Further  study  is  needed  to  optimize 
the  required  degree  of  compression  as  a  function  of  SNR  and 
image  contrast. 

Incoherent  scattering  and  coherent  multipath  affect  an  image 
differently.  Scattering  reduces  target  strength,  broadens  the  image 
lobe,  and  increases  background  level,  while  refraction  creates 
coherent  multipath  interference  which  produces  false  targets  or 
image  artifacts.  Adaptive  compensation  algorithms  developed  so 
far  are  useful  to  partially  remove  the  scattered  energy  and  coher¬ 
ently  add  it  back  to  the  coherent  image.  The  result  is  the  improve¬ 
ment  of  both  lateral  and  contrast  resolution  to  a  large  extent.  The 
ratio  of  image  lobe  to  artifacts  is  also  improved  because  of  the 
strengthening  of  the  coherent  filed  by  descattering  processes. 
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APPENDIX 

Let  5(w)  represent  the  complex  angular  source  distribution,  where 
u  =  sin(^)  and  0  is  the  angle  from  the  array  normal.  Its  complex 
radiation  field  in  the  axis  of  the  receiving  array  in  a  homogeneous 
medium  is  the  inverse  Fresnel  transform  of  s(w)  or  inverse  Fou¬ 
rier  transform  of  s{u)  after  correcting  for  near-field  curvature 
and  is  denoted  as  Eo{x).  In  an  inhomogeneous  medium,  the  radia¬ 
tion  field  is  i{x)  =  aW {x)Eo{x)M{x),  where  w(x)  is  the  aper¬ 
ture  weighting  function,  M(x)  =  is  the  medium- 

induced  distortion,  and  a  is  a  constant  to  reflect  any  signal  loss. 
The  complex  image  is 

Si{u)  =  aF[W  {x)Eo{x)M{x)] 

=  a/(M)*5(w)*F(M(x))  (Al) 

where  F(M(x))  =  //(«)  is  the  angular  impulse  response  of  the 
medium.  F  denotes  Fourier  transform.  The  subscript  i  denotes 
inhomogeneous  medium.  For  a  =  \  and  M{x)  =  1,  the  complex 
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image  corresponds  to  the  image  in  homogeneous  medium  and  is 
denoted  as  Sh{u). 

The  derivation  is  made  with  discrete  rather  than  continuous 
variables.  Thus,  jc  is  replaced  by  n  =  1,  2,  •  •  •  where  N  = 
number  of  elements  in  the  array.  The  following  assumptions  are 
made: 

1.  Medium-induced  amplitude  and  phase  errors  are  uncorre- 
lated. 

2.  The  phase  errors  of  (j)t(n) ,  n  =  1,2,  •  •  •  A^,  are  statistically 
uncorrelated  random  variables  with  a  common  probability 
density  function  (pdf). 

3.  The  amplitude  errors  of  A[(n),n  =  1,2,  *  •  *  A^,  are  statisti¬ 
cally  uncorrelated  random  variables  with  a  common  pdf. 

The  complex  image  Si(u)  can  be  written  as 

N 

Si  =  aF[Eo(n)At(n)]  =  a  X  Eo(n)At(n)oxp(jkndu)  (A2) 

n= ! 

where  the  aperture  taper  W {n)  has  been  chosen  as  unity. 

The  image  intensity  is  the  product  of  Equation  (A2)  and  its 
complex  conjugate: 

N  N 

Si(u)st(u)  =  X  Z  E„{n)Ef(m)A,(n)A,(m)exp(j(<t>,(n) 

n=l  m=l 

—  -  m)du) 


exp(y^,(?z))exp(“y<;)/(m))  =  ox^{j(i),{n))tx^{~  j4>Am)) 

because  of  assumptions  (2)  and  (3),  respectively. 

{N  N 

X  Ea{n)Et{n){crl)  +  S  Eo{n)E*{n) 

n=m  n—m 


X  (1  -  exp(7(0,)^)(A)^  +  A^exp(7(</>,)^  Z  Z  E^(n)E* 

n=l  m  =  l 

X  (m)oxp(jk{n  -  (A5) 


The  term  A,^exp(y(0r)^  Ea(n)Ef(n)  is  added  to  the  sec¬ 
ond  term  of  (A4),  resulting  in  the  third  term  of  (A5).  The  same 
term  is  subtracted  from  (A4),  resulting  in  the  second  term  of 
(A5). 


)(u)sf^(u)  =  aH  X  Eo(n)Ef(n)(aX)  +  X  E^(n)Ef{n) 


=  Z  Ea{n)E^(n)Aj(n) 

N 

+  S  Ea(n)E*{m)At(n)Af(m)txp(j((l)An) 


X  (1  -  exp(y(<^,)^)(A,)^  -h  iA,y&xp(j'((l>i)%Xu)stiu)\ 


(A6) 

Integration  with  respect  to  u  on  both  side  of  Equation  (A6)  yields 
I  Si(u)sf'{u)  du  +  I  Si{u)sf{u)  du 

V  main  lobe  v  side  lobe 


1 

=  a4  f 

1  (A3) 

v  main  lobe 

The  average  intensity  is 


{N  _  N  r  r  r  ^ 

X  Eo(n)Ef{n)A^(n)  +  X  E„(n)E^(m)  ^  (u)st(u)  du  +  a^i  I  X  E„{n)Et{n 

n=m  n*m  side  lohe  L 


X  (1  -  &Xp{j{<f>y){A, 

X  ( 


X  £,(«)£*(n)(al,)  +  S  EAn)E*{n) 

n=m 

y]  du+  ( 

J  V  me 


(A,)^exp(;(0,)  % 


0«) 


X  A,(n)A,(m)e\p(j(cj>,(n))e.xp(-j4>,(m)) 
X  cxp(jk(n  —  m)du) 


■} 


Z  EAn)Enn){y, 


+  X  Eo{n)Etin)i\  -  exp(7(</>,)^)(A,)2 

n=m 

I 


du 


+  I  {AAQxp{j{(l}t)%(u)sf(u)du 

side  lobe 


(A7) 


+  (Ay)  +  X  EMEfim) 


(At)  exp(y(</>r)  exp(y/:(/2  -  m)du)  >  (A4) 


A  Energy  ~  outside  main  lobe 
_ - 

Energy  —  inside  main  lobe 
Lideione  EAn)E*(n)  du} 


\aX+  -  exp(;(<^>,)  )A/ 


where 


and 


Lainlobe  (  A, )^eXp(y  (</>,))  (  M)#  (tt)  dU 


(A8) 


At(n)At(m)  =  At(n)At{m)  =  Ai^ 


where  (A,)^exp(;((^,) \{u)s^(u)  du  is  small  compared 

with  the  rest  of  the  side-lobe  terms  and  is  dropped  from  the 
numerator  and  Eo{n)E* {n){crX)  + 
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Eo{h)E^{n){  \  —  exp(7(</>/)  ^)(  du  is  small  compared  with 
anothe£jmain-lobe  term  and  is  dropped  from  the  denominator. 
Since  is  close  to  unity  and  exp(y‘(<^f)  is  also  close  to  unity 
when  the  phase  error  is  moderate,  the  right  side  of  Equation  ( A8) 
is  C  {  +  ( 1  -  exp(y(<^,)  } ,  where 

^  Ssi,e,o^^''^-^EAn)Enn)du 

which  is  related  to  the  source  properties  only.  When  the  phase  error 
is  moderate,  1  ~  exp(J((f>,)  ^  ^  al^  and  ER  «  C  { } . 
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Abstract — Techniques  for  the  correction  of  ultrasonic 
wavefront  distortion  are  compared  using  measured  pulse 
transmission  through  human  breast  specimens.  The  mea¬ 
sured  data  were  obtained  by  recording  the  pulse  signals  on 
each  element  of  a  linear  array  that  was  moved  in  elevation  to 
synthesize  a  two-dimensional  aperture.  The  one-way  point 
spread  functions  were  reconstructed  using  matched  filter¬ 
ing  technique,  phase  conjugation,  time-delay  compensation, 
and  backpropagation  followed  by  phase  conjugation,  each 
with  and  without  amplitude  compression.  Two  measures  of 
performance  were  compared,  the  mainlobe  diffraction  shape 
and  the  ratio  of  the  energy  outside  the  mainlobe  to  the  en¬ 
ergy  inside  the  mainlobe.  Matched  filtering,  which  compen¬ 
sates  the  phase  distortion  bat  also  increases  the  variation 
in  the  modulus,  performs  more  poorly  than  other  tech¬ 
niques  with  respect  to  both  measures.  Phase  conjugation 
and  time-delay  compensation,  which  leave  the  magnitude 
of  the  wavefront  unchanged,  have  similar  beamwidths  but 
phaise  conjugation  is  consistently  superior  with  respect  to 
energy  ratio.  The  backpropagation  method,  which  models 
wavefront  distortion  using  a  phase  screen  at  a  computed 
position  between  the  source  and  aperture,  is  shown  to  per¬ 
form  better.  The  use  of  a  novel  amplitude  compression  that 
approaches  inverse  filtering  improves  the  performance  of 
the  compensation  techniques  significantly.  This  is  because 
inverse  filtering  optimizes  image  fidelity,  in  contrast,  for  ex¬ 
ample,  to  matched  filtering,  which  optimizes  SNR.  With  the 
amplitude  compression  method,  the  results  of  the  one-way 
experiments  show  that  the  mainlobe  shape  can  be  recovered 
down  to  —30  dB. 

L  Introduction 

ULTRASONIC  WAVEFRONT  DISTORTION  which  includes 
phase  aberration  and  waveform  distortion  develops 
as  coherent  waves  propagate  through  an  inhomogeneous 
medium,  such  as  the  breast.  Wavefront  distortion  sources 
inside  the  breast  can  be  modeled  as  incoherent  scattering 
and  coherent  multipath  interference.  Scattering  reduces 
the  target  strength,  broadens  the  image  lobe  [1],[2],  raises 
the  background  level  and  therefore  lowers  the  image  con¬ 
trast  [3j.  Multipath  interference  produces  ghost  image  ar¬ 
tifacts  or  false  targets  in  addition  to  true  targets  in  the 
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image  [4].  The  interference  problem  is  more  severe  when 
the  aperture  is  large  for  high  spatial  resolution. 

Adaptive  compensation  of  wavefront  distortion  is  an 
active  area  of  current  research  [5], [6].  Two  basic  phase 
deaberration  algorithms  have  been  developed  in  [7]-[9]. 
The  algorithms  adjust  time  delays  of  received  waveforms 
at  different  elements  of  a  phased  array  to  compensate  de¬ 
lay  errors  due  to  an  inhomogeneous  medium  and  conse¬ 
quently  reduce  the  scattered  energy.  Although  these  two 
algorithms  have  been  derived  from  different  optimization 
criteria,  their  performances  are  similar;  they  minimize  time 
delay  errors  or  pha.se  errors  but  leave  waveform  distortion 
intact  [10].  Nevertheless,  phase  deaberration  algorithms 
are  useful  to  partially  remove  phase  distortion  and  build 
up  the  strength  of  the  coherent  field.  Waveform  distor¬ 
tion  remains,  as  does  residual  phase  distortion,  and  to¬ 
gether  these  produce  a  significantly  high  background  level 
in  the  image. 

The  compensation  of  wavefront  distortion  presents  a 
more  difficult  problem,  because  there  is  no  general  way 
to  reduce  waveform  distortion.  Liu  and  Waag  [11]  have 
improved  the  time-delay  type  compensation  algorithm  by 
backpropagating  the  received  waveforms  to  an  optimal 
distance  and  then  performing  time-delay  compensation 
at  this  distance.  Backpropagation  partially  compensates 
waveform  distortion  produced  by  wave  propagation  from 
the  backpropagation  distance  to  the  receiving  aperture. 
We  have  developed  a  novel  toward  inverse  filtering  ap¬ 
proach  (TIF)  which  is  a  spectrum- amplitude  compression 
technique  that,  when  used  in  combination  with  phase  de- 
aberration  and  backpropagation,  provides  a  general  means 
of  reducing  waveform  distortion  and  therefore  reducing  the 
incoherent  background  level  in  the  image. 

II.  Toward  Inverse  Filtering  Approach  (TIF) 

The  following  discussions  use  the  correction  of  a  point 
source  image  as  a  means  to  compare  different  correction 
methods  and  to  introduce  TIF.  But  TIF  is  not  limited  to 
point  source  targets  and  is  applicable  to  a  diffuse  scattering 
medium  (see  Section  V,  D),  such  as  the  breast.  In  the  re¬ 
maining  text,  wavefront  amplitude  and  phase  refer  specif¬ 
ically  to  the  complex  spectral  amplitude  at  each  discrete 
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DEABERRATION  TRANSFORMATIONS 


MF  =  matched  filter  (SNR) 

q5A  =  dominant  scuttcrcr  algorithm 

Y 

f 

S^^^UNIT  CIRCLE 

TRM  =  time  reversal  minrior^.^-'''' 

PC  =  phase  conjugation 

TDC  =time-delay  compjzmation 

BPT=  Backpropagatio/and 
time-delay  compcns^ion 

1 

SR  =  square  rooter  1 
(compromise  ampliiirde  weight) 

[F  =  inverse  filter  (fid^ty) 

DSATDCt  1 

TRM  BPT  / 

X 

o  complex  sampl 

^measured  in  array  f=  aexpCj  W 

•  following  phase  correction 

Fig  1  Deaberration  transformations.  /  is  distortion  vector  for  in¬ 
stantaneous  sample  of  the  radiation  field.  Transformations  to  the 
real  axis  all  correct  phase.  Left  (MF)  maximizes  SNR  but  increases 
amplitude  distortion.  Right  (IF)  maximizes  imaging  fidelity  but  m- 
creases  noise.  SR  is  good  compromise. 

frequency  within  the  signal  band.  The  correction  method 
is  developed  in  frequency  domain;  it  allows  for  frequency 
dependent  changes  in  both  phase  and  amplitude.  TIF  com¬ 
bines  a  spectrum  amplitude  compression  operation  with 
phase  deaberration.  ,  , 

Consider  complex  sigiiai  vectors  So(x)  =  Aoexpij'ipo) 
and  S(x)  =  Aexp(j'0)  received  across  an  array  in  the  ab- 
sence  and  presence  of  distortion,  respectively.  The  ratio  of 
S  to  So  is  a  multiplicative  distortion  vector  f  =  aexp{j(p). 
Fig.  1  shows  such  a  distortion  vector  f  for  some  instanta¬ 
neous  sample  of  the  radiation  field.  The  optimum  compen¬ 
sation  transformation  is  the  right-most  one,  which  carries 
the  distortion  component  of  the  complex  sample  to  the  in¬ 
tersection  of  the  x-axis  and  the  unit  circle,  for  then  both 
amplitude  and  phase  are  corrected.  The  left-most  trans¬ 
formation  is  matched  filtering  (MF),  which  squares  the 
amplitude  and  conjugates  the  phase.  MF  is  theoretically 
optimum  for  maximizing  SNR  on  a  white,  Gaussian  chan¬ 
nel.  It  is  an  optimum  solution  for  detection  (radar)  but  not 
so  when  fidelity  is  an  important  criterion,  as  in  imaging. 
But  it  does  satisfy  the  requirement  for  phase  correction. 
Examples  of  MF  algorithms  are  dominant  scatterer  algo¬ 
rithm  (DSA)  [12]  and  time  reversal  mirror  (TRM)  [13]. 
We  should  expect  that  MF  algorithms  are  useful  when 
wavefront  distortion  resides  principally  in  phase,  with  am¬ 
plitude  relatively  unaffected.  Otherwise  MF  algorithms 
would  be  less  effective  than  if  amplitude  were  ignored  and 
phase  only  were  corrected.  By  squaring  the  wavefront  am¬ 
plitude,  the  MF  algorithms  increase  amplitude  variance 
which  contributes  to  the  increase  of  energy  in  the  sidelobe 


Section  V,  A).  The  signal  amplitude  is  unchanged  and  the 
phcise  error  is  completely  compensated  by  conjugation. 
This  operator  is  represented  by  the  2nd  (from  the  left) 
transformation.  Time  delay  compensation  (TDC)  [7]-[9| 
is  in  the  same  class. 

The  right-most  transformation  which  corrects  both  am¬ 
plitude  and  phase  is  the  inverse  filter  (IF).  It  is  theoret¬ 
ically  ideal  for  fidelity  but  has  drawbacks.  Because  the 
IF  adjusts  the  channel  gain  to  be  the  reciprocal  ot  the 
signal  strength,  at  points  in  the  receiving  aperture  where 
signal  -Strength  is  weak  the  enhanced  channel  gain  raises 
the  noise  to  the  point  where  SNR  can  be  impaired.  This 
is  particularly  troublesome  when  there  is  coherent  refrac¬ 
tive  interference  in  the  receiving  array.  A  second  potential 
problem,  which  we  have  avoided  by  not  working  too  close 
to  the  IF  point,  is  that  an  IF  is  unstable  when  the  dis¬ 
torting  medium  has  zeroes  in  the  complex  plane.  Thus  the 
right-most  transformation  carries  the  correction  too  far. 

The  complex  weights  of  MF,  PC,  and  IF  are,  respec¬ 
tively,  a,  1,  and  a- \  each  multiplied  by  exp( -;</>).  The 
signal  after  weighting  is  a^,  a,  and  1.  The  first  and  last 
are  far  from  optimum  for  the  reasons  given.  The  4th  from 
the  left  transformation  SR  is  a  compromise.  The  correc¬ 
tion  vector  is  exp(-j<^).  The  signal  after  weight¬ 

ing  is  (called  later  rooter).  Another  transformation 
to  the  right  of  the  rooter  is  fourth  rooter  (signal  am¬ 
plitude  is  0^7“!  after  weighting).  In  general,  the  correc¬ 
tion  vector  is  and  signal  amplitude 

is  after  transformation  of  the  Mth  rooter;  this  oper¬ 

ation  is  a  toward  inverse  filtering  (TIF)  approach.  In  Sec¬ 
tion  IV,  we  demonstrate  that  low  order  rooters  (rooter  and 
fourth  rooter)  when  combined  with  phase  deaberration  at 
the  aperture  can  significantly  improve  contrast  resolution 
without  encountering  the  stability  problem. 

Backpropagating  the  received  wavefront  to  an  optimal 
backpropagation' distance  and  performing  phase  deaber¬ 
ration  at  this  optimal  distance  (BPT)  [11]  provides  first 
order  correction  of  wavefront  amplitude  distortion  due  to 
propagation  from  the  backpropagation  distance  to  the  re¬ 
ceiving  aperture  and  therefore  performs  better  than  phase 
deaberration  at  the  aperture.  In  Section  IV,  we  demon¬ 
strate  experimentally  that  BPT  transformation  is  some¬ 
where  to  the  right  of  TDC  and  PC,  and  left  of  PC  and 
rooter.  We  also  demonstrate  that  square  rooter  in  addition 
to  BPT  can  significantly  improve  the  performance  of  BPT. 

The  rooter  transformation  upon  wavefront  amplitude  is 
a  nonlinear  compression  operation.  Other  nonlinear  am¬ 
plitude  compression  functions  which  converge  toward  the 
inverse  filtering  solution  may  also  be  used  in  place  of  the 
Mth  rooter. 


region  [3]. 

A  phantom  consisting  of  a  random  phase  screen  placed, 
at  a  receiving  transducer  models  the  situation  where  the 
element  signal  samples  are  correct  in  amplitude  and  are 
distorted  only  in  phase  by  a  random  additive  component. 
Phase  aberration  correction  (PAC),  consisting  of  phase 
conjugation  (PC),  is  then  the  optimum  compensation  (see 


III.  Experiments 

A.  In  Vitro  Measurement 

The  2-D  wavefront  measurement  system  and  procedure 
were  described  in  detail  in  [14]  and  are  summarized  here 
for  convenience.  A  breast  tissue  specimen  was  placed  in 
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the  specimen  holder  with  the  skin  facing  the  direction  of 
the  receiving  transducer.  The  source  was  a  hemispheri¬ 
cal  transducer  emulating  a  virtual  point  source.  The  2-D 
array  consisted  of  a  92-mm  1-D  linear  array  translated 
46  mm  perpendicular  to  its  axis  to  form  a  synthetic  2- 
D  array  92  mm  x  46  mm.  The  element  pitch  in  the  re¬ 
ceiving  transducer  was  0.72  mm  and  a  reflecting  mask  re¬ 
duced  the  receiving  elevation  to  1.44  mm.  Transmission 
WcLS  one-way.  Frequency  was  3.7  MHz.  Bandwidth  was 
2  MHz.  Waveforms  were  measured  at  each  element,  from 
which  2-D  wavefronts  were  reconstructed  as  functions  of 
time. 

Breast  tissue  specimens  were  obtained  fresh  from  reduc¬ 
tion  mammoplasty  surgery  and  were  stored  frozen;  they 
were  later  defrosted  and  degassed  for  the  experiments. 
Each  specimen  was  essentially  planar  and  had  a  surface 
area  of  at  least  7  x  11  cm-.  The  average  thickness  was 
26.9  mm.  The  tissue  donors  were  women  ranging  in  age 
from  18  to  65  with  a  mean  age  of  34  years. 

Three  groups  of  measurements  were  made.  The  primary 
set  consisted  of  independent  measurements  using  nine  dif¬ 
ferent  specimens.  In  addition,  two  specimens  were  em¬ 
ployed  for  sequential  measurements  in  which  an  approx¬ 
imately  1-cm  layer  was  removed  from  the  bottom  of  the 
specimen  before  each  subsequent  measurement.  For  three 
other  specimens,  a  pair  of; measurements  was  made  with 
the  source  in  each  of  two  positions  located  12  mm  apart 
in  the  array  direction.  A  total  of  16  measurements  were 
made.  Water  measurement  was  also  made  for  compari¬ 
son.  In  the  following  analysis,  the  sequential  measurements 
were  treated  as  measurements  from  different  samples  be¬ 
cause  the  samples  resulting  from  approximately  1-cm  re¬ 
moval  from  the  original  samples  were  replaced  into  the 
sample  holder  after  cutting,  and  their  relative  positions  to 
the  2-D  transducer  aperture  changed  unpredictably.  The 
pair  measurements  resulting  from  12  mm  lateral  source 
translation  were  also  treated  as  measurements  from  dif¬ 
ferent  samples  because  the  measured  parameters  and  the 
point-source  image  properties  were  changed  due  to  the 
highly  inhomogeneous  nature  of  the  breast  tissue. 

B.  Data  Processing 

Pulse  arrival  time  at  each  element  was  calculated  by 
using  a  differential  edge  detector  to  locate  the  first  nega¬ 
tive  peak  of  the  pulse.  Geometric  effects  were  removed  by 
fitting  a  two-dimensional,  fourth-order  polynomial  to  the 
calculated  arrival  time  surface  and  subtracting  the  result. 
The  reasons  for  using  fourth-order  polynomial  fitting  were 
given  in  [11].  3-D  image  data  were  calculated  by  Fourier 
decomposition  of  the  temporal  waveforms  at  each  {x,y) 
position  in  the  aperture,  calculation  of  a  complex  CW  2-D 
image  at  each  Fourier  frequency  in  the  focal  plane  by  us¬ 
ing  the  angular  spectrum  technique  [15],  and  summation 
of  the  2-D  images  to  form  the  3-D  transient  image.  The 
final  2-D  image  used  in  this  report  was  obtained  by  detect¬ 
ing  the  peak  pressure  value  at  each  {xf,  yf)  position  in  the 
image  plane  within  the  transient  period.  Contour  maps  of 
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2-D  point  source  images  at  different  thresholds  are  used  to 
quantitatively  evaluate  image  quality. 

Before  Fourier  decomposition,  the  received  pulse  at  each 
element  was  temporally  weighted  by  a  10-point  cosine  ta¬ 
per  over  the  trailing  10  points  of  a  40-point  time  window 
that  contained  the  main  arrival  pulse.  The  original  40- 
point  interval  was  padded  with  zeros  to  a  128-point  inter¬ 
val. 

In  step  one,  a  frequency  range  (0.31-6.25  MHz)  that 
contained  most  of  the  signal  energy  was  chosen.  This  range 
corresponds  to  the  3rd  through  41st  harmonic  of  the  FFT 
of  a  signal  sampled  at  20  MHz  over  a  6.4  /.is  interval  con¬ 
taining  128  points. 

In  step  two,  the  complex  signals  in  the  128  x  32  spa¬ 
tial  aperture  were  zero  padded  symmetrically  to  256  x  64 
before  each  harmonic  component  in  the  2-D  aperture  was 
focused  in  the  image  plane  via  the  angular  spectrum  tech¬ 
nique.  Spatial  raised  cosine  taper  was  employed.  The  sys¬ 
tem  bandwidth  is  about  50%;  therefore,  Q  =  ^  2, 

where  /o  is  the  central  frequency.  For  such  Q,  the  first 
sidelobes  (due  to  finite  aperture  size)  in  both  principal 
axes  remain  [16]  and  are  reduced  to  about  —30  dB  level 
by  employing  the  raised  cosine  taper  [17]. 

C.  Implementation  of  Correction  Methods 

Two  basic  phase  deaberration  algorithms,  time-delay 
compensation  (TDC)  and  phase  conjugation  (PC),  and 
two  wavefront  compensation  algorithms,  wideband  DSA 
and  BPT,  were  applied  to  each  measured  2-D  wavefront 
set  and  their  effects  upon  focusing  were  compared.  TIF 
was  implemented  by  taking  rooter  and  fourth  rooter  on 
wavefront  amplitude  in  combination  with  PC,  BPT,  and 
TDC;  its  effects  upon  focusing  were  compared  with  the 
other  four  deaberration  algorithms. 

TDC:  An  arrival  time  fluctuation  profile  at  the  aperture 
is  calculated  in  three  steps. 

Step  1 — A  reference  waveform  with  good  signal- to- noise 
ratio  is  first  selected  from  geometrically  corrected  wave¬ 
forms.  The  selected  waveform  is  continuously  modified  to 
reduce  the  effect  of  the  arbitrariness  in  the  initial  selec¬ 
tion  as  waveforms  in  the  aperture  are  cross-correlated  with 
it.  The  modification  is  accomplished  by  time-shifting  and 
adding  a  new  waveform  to  the  reference  pulse  if  the  peak 
value  of  the  correlation  function  between  the  new  wave¬ 
form  and  the  current  reference  waveform  is  greater  than 
0.8.  The  average  total  number  of  waveforms  that  are  in¬ 
corporated  into  the  final  reference  waveform  selection  is 
2690  (std.  391),  which  accounts  for  66%  (std.  9.5%)  of  the 
total  pulses.  The  average  is  calculated  upon  a  total  of  16 
in  vitro  samples.  This  selection  procedure  is  similar  to  the 
one  developed  in  [11]. 

Step  2 — The  arrival  time  fluctuation  profile  is  first  cal¬ 
culated  from  cross-correlation  peaks  of  geometrically  cor¬ 
rected  waveforms  and  the  reference  waveform  when  the 
peak  values  of  the  correlation  function  between  the  wave¬ 
forms  and  the  reference  waveform  are  greater  than  0.7.  At 
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this  step,  the  average  total  number  of  waveforms  that  are 
incorporated  into  the  calculation  of  the  arrival  time  profile 
is  3518  (std.  173),  which  accounts  for  86%  (std.  4.2%)  of 
the  total  pulses. 

Step  3 — For  the  bad  waveforms  with  cross-correlation 
peaks  less  than  0.7,  the  arrival  time  fluctuations  at  these 
positions  are  calculated  as  follows:  one,  calculate  cross¬ 
correlation  peaks  between  each  bad  waveform  and  neigh¬ 
bor  waveforms  (within  a  5  by  5  spatial  window)  that  are 
incorporated  in  Step  2;  two,  select  those  neighbor  wave¬ 
forms  having  cross-correlation  peaks,  calculated  in  one, 
greater  than  0.5;  three,  calculate  the  average  of  arrival 
time  fluctuation  from  arrival  time  differences  between  the 
bad  waveform  and  the  selected  neighbor  waveforms,  and 
the  arrival  time  fluctuations  of  the  corresponding  neigh¬ 
bor  positions.  This  average  is  assigned  to  the  bad  wave¬ 
form  position.  A  justification  of  this  step  is  that,  for  badly 
distorted  waveforms,  they  are  dissimilar  to  the  reference 
waveform  to  a  large  extent,  but  they  are  somewhat  similar 
to  the  waveforms  in  the  closest  neighborhood  which  are 
reasonably  correlated  with  the  reference  waveform.  After 
this  step,  the  average  total  number  of  pulses  that  are  in¬ 
corporated  into  the  calculation  of  the  arrival  time  profile 
is  3966  (std.  41),  which  accounts  for  97%  (std.  1%)  of  the 
total  pulses. 

For  the  waveforms  that  have  not  been  incorporated  into 
the  calculation  of  the  arfwal  time  profile  after  Step  3,  ar¬ 
rival  time  fluctuations  at  these  positions  are  assigned  to 
zero.  The  calculated  average  arrival  time  fluctuation  of 
the  entire  16  profiles  is  37,6  ns  with  a  10.5  ns  standard 
deviation. 

Except  for  TDC,  the  following  correction  procedures 
were  implemented  upon  the  complex  wavefront  at  each 
frequency  after  Fourier  decomposition.  The  image  forma¬ 
tion  procedures  after  using  each  correction  method  were 
the  same  as  the  procedures  described  in  part  B  steps  2 
and  3. 

PC:  Let  Af{x,y)  exp{j(l)f{x,  y))  represent  the  complex 
signal  at  each  frequency  after  Fourier  decomposition, 
where  Af  and  0/  are  functions  of  positions  {x,y)  in  the 
2-D  array  at  frequency  /.  The  complex  weight  vector  of 
phase  conjugation  Is  Wf  ^  exp(— j0/)  and  the  signal  after 
weighting  is  A/.  The  image  formation  after  phase  conju¬ 
gation  is  the  same  as  described  in  part  B,  steps  2  and 
3.  Phase  conjugation  at  each  frequency  within  the  signal 
band  (without  amplitude  weight)  provides  perfect  phase 
compensation  except  for  a  27r  phase  jump.  The  result  is 
better  than  TDC  because  of  the  removal  of  residual  phase 
error.  The  wavefront  amplitude  distortion  remains,  how¬ 
ever.  The  objective  here  is  to  evaluate  the  limit  of  phase 
compensation  and  the  extent  of  the  wavefront  amplitude 
distortion  effect. 

Wideband  DSA:  The  complex  weight  vector  of  DSA  at 
each  frequency  is  Wf  =  A/  Qxp{—j(j)f)  and  the  signal  after 
weighting  is  Aj. 

Wideband  DSA  is  an  approximation  of  the  TRM  receive 
mode  developed  in  [18].  The  difference  between  wideband 
DSA  and  TRM  receive  mode  is  the  scanning  procedure. 
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Fig.  2.  Scanning  configuration. 


At  target  position  A  (Fig.  2),  wideband  DSA  and  TRM 
receive  mode  produce  exactly  the  same  responses,  because 
of  the  reciprocity  principle.  At  nontarget  position  B,  wide¬ 
band  DSA  uses  linear  geometry  between  A  and  B  to  scan 
the  beam  while  the  TRM  receive  mode  could  backprop- 
agate  the  wave  to  the  aberrator  position,  and  then  use 
linear  geometry  to  scan  the  beam  [18],  However,  for  points 
that  are  within  the  correlation  distance  of  A,  wideband 
DSA  and  the  TRM  receive  mode  should  produce  similar 
responses.  The  correlation  distance  of  the  breast  is  about 
4  to  5  mm  [14],  [19]-[21].  Therefore  within  a  4  to  5  mm  ra¬ 
dius,  we  should  expect  similar  performances  of  wideband 
DSA  and  TRM  receive  mode.  The  objective  of  using  wide¬ 
band  DSA  here  is  to  evaluate  the  possible  improvement  of 
deaberration  due  to  the  additional  SNR  gain  provided  by 
amplitude  weighting. 

BPT:  BPT,  developed  in  [11],  is  implemented  by  back- 
propagating  the  complex  wavefront  at  each  frequency  after 
Fourier  decomposition  to  the  optimal  distance  and  then  us¬ 
ing  phase  conjugation  at  that  distance.  The  use  of  phase 
conjugation  instead  of  TDC  at  the  backpropagation  dis¬ 
tance  reduces  the  complexity  in  the  calculation.  The  back- 
propagation  distances  are  taken  from  [14]. 

TIP: 

PC  and  rooter  (PC-rooter) — The  weight  vector  of  PC 

-I/O 

and  rooter  at  each  frequency  is  Wf  =  Aj^  '  exp(-j(p/) 
and  the  signal  after  weighting  is  Aj  " . 

PC  and  fourth  rooter  (PC-fourth  rooter) — The  weight 
vector  at  each  frequency  is  u;/  = 
signal  after  weighting  is 

BPT  and  rooter  (BPT-rooter) — BPT  and  rooter  was 
implemented  the  same  as  PC  and  rooter,  but  the  procedure 
was  done  at  the  optimal  back  propagation  distance  after 
backpropagating  each  complex  wavefront  to  the  optimal 
distance. 

TDC  and  rooter  (TDC-rooter) — TDC  and  rooter  was 
implemented  by  correcting  delay  errors  in  time  domain 
using  TDC,  and  taking  rooter  on  wavefront  amplitude  at 
each  frequency. 
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IV.  Experimental  Results 
A.  Correction  of  Scattered  Energy 

1.  Coirection  Upon  Scattering  Samples  Figs.  3(a)  and 
(b)  show  -6  dB  contour  plots  of  images  obtained  from  the 
water  path  and  the  3.5-cm  breast  tissue  path  (brsOOS). 
Contour  spacing  is  1  dB.  Abscissa  and  ordinate  are  az¬ 
imuth  and  elevation  in  millimeters  in  the  image  plane. 
Fig.  3(a)  shows  the  system  diffraction  pattern  at  -6  dB. 
The  azimuthal  beamwidth  is  about  1.5  mm  (at  a  distance 
of  180  mm),  meaning  that  the  point  (or  lateral)  resolution 
is  1.5/180  or  8.3  mrad.  Because  the  size  of  the  receiving 
array  in  elevation  is  half  the  size  in  the  array  direction, 
the  width  of  the  image  in  elevation  is  twice  as  large  as  in 
azimuth.  When  tissue  is  present  Fig.  3(b),  the  image  lobe 
is  broadened.  In  Fig.  3(b)  the  beamwidth  has  grown  from 
1.5  mm  to  3  mm  in  azimuth  and  3.1  mm  to  9.3  mm  in 
elevation  and  the  point  resolution  has  worsened  approxi¬ 
mately  2:1  in  azimuth  and  3:1  in  elevation.  A  symmetrical 
scattering  pattern  appears  when  the  threshold  is  reduced. 
Fig.  3(c)  shows  a  -20  dB  contour  plot  of  the  sample  image 
brs005.  Contour  spacing  is  2  dB.  Energy  is  spread  out  over 
a  large  area.  No  distiiict  structure  or  multipath  lobe  can 
be  identified  from  the  irnage. 

Fig.  4(a)-(g)  shows  the  correction  result  of  the  sample 
image  brs005  by  using  different  methods.  Part  (h)  is  the 
-30  dB  level  contour  of  a  water  path  image.  The  contour 
spacings  in  Fig.  4  are  the  outer  contour  levels  divided  by 
10.  Fig.  4(a)-(c)  shows  TDC  at  —6,  —20,  and  —30  dB  lev¬ 
els  (see  Fig.  3(b)  and  (c)  for  comparison).  Image  quality 
improved  significantly  and  the  system  diffraction  pattern 
is  restored  up  to  -20  dB  level.  The  symmetric  scattering 
pattern  outlined  by  outer  contours  at  the  -30  dB  level  is 
caused  by  residual  phase  errors  after  phase  deaberration 
and  incoherent  amplitude  distortion.  Fig.  4(d)  is  the  cor¬ 
rection  result  of  PC  at  the  -30  dB  level.  Because  of  the 
complete  removal  of  phase  errors  except  for  a  27r  jump, 
the  scattering  pattern  is  further  improved  as  compared  to 
Fig.  4(c).  Fig.  4(e)  is  the  result  of  BPT,  and  the  scattering 
pattern  is  much  better  than  that  of  PC  at  the  aperture. 
The  optimum  backpropagation  distance  of  this  tissue  sam¬ 
ple  is  30  mm.  Fig.  4(f)  is  the  correction  result  of  PC-rooter 
at  -30  dB.  The  pattern  is  very  close  to  the  system  diffrac¬ 
tion  pattern  (see  (h)  for  comparison).  Fig.  4(g)  is  the  result 
of  wideband  DSA  at  the  -30  dB  level.  The  correction  re¬ 
sult  is  worse  than  those  of  other  methods  (see  Section  V , 

C). 

2.  Correction  Upon  Samples  with  Well-defined  Multi- 
path  In  most  breast  samples,  scattering  and  coherent  in¬ 
terference  are  both  present.  Fig.  .5(a)  and  (b)  are  —10  dB 
contour  plots  of  images  obtained  from  the  water  path  and 
the  4-cm  breast  tissue  path  (brsOOO).  Contour  spacing  is 
2  dB-.  In  tissue  image  Fig.  5(b),  three  lobes  instead  of  one 
central  image  lobe  appear.  Inner  contours  of  the  image 
lobe  are  close  to  the  system  diffraction  pattern  (Fig.  5(a)) 
while  the  outer  contours  are  highly  irregular  due  to  scat¬ 
tering.  The  point  resolution  is  worsened  by  approximately 


Fig.  3.  2-D  contour  maps  of  point  source  images,  (a)  Water  data 
showing  diffraction  pattern  of  system.  Outer  contour  is  -6  dB  level, 
(b)  Measured  through  the  3.5-cm  tissue  (brsOOS),  showing  asymmet¬ 
ric  scattering  effect.  Outer  contour  is  -6  dB  level,  (c)  Image  brsOOS 
at  -20  dB  contour  level.  The  outer  contour  shows  symmetrical  scat¬ 
tering  pattern.  No  distinct  structure  or  multipath  is  observed.  Note 
the  change  of  scale  from  (a)  and  (b). 
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(b)  -lOdB  (brs006) 


(c)  -16  dB  (brs006) 


Fig.  5.  2-D  contour  maps  of  point  source  images,  (a)  Water  data 
showing  diffraction  pattern  of  system.  Outer  contour  is  -10  dB  level, 
(b)  Measured  through  the  4-cm  tissue  (brs006),  showing  highly  asym¬ 
metric  interference  pattern.  Outer  contour  is  -10  dB  level,  (c)  Image 
brsOOe  at  -16  dB  contour  level.  The  outer  contour  shows  more  sym¬ 
metrical  scattering  pattern.  Note  the  change  of  scale  from  (a)  and 

(b). 


Fig.  6.  (a)  45  degree  cuts  of  2-D  images;  water,  (b)  brs006. 

2:1  in  azimuth  and  3:1  in  elevation.  Two  refractive  lobes 
appear  at  (-4,  ~2.5)  and  (3.5, 1.5),  and  their  strengths  are 
—6  and  ~8  dB.  A  closer  examination  is  obtained  from  45 
degree  cuts  (Fig.  6)  of  the  Fig.  5(a)  and  (b)  images.  In 
Fig.  6(a),  the  peak  sidelobe  is  -36  dB,  which  is  the  best 
we  can  obtain  from  the  measurement  system.  In  the  tis¬ 
sue  image,  two  large  sidelobes  (— 6  dB  and  —  8  dB)  appear 
on  either  side  of  the  mainlobe.  The  two  nulls  between  the 
lobes  are  deep  (left:  -13  dB  below  peak,  right:  -15  dB), 
which  is  a  typical  coherent  interference  phenomenon. 

A  symmetric  scattering  pattern  generally  appears  when 
the  threshold  is  reduced,  as  shown  Fig.  6(c)  with  a  -16  dB 
threshold  contour  and  contour  spacings  of  approximately 
1.5  dB.  The  outer  contours  show  a  roughly  symmetric  pat¬ 
tern  typical  of  scattering  while  the  inner  contours  show  an 
asymmetric  lobular  interference  pattern. 

Fig.  7  shows  correction  results  of  the  sample  image 
brs006  after  applying  different  methods.  The  contour  spac¬ 
ings  are  the  outer  contour  levels  divided  by  10.  Fig.  7(a)  is 
the  -16  dB  contour  map  of  applying  TDC  (see  Fig.  5(c)  for 
comparison).  Image  quality  improved  significantly  because 
of  the  minimization  of  the  phasefront  distortion  caused  by 


elevation  (mm) 
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scattering.  The  area  within  the  —16  clB  inner  contour  is 
reduced  by  a  factor  of  11.  The  mainlobe  pattern  is  close 
to  the  system  diffraction  pattern,  although  the  area  is  still 
1.6  times  larger  than  the  area  in  water  at  the  same  level. 
Two  large  lobes  remain  as  expected  and  appear  as  two 
additional  sources.  Fig.  7(h)  is  the  correction  result  at  a 
-16  dB  level  of  applying  wideband  DSA.  It  is  evident  that 
the  procedure  enhances  interference  in  this  example;  the 
artifacts  are  5  dB  higher  than  in  Fig.  7(a).  The  reason  is 
given  in  [3].  Figs.  7(b)  and  (c)  show  -25  and  -30  dB  con¬ 
tour  plots  of  TDC.  In  Fig.  7(b)  the  outer  contours  show 
an  asymmetric  interference  pattern  with  scattered  energy 
distributed  around  the  two  interference  lobes.  In  Fig.  7(c) 
the  outer  contours  show  the  more  symmetric  scattering 
pattern  caused  by  incoherent  wavefront  amplitude  distor¬ 
tion  and  residual  phase  distortion.  Fig.  7(d)  is  the  result  of 
PC  at  -30  dB.  The  approximately  symmetric  scattering 
pattern  outlined  by  outer  contours  in  (c)  is  improved  due 
to  better  phase  deaberration.  The  energy  is  more  concen¬ 
trated  around  the  image  lobe  and  interference  lobes.  Part 
(e)  is  the  result  of  BPT,  and  the  scattering  pattern  is  fur¬ 
ther  improved  as  compared  with  the  result  of  PC  at  the 
aperture.  The  backpropagation  distance  of  this  sample  is 
50  mm.  Part  (f)  is  the  result  of  PC-rooter  at  -30  dB.  The 
pattern  is  very  close  to  the  system  diffraction  pattern.  PC- 
fourth  rooter  shown  in  Fig.  7(g)  restores  the  diffraction- 
limited  image  to  a  -30  dB  level. 

It  is  interesting  to  note  that  the  strengths  of  the  two 
interference  lobes  are  also  reduced  by  8  dB  after  TDC 
and  PC.  They  are  further  reduced  by  5,  7,  and  15  dB 
after  BPT,  PC-rooter,  and  PC-fourth  rooter,  respectively. 
The  improvement  comes  from  the  scattered  energy  which 
is  originally  distributed  around  the  mainlobe  and  each  of 
the  refracted  multipath  lobes.  This  energy  is  coherently 
added  back,  by  the  deaberration  process,  to  the  dominant 
lobe,  which  in  this  case  is  the  mainlobe.  The  overall  result 
is  an  improvement  of  the  image  lobe  to  the  interference 
lobe  ratio. 

Determination  of  whether  interference  was  coherent  or 
scattered  was  based  upon  whether  or  not  large  lobes  could 
be  removed  by  applying  phase-deaberration  algorithms. 
Another  indicator  is  whether  or  not  a  lobe  location  moves 
(scattering)  or  remains  relatively  the  same  (refraction)  af¬ 
ter  the  phase  deaberration  process.  Both  means  were  used 
to  identify  whether  these  lobes  were  produced  by  refrac¬ 
tion  or  scattering.  Large  lobes  in  nine  samples  are  believed 
to  be  caused  by  the  interference  process.  In  these  samples, 
we  have  found  an  average  of  a  10  dB  (std.  4dB)  improve¬ 
ment  in  the  range  of  5  to  16  dB  after  TDC,  due  to  the 
coherent  strengthening  of  the  mainlobe. 

B.  Statistics 

1.  Point  Resolution  and  Refractive  Artifacts  We  ana¬ 
lyzed  images  at  different  contour  levels  and  classified  the 
total  of  16  samples  into  three  categories:  primarily  scatter 
(4  samples),  more  scatter  than  refraction  (3  samples),  and 
scatter  plus  well-defined  refraction  (9  samples).  Statistics 


TABLE [ 

Statistics  of  Increase  in  -6  dB  Beam  width  (Percent). 


Sample  size 

Classification 

Azimuth 

Elevation 

4 

A 

130  (std.  131.3) 

100  (std.  96.5) 

3 

B 

•50  (std.  .53.3) 

.30  (std.  .30.0) 

9 

C 

.30  (std,  33.3) 

40  (std.  42.2) 

A.  Primarily  scatter 

B.  More  scatter  than  refraction 

C.  Scatter  plus  well-defined  refraction 


TABLE  n 

Statistics  of  Refracted  Lobes  From  9  M  Vitro 
Breast  Samples. 


Average 

Std.  Dev. 

Number  of  refractive  lobes 

2 

0.53 

Relative  peak  strength  (dB) 

10 

4.7 

Distance  from  primary  lobe  (deg) 

1.8 

0.57 

of  -6  dB  point  resolution  of  three  groups  is  given  in  Table 
I.  The  average  beam-broadening  as  compared  with  water 
is  70%  in  azimuth  and  57%  in  elevation.  These  numbers 
are  larger  than  the  result  reported  in  [21]  where  an  18% 
increase  in  beamwidth  at  the  same  threshold  is  observed. 
This  is  likely  due  to  the  different  aperture  sizes  used.  For 
the  samples  with  refractive  artifacts,  the  point  resolution 
is  evaluated  upon  the  central  image  lobe.  The  percentage 
increase  in  beamwidth  of  the  scatter  group  is  significantly 
larger  than  that  of  the  other  two  groups.  It  is  obvious  that 
point  resolution  is  not  enough  to  evaluate  image  distor¬ 
tion  when  multiple  lobes  are  present.  Statistics  of  refrac¬ 
tive  lobes  of  the  third  group  is  given  in  Table  II.  In  this 
group,  we  found  1-3  refractive  lobes  at  levels  between  —4 
and  -17  dB  (peak  value),  all  within  a  radius  of  1-3°. 

2.  A  Measure  of  Contrast  Contrast  is  the  measure  used 
to  quantitatively  evaluate  improvement  upon  image  qual¬ 
ity  through  the  use  of  different  correction  algorithms.  Im¬ 
age  contrast  is  the  metric  of  the  brightness  point  in  an 
image  to  the  peak  background,  i.e.,  it  defines  the  dynamic 
range  window  between  target  strength  and  background 
(noise,  scatter,  refractive  artifacts).  It  is  similar  to  sidelobe 
level  as  a  measure  of  quality  of  a  radiation  pattern  [17] . 
In  a  point  source  image,  the  measurement  can  be  made 
by  direct  observation  of  the  background  peak.  However, 
in  most  cases  the  background  peak  occurs  on  the  skirts 
of  the  main  diffraction  lobe  of  the  target  and  its  ampli¬ 
tude  cannot  be  directly  discerned.  Its  effect  is  to  distort 
the  mainlobe.  In  these  cases,  contrast  is  measured  down 
to  the  lowest  level  at  which  mainlobe  diffraction  shape  is 
maintained.  The  mainlobe  diffraction  shape  is  identified 
when  the  elevation- image  width  to  the  horizontal  width 
is  approximately  a  2  to  1  ratio,  and  the  outermost  image 
contour  has  an  elliptical  shape.  The  evaluation  is  done  by 
visual  inspection  at  each  contour  level  until  the  criteria  are 
met.  At  the  levels  selected,  65%  of  ratios  are  within  0.1  of 
2  to  1,  27%  within  0.2  of  2  to  1,  and  8%  within  0.3. 


584 


IEE0  TRANSACTIONS  ON  ULTRASONICS,  FERROELECTRICS,  AND  FREQUENCY 

TABLE  in 
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Diffraction  pattern  (dB) 
w/o  correction 

TDC 

PC 

Diffraction  pattern  after  correction  (dB) 

Wideband  BPT  PC-rooter  BPT-rooter 

PC- fourth  rooter 

A  -0(0) 

B  -.5.8(1.3) 

C  -6.1(3.5) 

-18.8(1.9) 
-19.7(2. .3) 
-19. .3(2. 7) 

-19.0(2.2) 

-19.7(3.5) 

-19.6(2.5) 

-14.8(2.6) 

-14.7(1.5) 

-14.8(2.5) 

-24.3(2.1) 

-23.7(0.6) 

-23.9(2.4) 

-25.3(3.0) 

-25.3(2:1) 

-25.4(2.5) 

-29.5(2.7) 

-29.3(2.3) 

-29.3(3.0) 

-.32.3(2.8) 

-.30.7(1.2) 

-30.7(1.4) 

The  average  levels  (std.  dev)  of  the  restored  mainlobe 
diffraction  shape  before  and  after  corrections  of  the  three 
groups  classified  earlier  are  given  in  Table  III.  The  main- 
lobe  shape  is  restored  up  to  -19.3  (2.4)  and  -19.5  dB 
(2.6  dB)  by  using  two  phase-deaberration  algorithms,  TDC 
and  PC,  respectively.  The  fact  that  the  correction  result 
of  TDC  is  very  close  to  that  of  PC  suggests  that  the  limit 
of  phase  correction  with  respect  to  the  selected  metric  has 
been  achieved  by  TDC.  Further  improvement  requires  cor¬ 
rection  algorithms  that  can  take  wavefront  amplitude  dis¬ 
tortion  into  account.  Wideband  DSA  is  able  to  restore  the 
mainlobe  shape  to  -14.8  dB  (2.3  dB).  This  result  is  about 
5  dB  worse  than  that  of  using  two  phase-deaberration  al¬ 
gorithms.  The  cause  is  the  increase  of  wavefront  amplitude 
variance  by  squaring  the  amplitude.  Wavefront  amplitude 
variance  contributes  to  the  energy  in  the  sidelobe  region 
[3j.  BPT  at  an  optimal  baickpropagation  distance  improves 
the  performance  of  PC  at  aperture  by  5  dB.  PC-rooter  re¬ 
stores  the  mainlobe  diffraction  shape  to  -25.4  dB  (2.6  dB). 
The  correction  result  of  BPT  is  similar  to  PC-rooter.  PC- 
fourth  rooter  achieves  a  -31.1  dB  (1.7  dB)  level,  which 
is  about  a  12  dB  improvement  beyond  two  phase  deaber¬ 
ration  algorithms.  This  level  is  very  desirable  for  a  high 
quality  ultrasound  echo  scanner.  BPT-rooter  further  im¬ 
proves  BPT  by  another  5  dB. 

The  restored  mainlobe  diffraction  levels  among  different 
groups  are  similar  when  different  correction  methods  are 
used.  But  the  improvements  are,  in  general,  6  dB  larger 
in  the  scatter  group  than  that  of  the  other  two  groups 
because  the  mainlobe  shape  is  distorted  at  zero  dB  in  the 
scatter  only  group. 

3.  Depth- dependent  Distortion  and  Correction  Scat¬ 
tered  energy  increases  with  the  propagation  depth  D.  Us¬ 
ing  first  order  approximation,  it  is  shown  analytically  in  [3] 
that  wavefront  amplitude  variance  ,  and  phase  variance 
or  arrival  time  variance  al,  linearly  increase  with  depth. 

is  calculated  as  follows;  calculate  energy  at  each  (x,  y) 
position  at  the  aperture  by  summing  the  energy  over  the 
time  window,  calculate  amplitude  at  each  (x,  y)  position  by 
taking  the  square  root  of  the  energy,  calculate  the  reference 
amplitude  profile  which  is  the  fourth  order  polynomial  fit 
of  the  measured  amplitude  profile,  and  calculate  normal¬ 
ized  amplitude  variance  of  the  measured  amplitude  profile 
divided  by  the  reference  amplitude  profile.  The  calculation 
of  cr^  is  described  in  Section  III,  C.  Fig.  8  shows  the  linear 
regression  plots  of  log(0'J)  and  log((7];^^)  vs.  log(D).  A  to¬ 
tal  of  16  independent  measurements  with  sample  thickness 
ranging  from  1  to  4  cm  (see  Section  III,  A)  are  included 


y  =  -  1.40  +  0.73X  R''2  =  0.74 


Fig.  8.  Linear  regression  of  (a)  logfc?)  vs.  log(r>)  and  (b)  logfcr;^^) 
vs.  log(O).  D  is  in  millimeters,  ajrival  time  variance  o'?  is  in  ns-  and 
0-2  is  the  normalized  amplitude  variance  (no  dimension). 
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pio*.  9.  ER  measurement.  Mainlobe  region  of  distorted  image  is  iden¬ 
tified  by  extrapolating  between  Image  peak  and  the  lowest  level  of 
diffraction  pattern  (see  the  two  dashed  lines).  Mainlobe  and  side- 
lobe  energies  of  the  distorted  image  are  obtained  by  integrating  the 
energy  inside  and  outside  the  mainlobe  region. 

in  the  regression  study.  Two  outliers  are  removed  from 
pig  8(a).  The  correlation  coefficients  are  0.55  and  0.74, 
respectively.  The  slopes  are  in  the  order  of  unity.  In  both 
plots,  the  linear  relationships  are  statistically  significant  . 
The  low  correlation  coefficient  in  Fig.  8(a)  is  due  to  the 
difficuty  of  estimating  arrival  time  profiles  when  the  wave¬ 
forms  are  highly  distorted.  However,  the  general  growth  of 

with  depth  D  indicates  that  scattered  energy  increases 

with  D. 

A  quantitative  measure  of  scattered  energy  at  the  image 
plane  is  the  energy  ratio  ER,  which  is  defined  as  the  en- 
ergy  outside  the  main  image  lobe  to  the  energy  inside  the 
main  image  lobe.  It  is  shown  analytically  in  [3]  that  ER 
linearly  grows  with  the  propagation  distance  D,  ER  is  di¬ 
rectly  related  to  image  contrast  and  is  calculated  by  taking 
the  ratio  of  energies  measured  outside  the  mainlobe  region 
and  inside  the  mainlobe  region.  The  mainlobe  region  is 
identified  by  linear  extrapolation  between  the  peak  image 
value  and  the  outer  contour  of  the  mainlobe  diffraction 
pattern  at  the  lowest  recognizable  level.  Fig.  9  illustrates 
the  measurement  procedure.  Fig.  10(a)  upper  curve  is  the 
linear  regression  plot  of  ER  vs.  D  without  correction.  One 
outlier  was  removed  from  the  data.  ER  in  water  image  is 
0.088,  which  is  included  in  the  data  points  to  represent 
system  performance  at  zero  depth.  The  correlation  coeffi¬ 
cient  is  0.49  and  the  linear  relationship  is  significant.  The 
low  correlation  coefficient  is  due  to  the  difficulty  of  iden¬ 
tifying  the  diffraction  pattern  for  some  severely  distorted 
images.  In  Fig.  10,  all  regression  lines  are  extrapolated  to 
.50  mm  propagation  depth,  which  is  considered  the  average 
penetration  depth  of  breast  under  compression. 

The  performances  of  descattering  algorithms  deterio¬ 
rate  with  propagation  depth.  The  rest  of  the  regression 
curves  in  Fig.  10(a)  and  all  regression  curves  in  Fig.  10(b) 
show  ER  vs.  propagation  depth  when  seven  different  cor¬ 
rection  methods  are  used.  The  regression  curve  obtained 
from  PC  show  in  Fig.  10(a)  is  also  included  in  Fig.  10(b) 
for  comparison.  ERs  in  water  alone  data  provided  the  val¬ 
ues  at  zero  depth.  The  correlation  coefficients  of  regression 
curves  are  indicated  in  Fig.  10  by  R.  The  linear  relation¬ 
ships  are  statistically  significant  for  all  curves.  The  average 
values  of  ER  (std.  dev.)  are  given  in  Table  IV. 


iThe  linear  relationship  is  statistically  significant  (a  =  0.05)  if  the 
correlation  coefficient  is  greater  than  0.53  for  sample  size  14  and 
0.497  for  sample  size  16. 


propagation  depth  (mm) 

(a) 


propagation  depth  (mm) 

(b) 

Fi*^.  10.  (a)  Linear  regression  curves  of  ER  vs.  propagation  depth 
obtained  from  without  correction  and  three  correction  procedures, 
(b)  Results  of  ER  vs.  propagation  depth  obtained  from  five  correction 
procedures.  The  vertical  scales  of  (b)  is  the  predicted  image  contrast 
(see  text). 


The  perforinarice  of  wideband  DSA  is  the  lowest  among 
all  correction  procedures.  Data  sets  of  wideband  DSA,  PC 
and  PC-rooter  are  clearly  separated  with  averages  103% 
(std.  33.9%),  25%  (6.9%).  and  6.8%  (1.6%).  respectively. 
The  ER  ratios  of  wideband  DSA  and  PC,  and  PC  and 
PC-rooter  are  4.1  and  3.7,  respectively.  These  results  agree 
very  well  with  the  factor  of  4  prediction  given  by  analysis 
in  [3].  The  average  ER  after  TDC  is  19%  more  than  that  of 
PC  while  the  averages  of  the  lowest  level  at  which  mainlobe 
diffraction  shape  is  maintained  are  the  same  when  TDC 
and  PC  are  used.  This  implies  that  the  residual  phase  error 
due  to  imperfect  phase  correction  contributes  primarily  to 
energy  in  the  sidelobe  lobe  region.  BPT  reduces  the  aver¬ 
age  ER  by  12.2%  below  that  of  PC  and  is  the  best  among 
tL  linear  operations.  The  performance  of  PC-rooter  at 
the  aperture  is  better  than  that  of  BPT.  At  50  mm  depth, 
the  PC-rooter  curve  reaches  0.12  which  is  about  the  aver¬ 
age  ER  level  with  BPT  correction  at  the  average  2.7  cm 
propagation  depth.  Therefore,  this  ER  is  likely  to  result 
in  24  dB  image  contrast  at  this  depth.  At  50  mm,  the  re¬ 
gression  curves  of  BPT-rooter  and  PC- fourth  rooter  reach 
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TABLE  IV 

Statistics  of  ER  Before  .and  After  Corrections. 


vv/o  Corr. 

TDC 

PC 

VViclebancl  DS.A 

BPT 

PC-rooter 

BPT-rooter 

PC- 

fourth  rooter 

•2.14.7% 

(std.  130.1%) 

44.2% 

(13.4%) 

2.5.1% 

(6.9%) 

10.3. .5% 
(.33.9%) 

12.8% 

(2.9%) 

6.8% 

(1.6%) 

3.7% 

(0.7%) 

2. 1% 
(0.4%) 

0.06  and  0.03,  respectively.  These  ERs  are  about  the  aver- 
age  levels  of  PC-rooter  and  BPT-rooter  obtained  from  the 
average  2:T  cm  depth.  These  results  indicate  that  focusing 
with  image  contrast  of  25  to  29  dB  throughout  the  50  mm 
propagation  depth  may  be  achievable  by  taking  the  low 
order  root  on  wavefront  amplitude  in  combination  with 
phase  deaberration  at  the  aperture  or  with  BPT. 

Image  contrast,  the  lowest  level  of  the  mainlobe  diffrac¬ 
tion  shape,  increases  monotonically  with  ER.  The  relation 
of  the  two  can  be  predicted  by  using  ER  and  image  con¬ 
trast  obtained  with  PC  (or  any  correction  method)  as  a 
base  line  (bold  numbers  in  Table  V).  For  example,  the  mea¬ 
sured  image  contrast  with  PC  is  19.5  dB.  The  ER  ratio  of 
PC  alone  to  PC-rooter  is  25.1/6.8  =  3.T,  corresponding  to 
an  increase  of  5.7  dB.  Thus  the  predicted  image  contrast 
is  19.5  dB  -h  5.7  dB  =^,‘25.2  dB,  which  closely  agrees  with 
the  measured  value  —25.4  dB.  Table  V  lists  the  measured 
ER  (first  row),  predicted  image  contrast  from  the  mea¬ 
sured  ER  (second  row),  and  the  measured  image  contrast 
(third  row).  The  predicted  and  measured  image  contrasts 
are  in  reasonably  good  agreement  except  for  the  pair  with¬ 
out  correction.  The  reason  is  due  to  ER  measurement  error 
in  images  without  correction.  For  some  severely  distorted 
images,  it  is  often  difficult  to  identify  the  mainlobe  region. 

Refracted  energy  exhibits  a  different  depth  distortion 
effect  than  scattered  energy.  The  strengths  and  the  num¬ 
bers  of  interference  lobes  depend  on  orientations,  curva¬ 
tures,  sizes,  and  numbers  of  refractive  bodies  in  the  insom- 
fied  medium,  and  may  not  be  directly  related  to  propaga¬ 
tion  depth.  Although  the  chance  of  incident  rays  passing 
through  more  refractive  bodies  increases  with  the  propa¬ 
gation  depth,  the  strengths  of  interference  lobes  resulting 
from  multiple  refraction  may  not  be  significant  compared 
with  those  resulting  from  single  refraction.  One  quantita¬ 
tive  measure  of  coherent  interference  phenomenon  is  the 
ratio  of  a  large  interference  lobe  or  artifact  to  image  jobe. 
For  nine  samples  that  were  determined  to  have  significant 
levels  of  coherent  interference,  the  ratio  varied  randomly 
from  sample  to  sample  and  showed  no  increase  with  thick¬ 
ness.  The  average  strength  (peak  value)  is  10  dB  (std. 
4  dB).  Another  quantitative  measure  is  the  total  number 
of  large  refractive  lobes  in  each  image  (Table  II).  The  av¬ 
erage  is  2  and  standard  deviation  is  0.53.  The  lobe  number 
is  independent  of  propagation  depth.  The  refractive  lobes 
are  tightly  clustered  around  the  main  lobe  with  an  average 
radius  of  5.6  mm  (std.  1.8  mm)  which  is  about  2  degrees. 
This  is  because  the  speed  variations  across  glandular  tissue 
and  fat  boundaries  are  small  (5-10%). 


V.  Discussion 


A.  Comparison  of  TDC  and  PC 


Experimentally  we  found  that  PC  and  TDC,  which 
leave  the  magnitude  of  the  wavefront  unchanged,  have 
similar  performance  with  respect  to  beamwidth  but  PC 
is  consistently  superior  with  respect  to  energy  ratio.  The 
following  analysis  explains  the  reason. 

We  model  the  received  pulse  as  the  transmitted  pulse, 
p{t).  convolved  with  the  medium  induced  distortion,  i.e., 
r{t,xi)  =  p{t)*  m{t,Xi),  where  m{t,Xi)  is  the  impulse  re¬ 
sponse  of  the  medium  from  the  source  position  to  the  re¬ 
ceiver  element  Xi-  Assume  that  system  noise  is  negligible. 
The  estimation  of  the  arrival  time  profile  of  r{t,x)  across 
the  receiving  aperture  is  calculated  by  cross  correlation  of 
neighboring  signals  as: 


0  =  max 

r 


*  m{t,  Xi)  X  p(f)  *  m(t  —  T ,  Xj)dt 


The  quality  of  the  estimated  profile  is  highly  affected 
by  the  waveform  distortion  m(t,  Xi)  across  the  receiv¬ 
ing  aperture,  i.e.,  the  estimation  of  arrival  time  profile  is 
coupled  with  the  wavefront  amplitude  distortion.  In  the 
frequency  domain,  R{uj,Xi)  =  P{ai)  x  M\pj,Xi),  where 
the  medium  distortion  M{uj,Xi)  has  amplitude  a{uj,Xi) 
and  phase  exp{jujto{xi)).  Since  exp(jwto(^))  is  decoupled 
from  the  amplitude  of  R{uj,  Xi)  by  multiplication,  the  phase 
conjugation  operation  at  each  frequency  completely  com¬ 
pensates  to{xi)  independent  of  amplitude  error.  Therefore, 
PC  is  better  than  TDC  in  correcting  residual  phase  error, 
which  primarily  contributes  to  the  energy  in  the  sidelobe 
region. 


B.  Importance  of  Better  Phase  Correction 

We  have  shown  in  Table  IV  that  PC,  the  best  phase  de¬ 
aberration  procedure  at  the  aperture  and  ideal  for  a  point 
target,  improves  ER  19%  more  than  TDC,  while  the  aver¬ 
ages  of  the  lowest  level  at  which  mainlobe  diffraction  shape 
is°maintained  are  the  same  by  using  both  procedures.  This 
implies  that  residual  phase  error  due  to  imperfect  phase 
correction  of  TDC  primarily  contributes  to  the  energy  in 
the  sidelobe  region.  Therefore,  better  phase  deaberration 
procedures  that  improve  the  performance  of  TDC  are  valu- 
able  to  improve  the  contrast  resolution. 

The  procedure  to  calculate  TDC  profile  used  in  this  pa¬ 
per  is  ad  hoc.  Steps  2  and  3  in  the  TDC  calculation  (Sec¬ 
tion  III)  smooth  the  arrival  time  fluctuation  estimate  to  a 
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TABLE  V 

Relation  of  ER  and  Image  Contrast  (Predicted  and  i\[easured). 


w/o  Corr. 

TDC 

PC 

Wideband  DSA 

BPT 

PC-  rooter 

BPT- rooter 

PC- 

fourth  rooter 

ER 

214.7% 

44.2% 

25.1% 

103.5% 

12.8% 

6.8% 

3.7% 

2.1% 

(predicted) 

-10.2  (dB) 

-17.1 

-19.5 

-13.3 

—22.4 

-25.2 

-27.8 

-30.3 

(measured) 

-4.2  (dB) 

-19.3 

-19.5 

-14.8 

-24.0 

-25.4 

-29.4 

-31.1 

certain  extent.  To  test  the  effect  of  the  ad  hoc  smoothing, 
another  TDC  estimate  procedure  was  used.  We  call  this 
procedure  nonsmoothing  TDC.  In  this  procedure,  the  ref¬ 
erence  waveform  selection  in  the  previous  TDC  estimate 
was  unchanged,  but  steps  2  and  3  were  replaced  by  simply 
calculating  the  cross-correlation  peaks  of  correlation  func¬ 
tions  betw'een  geometrically  corrected  waveforms  and  the 
reference  waveform.  The  image  formation  was  the  same  as 
before.  Contrast  and  ER  after  using  nonsmoothing  TDC, 
were  measured.  The  average  level  of  the  lowest  mainlobe 
diffraction  shape  is  —18.9  dB  (std.  2.4  dB)  and  the  average 
ER  is  41.4%  (11.2%).  The  corresponding  results  of  using 
smoothing  TDC  are  —19.3  dB  (std.  2.4  dB)  and  44.2% 
(13%).  Thus,  smoothing  does  not  significantly  improve  im¬ 
age  quality.  The  reasons  are  an  average  86%  of  pulses  are 
highly  correlated  with  the ‘reference  waveform  (correlation 
coefficient  0.7)  and  arrival  time  fluctuations  calculated  at 
these  positions  are  reasonably  good,  and  for  14%  of  pulses 
that  are  not  correlated  well  with  the  reference  waveform, 
estimation  errors  of  arrival  time  fluctuation  occur  at  these 
positions.  But  signals  at  these  positions  are,  in  general, 
weak  and  therefore  arrival  time  estimation  errors  at  these 
positions  contribute  little  to  image  quality. 

C.  Effects  of  Amplitude  Weight  of  Compensation  Vector 


Fig.  11.  Effect  of  refraction  upon  puise  echo  imaging  (simplest  case 
with  no  interference  echo  signal  coming  back  from  targets  insonified 
by  subbeam). 


treated  exactly  as  in  one-way  propagation,  i.e.,  they  are 
decomposed  using  Fourier  decomposition  technique  and 
TIF  can  be  implemented  at  each  discrete  frequency.  TIF 
does  not  require  a  point  target  in  the  medium.  Our  ex¬ 
periments  with  square  root  amplitude  weighting  and  the 
time  domain  correction  algorithm  TDC  have  shown  that 
ER  can  be  reduced  by  15.6%  on  average  relative  to  TDC 
alone. 

E.  Significance  of  Coherent  Interference 
in  Pulse-echo  Imaging 


The  complex  weights  of  MF,  PC,  and  PC-rooter  are,  re¬ 
spectively,  Af,  1  and  each  multiplied  by  exp(-j0). 

The  signal  after  weighting  is  A^,  Af  and  A^p .  By  using 
first  order  approximation,  it  is  shown  in  [3|  that  the  vari¬ 
ance  of  A^j  is  about  four  times  larger  than  cr\^  while  the 

variance  of  “  is  about  four  times  smaller  when  a  narrow- 
band  waveform  is  used,  where  cr\^  denotes  the  variance  of 
Af.  These  increases  and  decreases  in  amplitude  variance 
account  for  the  linear  increase  and  decrease  in  sidelobe 
energy  ER  [3],  Experimentally,  we  measured  a  4.1  and  3.7 
increase  and  decrease  in  ER  compared  with  the  ER  of  us¬ 
ing  PC  w'hen  wideband  DSA  and  PC-rooter  are  used,  re¬ 
spectively  (Table  IV).  The  results  agree  with  the  analysis 
and  indicate  that  intelligent  use  of  amplitude  weighting  is 
crucial  for  improving  image  contrast. 

D.  Application  of  TIF  to  Pulse-echo 
and  Diffuse  Scattering  Medium 

TIF  is  applicable  to  wavefront  compensation  on  re¬ 
ceive  in  pulse-echo.  The  received  wideband  waveforms  are 


A  coherent  interference  phenomenon  was  demonstrated 
in  the  1-way  transmission  model  and  experiments.  Its  im¬ 
pact  upon  2-way  echo  scanning  is  illustrated  in  Fig.  11.  A 
transmitting  beam  illuminates  a  tumor  or  target  T  (part 
a).  Because  of  refraction,  a  subbeam  is  split  from  the  orig¬ 
inal  beam  and  insonifies  targets  along  di.  On  reception,  if 
no  echo  signal  is  coming  back  from  the  targets  illuminated 
by  the  subbeam  (simplest  case),  the  receiving  beam  upon 
T  (part  b)  may  split  in  the  same  fashion  as  the  transmit¬ 
ting  beam  (reciprocity)  and  a  dual  image  may  result.  This 
simplest  case  corresponds  exactly  to  the  situation  in  the  1- 
way  transmission  experiments  reported  in  this  paper,  and 
in  [13]-[15],  where  an  active  point  source  instead  of  a  pas¬ 
sive  point  target  is  used.  In  reality,  refracted  energy  from 
the  targets  insonified  by  the  subbeam  will  arrive  at  angles 
of  1  to  3°  from  the  target  direction  and,  therefore,  in  the 
sidelobe  region  of  array.  The  sidelobe  level  of  the  receiv¬ 
ing  beam  will  be  ~  —30  dB  as  achieved  by  the  procedure 
described  in  this  paper.  The  average  level  of  refracted  mul¬ 
tipath  is  found  to  be  —10  dB  relative  to  the  target  image 
and  therefore  the  average  multipath  signal  level  entering 
the  system  is  ~  (10  +  30)  dB  below  the  echo  strength  of  the 
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illuminated  target.  Cancellation  algorithms  [22]  must  fu^ 
ther  suppress  these  multipath  signals  to  achieve  a  -60  clB 
contrast  resolution,  which  is  essential  for  preventing  con¬ 
tamination  of  otherwise  black  cysts  with  scattered  and/or 
refracted  echo  energy  and  thereby  causing  them  to  look 
like  speckled  tumors. 


VI.  Conclusions 


comments  and  suggestions  on  the  manuscript.  The  authors 
would  also  like  to  thank  the  reviewers  whose  comments  and 
suggestions  have  improved  the  writing  of  the  paper.  The 
in  vitro  measurements  were  done  with  the  collaboration  of 
Ms.  Laura  M.  Hinkelman  and  Professor  Waag.  Invaluable 
consultation  was  provided  by  Dr.  Kai  Thomenius,  Direc¬ 
tor  of  Research,  Interspec,  Inc.,  Ambler,  PA.,  a  division  ot 
ATL,  in  the  course  of  this  study. 


The  evidence  of  both  scattering  and  refraction  is  ob¬ 
served  from  point  source  images  obtained  from  in  vitro 
breast  samples  with  thickness  varying  from  1  to  4  cm. 
Scattering  reduces  target  strength,  broadens  the  mainlobe 
and  increases  background  level  and,  therefore,  lowers  the 
image  contrast.  Refraction  creates  coherent  multipath  in¬ 
terference  that  produces  false  targets  in  addition  to  the 
true  target  in  the  image.  Scattered  energy  increases  with 
the  propagation  depth  while  refracted  energy  does  not  ap¬ 
pear  to  increase  with  the  depth. 

Phase  deaberration  algorithms  are  useful  to  partially 
remove  scattered  energy  and  build  up  the  strength  of  the 
coherent  image.  As  a  result,  the  diffraction-limited  image 
is  restored  up  to  -20  dB  level  for  average  thickness  of 
2.7  cm.  Better  phase  deabecration  procedures  are  valuable 
to  improve  energy  ratio  and  therefore  image  contrast.  Im¬ 
provement  of  a  diffraction-limited  image  to  the  -30  dB  re¬ 
gion  requires  wavefront-deaberration  algorithms  that  can 
further  remove  scattered  energy  caused  by  wavefront  am¬ 
plitude  distortion  and  consequently  strengthen  the  coher¬ 
ent  image.  BPT  reduces  wavefront  amplitude  variance  and 
provides  5  dB  improvement  in  addition  to  that  of  using 
phase  deaberration  procedures  at  the  aperture.  Wideband 
DSA  increases  wavefront  amplitude  variance  and  the  result 
is  5  dB  worse  than  phase-deaberration  procedures  alone 
for  the  tissue  thickness  studied.  TIF,  which  combines  low- 
order-root  spectrum-amplitude  compression  with  phase- 
deaberration  or  BPT,  can  further  reduce  the  amplitude 
variance  and  restore  the  diffraction-limited  image  to  the 
-30  dB  region.  Linear  regression  analysis  of  ER  vs.  prop¬ 
agation  depth  indicates  that  high  quality  focusing  as  far 
as  50  mm  propagation  depth  is  possible. 

The  ratio  of  multipath  artifact  to  image  lobe  is  im¬ 
proved  by  an  average  of  10  dB  when  phase  deaberration  is 
used  and  further  improved  by  another  10  dB  when  the  low 
order  root  upon  wavefront  amplitude,  in  addition  to  phase 
deaberration  or  in  addition  to  backpropagation  and  phase 
deaberration,  is  used.  The  improvement  comes  from  the 
coherent  build-up  of  the  target  strength  by  descattermg. 

TIF  is  applicable  to  wavefront  compensation  on  receive 
in  pulse-echo  and  experiments  are  currently  underway. 
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Abstract  - The  objective  of  this  research  is  to  investigate  the  potential  use  of 

combining  near  infrared  (NIR)  diffuse  light  and  ultrasound  methods  to  increase  the 
detection  sensitivity  and  specificity  of  small  breast  lesions  by  conventional  ultrasound.  A 
novel  imaging  system  which  combines  these  two  methods  has  been  introduced.  A  NIR 
phased  array  positioner  capable  of  providing  the  line-of-sight  to  a  small  target  in  a 
scattering  medium  is  incorporated  into  an  ultrasound  imaging  system.  This  combined 
system  has  been  tested  using  models  consisting  of  a  single  emulated  tumor  of  a  few  mm 
size  buried  in  a  scattering  phantom.  In  each  test,  the  combined  system  was  mechanically 
translated  to  scan  for  the  small  object.  At  each  scanning  position,  optical  amplitude  and 
phase  data  and  an  ultrasound  image  were  acquired.  The  object  line-of-sight  information 
was  coded  into  ultrasound  images  after  scanning.  Ultrasound  images  without  and  with 
line-of-sight  information  were  graded  on  a  five-point  scale  by  two  radiologists.  Receiver- 
Operator  Characteristic  (ROC)  curves  of  ultrasound  and  combined  methods  were 
measured.  ROC  curves  of  the  combined  method  were  much  higher  than  ultrasound  alone. 
Average  areas  under  ROC  curves  of  ultrasound  alone  and  combined  method  were  51.7% 
and  82.5%,  respectively.  This  difference  in  area  is  statistically  significant  (p<0.05).  In 
conclusion,  the  combined  method  affords  a  significant  improvement  in  detection  sensitivity 
and  specificity  compared  with  conventional  ultrasound  in  this  phantom  study. 
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I.  INTRODUCTION 


Ultrasound  is  frequently  used  in  conjunction  with  mammography  to  differentiate 
simple  cysts  from  solid  lesions.  When  all  criteria  for  a  simple  cyst  are  strictly  adhered  to, 
the  accuracy  of  ultrasound  is  96%-100%  [1].  However,  results  with  respect  to  early 
detection  of  small  solid  lesions  remain  disappointing  [2].  Ultrasound  sensitivity  in  the 
detection  of  palpable  solid  masses  less  than  2  cm  in  size  is  71%  which  is  about  10%  less 
than  the  sensitivity  of  physical  examination  [3].  The  sensitivity  in  the  detection  of  small 
invasive  lobular  carcinomas  (<1  cm)  is  only  25%  [4].  Ultrasound  specificity  is  low  due  to 
the  considerable  overlap  of  benign  and  malignant  lesion  characteristics  [5-6].  This  may  be 
caused  by  the  fact  that  even  when  tissues  are  pathologically  different,  their  ultrasound 
properties  may  be  quite  similar  [7]. 

Medical  Optical  Imaging  and  Spectroscopy  is  a  rapidly  advancing  field  [8].  The 
emphasis  is  on  diffusion  waves  in  the  near  infrared  (NIR)  region:  essentially  750  to  850 
nm  are  required  where  a  window  of  the  lowest  tissue  absorption  permits  the  penetration  of 
light  through  the  necessary  10-20  cm  of  a  large  human  breast  [9].  NIR  light  is  expected  to 
have  high  sensitivity  and  specificity  because  certain  characteristics  differ  in  tumors  and 
normal  breast  tissue:  1)  The  amount  of  blood  needed  to  serve  the  tumor's  metabolic  needs 
is  increased  over  that  of  normal  background  tissue;  2)  Hemoglobin  desaturation  in  tumors 
is  increased  due  to  the  high  oxygen  demand  of  cancers;  3)  Light  scattering  properties  of 
the  tumor  are  expected  to  be  enhanced  due  to  the  increased  mitochondria  population  with 
respect  to  the  background  of  normal  cells  [10];  4)  The  angiogenesis  phenomenon  [11] 
allows  tumor  identification  by  observing  the  delivery  of  contrast  agents,  such  as  the  well- 
characterized  NIR  absorber  indocyanine  green  (ICG),  through  the  permeable  vascular  beds 
to  the  extravascular  space  surrounding  the  tumor.  The  intense  absorption  band  of  ICG  at 
800  nm  affords  a  high  contrast  of  the  tumor  extravascular  space  as  compared  with  the 
background  capillary  network  of  the  normal  tissue.  In  addition,  ICG  has  an  interesting 
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redshifted  fluorescence  property  [12].  It  absorbs  radiation  at  the  source  wavelength  of  780 
nm  and  reradiates  photons  at  830  nm.  In  this  process,  the  tumor  extravascular  space  is 
converted  into  a  source  of  secondary  diffuse  waves  with  reduced  background.  Finally,  the 
contrast  agents  could  be  linked  to  tumor-specific  antibodies  or  to  specific  sites  in  the  tumor 
such  as  their  mitochondria  membranes. 

A  simple  optical  phased  array  system  (OPAS)  for  object  positioning  has  been 
constructed  [13]  that  measures  the  interference  of  intensity-modulated  diffuse  light  signals 
at  detector  positions  equidistant  from  two  spatially  separated  sources.  In  the  absence  of 
heterogeneities  (Fig.l.  a)  the  signals  will  cancel,  and  the  null  plane  or  line  is  always  at  the 
bisecting  plane  or  line  and  the  deviation  between  null  and  biosecting  lines  is  zero.  If  a 
hidden  absorber  is  present  in  the  medium  (Fig.  lb),  the  balance  of  the  sources  on  the 
bisecting  line  is  perturbed  and  the  null  line  will  deviate  from  it.  This  deviation,  depending 
on  the  relative  positions  of  the  sources  and  the  absorber,  is  used  as  an  indicator  in  the 
scanner  to  detect  an  inhomogeneous  object  with  respect  to  the  background.  The  line-of- 
sight  of  the  absorber  is  obtained  from  this  newly  formed  null  line  through  scanning.  This 
type  of  "pointing"  device  affords  highly  sensitive  localization  of  a  few  mm  diameter  tube 
obscured  deep  within  a  highly  scattering  material  to  an  accuracy  of  better  than  1°.  In  vivo 
experiments  have  demonstrated  clear-cut  amplitude  nulls  and  sharp  phase  transitions 
[13],[34].  Complete  information  about  object  shapes  and  optical  properties  requires  time- 
consuming  tomographic  approaches  to  data  acquisition  and  inversion  [14-15].  Because  of 
the  diffuse  nature  of  light  propagation  in  human  tissue,  high  resolution  optical  images  are 
intrinsically  difficult  to  obtain.  The  NIR  light  is  strongly  scattered  by  human  tissues. 
Typical  reconstructions  can  distinguish  simple  structures  of  approximately  1  cm  in  size; 
sharp  edges  are  typically  blurred  by  a  few  millimeters  [7]. 
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Ultrasound,  on  the  other  hand,  has  real  time  imaging  capability  and  the  imaging 
capability  of  ultrasound  equipment  has  been  steadily  improving.  More  recently,  the  result 
reported  by  Stavros  et  al.  is  very  encouraging  [16].  Although  no  lesion  size  was 
documented  in  [16],  of  125  malignant  lesions,  123  were  correctly  classified  as  intermediate 
or  malignant  (98.4%  sensitivity).  Coupling  OPS  A  into  an  ultrasound  system  will  increase 
ultrasound  detection  sensitivity  and  specificity  while  maintaining  ultrasound  real  time 
imaging  capability  and  resolution.  Furthermore,  the  combined  system  is  low  cost  and 
efficacious,  which  meets  the  challenges  of  health  care  reform. 


II.  BASIC  PRINCIPLES 


The  principle  of  object  detection  by  OP  AS  is  described  briefly  below.  The  density 
of  photons  at  position  r  of  intensity-modulated  light,  U(r),  traveling  in  an  infinitely-large 
homogeneous  scattering  medium  can  be  described  by  diffusion  theory  and  expressed  in  the 
frequency-domain  as  [17] 

U{r) - (1) 

AnvOr 


where  D 


k  = 


t’-a-ico/v 


3(aa  +  Hj' )  I  D 

D  is  the  diffusion  coefficient,  pa  atid  ps'  are  absorption  and  reduced  scattering  coefficients, 
V  is  the  speed  of  light  in  the  medium,  and  m  is  the  modulation  frequency  of  the  light 


source.  Since  k  is  a  complex  number,  Eq.  1  can  be  rewritten  as 


U{r)  =  ■ 


I  -k^r  -ik:r 
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where  kir  corresponds  to  a  phase  shift  produced  by  traveling  in  the  scattering  medium  with 
respect  to  the  source. 


A  phased  array  system  consists  of  one  in-phase  (i.e.  initial  phase=0r>)  and  one  out- 
phase  source  (i.e.  initial  phase=180O)  as  shown  in  Fig.  2.  The  total  photon  density  at 
position  r  is 
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where  r,  ri,  and  rj  are  related  by  rj  =  +a^  -  2racos(e)  and  r2  =  '^r'^  +a^  +  2racos(0) .  In 

this  case,  the  intensities  for  the  two  in-phase  and  out-phase  sources  are  assumed  the  same, 
and  the  separation  between  the  two  sources  is  2a.  Eq.  3  indicates  that  the  photon  density  at 
the  bisecting  line  of  the  array,  i.e.  ri=r2,  is  equal  to  zero  due  to  the  cancellation  of  the  two 
sources;  the  phase  shift  relative  to  the  extremes  of  0°  an  ISO^  is  expected  to  be  90°  and  the 
amplitude  is  theoretically  zero  but  practically  determined  by  the  system  noise.  When  a 
small  heterogeneity  is  present  in  the  medium,  the  balance  of  the  sources  on  the  bisecting 
line  will  be  destroyed  and  the  null  line  will  bias  toward  the  object.  The  object  line-of-sight 
is  obtained  from  the  newly  formed  null  line  through  scanning. 

III.  METHODS 

A.  Experimental  set-up 

Mechanical  assembly  of  the  optically  guided  ultrasound  system  is  shown  in  Fig.3. 
A  coaxial  probe  is  assembled  by  mounting  two  laser  diode  sources  (in-  and  out-phase 
sources)  and  an  ultrasound  transducer  on  a  plastic  holder.  A  Photo-Multiplier  Tube  (PMT) 
optical  detector  is  mounted  on  the  opposite  side  of  the  coaxial  probe.  The  separation  of  the 
source  and  the  detector  can  be  changed  from  several  cms  to  15  cm  depending  on  the  size  of 
a  testing  phantom.  The  coaxial  probe  and  optical  detector  assembly  is  translated  laterally  by 
a  stepping  motor  in  the  X  direction  to  scan  an  object.  The  coupling  between  the  coaxial 
probe  and  a  scattering  phantom  described  later  was  made  using  the  standard  ultrasound 
coupling  gel.  In  a  typical  experiment,  the  step  size  used  is  1  mm.  The  optical  wavelength 
is  780  nm  and  the  modulation  frequency  is  5  KHz.  The  ultrasound  probe  is  a  64-element 
phased  array  operating  at  5  MHz  center  frequency  and  the  ultrasound  system  used  is  a  GE 
RT  3600  commercial  scanner  located  at  the  Angiography  Section  of  the  Radiology 


5 


Department.  The  focal  depth  of  the  scanner  can  be  changed  from  1  cm  to  4  cm  and  we 
used  4  cm  in  all  experiments. 

B.  Data  acquisition 

The  detailed  description  of  OPS  A  can  be  found  in  Ref.  [18]  and  data  acquisition  is 
described  here.  The  optical  amplitude  and  phase  data  at  each  lateral  position  were  acquired 
through  a  12  bit  A/D  converter  (AD  2210,  Real  Time  Devices,  Inc)  and  stored  on  a  PC. 
The  ultrasound  B-scan  image  at  each  lateral  position  was  captured  using  a  video  board 
(ComputerEyes/RT,  Digital  Vision,  Inc).  The  resolution  of  the  image  grabber  is  200  x  320 
with  8  bit  256  gray  levels. 

C.  Scattering  phantom 

A  scattering  phantom  was  made  of  an  acrylamide  gel  (Boehringer  Mannheim 
Biochemicals  Corporation)  which  is  typically  used  as  an  electrophophresis  gel  for  DNA 
studies  and  is  also  used  to  construct  ultrasound  phantoms  [19-20] .  Scattering  particles  at 
optical  and  ultrasound  wavelength  scales  were  added  to  the  acrylamide  solution  before  the 
gel  was  formed.  First,  the  premixed  acrylamide  and  bis-acrylamide  powder  was  dissolved 
in  distilled  and  degassed  water.  20%  concentration  of  Intralipid  (Kabi  Pharmatica)  was 
added  to  the  acrylamide  solution  to  dilute  the  solution  to  a  0.5%  intralipid  concentration. 
The  optical  scattering  parameter  |is  at  this  concentration  closely  matches  fxs  of  human 
tissue.  Second,  scattering  particles  made  of  glass  spheres  (Jenson  Scientific)  were  added 
to  the  solution.  The  mean  diameters  of  the  spheres  were  112,  150  and  205  |im  and  the 
estimated  particle  density  based  on  the  calculation  given  in  [21]  should  produce  Rayleigh 
amplitude  distribution.  In  order  to  produce  polymerization,  components  of  ammonium 
persulfate  and  tetramethylethylenediamine  (known  as  TEMED)  were  added.  The  outer 
dimensions  of  the  scattering  phantom  were  7.5  cm  x  12  cm  x  14  cm. 
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Objects  of  several  mm  sizes  were  inserted  in  the  scattering  phantom  during 
experiments. 

D.  Experimental  procedure 

Since  we  are  to  compare  the  two  detection  methods,  ultrasound  alone  and 
ultrasound-light,  it  is  necessary  to  select  a  metric  for  comparison.  Receiver-Operator 
Characteristic  (ROC)  curves  are  commonly  used  to  depict  figures  of  merit  for  clinical 
methods  of  detecting  lesions  or  differentiating  among  diseases  [22-23].  Each  ROC  curve 
indicates  the  trade-offs  between  sensitivity  and  specificity  that  are  available  from  a 
diagnostic  method,  and  thus  describes  the  inherent  discrimination  capacity  of  that  method. 
A  higher  ROC  curve  would  indicate  greater  discrimination  capacity,  because  a  larger  value 
of  sensitivity  at  each  value  of  false-positive-fraction  (EPF  =  1 -specificity)  or,  equivalently, 
a  smaller  value  of  EPF  at  each  sensitivity,  can  be  achieved  on  a  higher  curve  if  an 
appropriate  confidence  threshold  is  used.  Similarly,  a  lower  ROC  curve  would  indicate 
less  discrimination  capacity. 

Three  different  experimental  approaches  can  be  used  to  measure  ROC  curves, 
namely,  yes-no,  rating  and  forced  choice  methods  [24].  The  rating  method  is  more 
efficient  and  is  almost  always  used  in  practice  [23].  In  this  approach  the  observer  is 
required  to  select  one  of  several  categories  of  confidence,  usually  represented  by  numerical 
ratings,  to  report  his  impression  of  the  likelihood  that  lesions  or  targets  are  present  or 
absent.  These  categories  can  be  given  qualitative  labels  such  as  (1)  definitely  or  almost 
definitely  absent,  (2)  probably  absent,  (3)  may  be  present,  (4)  probably  present,  and  (5) 
definitely  or  almost  definitely  present.  It  is  not  necessary  for  different  observers  to 
interpret  these  category  labels  in  the  same  way;  the  essential  requirement  is  that  the  labels 
provide  an  unambiguously  ordered  set  of  categories  for  eaeh  observer  s  relative  confidence 
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in  the  two  states  of  truth.  We  have  adapted  the  rating  method  to  measure  ROC  curves  of 
ultrasound  and  ultrasound-light  methods. 

It  is  well  known  that  weak  targets  in  a  sequence  of  images  are  more  readily  detected 
than  in  a  single  image.  Thus  triple  images  taken  in  1  mm  steps  at  the  neighborhood  of  the 
optical  phase  transition  and  amplitude  null  were  selected  in  each  test.  The  criteria  of  triple 
image  selection  are;  (1)  identify  the  minimum  point  of  the  amplitude  null  from  the  optical 
amplitude  curve,  or  in  the  case  that  the  minimum  is  not  unique,  choose  the  point  which  is 
closest  to  90®  from  the  optical  phase  transition  curve;  (2)  select  the  corresponding 
ultrasound  image  at  the  same  lateral  position  as  the  second  image  in  the  triple;  (3)  select 
images  taken  1  mm  left  and  right  of  the  second  image  as  the  first  and  third  images  in  the 
triple.  Since  the  optical  assembly  and  ultrasound  probe  are  coaxially  mounted,  the  object  in 
the  triple  images  always  appears  along  the  ultrasound  scan  axis.  To  avoid  giving  any  prior 
knowledge  to  human  observers  about  the  object  location,  we  chose  different  image 
windows  in  different  tests  to  randomize  the  object  location.  In  each  test,  the  image 
windows  used  for  triple  images  were  the  same  and  the  object  position  shift  was  caused  by 
the  1  mm  lateral  translation  of  the  imaging  system. 

All  triple  images  without  and  with  optical  information  were  evaluated  by  2 
experienced  radiologists  in  the  Breast  Imaging  Section  of  the  Radiology  Department  using 
the  evaluation  form  shown  in  Table  I.  They  were  informed  that  the  images  were  obtained 
at  slightly  different  lateral  positions  (1mm)  and  that  they  could  mentally  average  these 
images  to  filter  out  artifacts  and  noise.  For  ultrasound-alone  images,  the  observers  were 
told  that  they  should  look  for  possible  lesions  or  targets  in  the  entire  image.  For  the  images 
coded  with  optical  line-of-sight  information,  the  observers  were  instructed  to  look  for 
possible  lesions  or  targets  in  the  neighborhood  along  the  line-of-sight  direction.  Observers 
were  also  told  that  targets  might  be  hypoechoic  or  hyperechoic.  Triple  images  with  no 
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lesion  or  target  present  were  also  presented  to  observers  during  the  trials.  The  positive  and 
negative  sets  of  images  were  presented  to  observers  in  random  order. 

IV.  RESULTS 

Fig.4  shows  optical  amplitude  and  phase  responses  when  concentric  black-plastic 
tubes  of  1.9  and  1.7  mm  diameters  filled  with  0.5%  intralipid  were  inserted  into  the 
scattering  phantom  at  approximately  6  cm  cm  depth.  The  testing  object  is  an  optical 
absorber  which  also  generates  an  inhomogeneous  ultrasound  echotexture,  as  some  tumors 
may  do,  as  shown  in  Fig. 5.  The  sharp  phase  transition  and  amplitude  null  correspond  to 
the  object  line-of-sight.  Three  ultrasound  images  taken  at  phase  transition  and  amplitude 
null  positions  (16,  17,  18  mm)  are  shown  in  Fig. 5.  The  image  of  the  object,  indicated  by 
an  arrow,  is  highly  obscured  by  speckle  noise.  However,  one  observer  picked  up  the 
target  and  marked  category  4  on  the  evaluation  form.  When  the  same  images  coded  with 
optical  line-of-sight  information,  shown  by  dashed  lines  in  Fig.6,  were  presented  to  the 
observer,  his  detection  confidence  increased  from  category  4  to  5.  Another  observer 
completely  missed  the  target  in  ultrasound-alone  images  and  marked  category  1  on  the 
form.  When  the  images  with  object  line-of  sight  information  were  presented,  the  observer 
picked  up  the  target  and  selected  category  3  on  the  form.  This  example  shows  the 
improvement  in  detection  sensitivity  when  the  combined  method  is  used. 

Figs.7-9  show  another  example.  A  4.5  mm  white  Teflon  bar  was  inserted  into  the 
scattering  phantom  and  was  imaged  using  the  combined  system.  The  testing  object  is  an 
optical  scatterer  which  generates  a  slightly  hyperechoic  echotexture  as  shown  in  Fig.8.  The 
optical  amplitude  and  phase  responses  are  shown  in  Fig.7.  The  triple  images  taken  at  the 
neighborhood  of  optical  phase  transition  and  amplitude  null  (20,  21,  22  mm)  are  shown  in 
Fig.8  with  the  object  indicated  by  an  arrow  (upper  one)  in  the  images.  One  observer 
missed  the  target  entirely  and  marked  category  1  on  the  evaluation  form.  When  the  images 
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with  object  line-of-sight  (Fig.9)  were  presented,  the  observer  identified  the  target  and 
circled  category  3  on  the  form.  Another  observer  missed  the  target  completely  and  selected 
an  artifact,  pointed  by  an  arrow  (lower  one),  and  circled  category  4  on  the  form.  Once  the 
images  with  object  line-of-sight  were  presented,  the  observer  identified  the  target  and 
circled  category  3  on  the  form.  This  example  shows  the  combined  method  can  significantly 
improve  detection  sensitivity  and  specificity  of  ultrasound. 

We  have  studied  a  total  of  30  positive  trials  and  21  negative  trials.  In  the  positive 
trials,  the  testing  objects  used  are  black  Delrin  rods,  black  rubber  sticks  of  different 
stiffness,  black  plastic  tubes  or  concentric  tubes  filled  with  2.5|ig/nil  ICG  or  5%  intralipid, 
and  white  Teflon  rods.  The  sizes  of  the  testing  objects  range  from  1.6  mm  to  8.2  mm,  and 
the  average  is  3.64  mm  with  1.97  mm  standard  deviation.  In  the  case  of  non  cylindrical 
object,  geometric  mean  is  used  in  the  calculation.  These  testing  objects  are  either  optical 
absorbers  or  scatterers  which  emulate  to  a  large  extent  tumor  optical  properties.  The 
acoustics  reflectivity  of  these  materials  varies  from  hyperechoic,  heterogeneous  to  slightly 
hypoechoic.  In  the  negative  trials,  images  of  the  scattering  phantom  obtained  at  three 
consecutive  lateral  positions  in  1  mm  steps  were  presented.  ROC  data  of  positive  and 
negative  trials  of  two  independent  readings  obtained  from  the  two  observers  were  shown  in 
Table  II  and  III.  The  first  section,  consisting  of  the  first  row  of  the  table,  lists  the  5 
categories  used  in  the  evaluation.  The  second  section  shows  four  quantities  based  on 
positive  trials;  and  the  third  shows  the  corresponding  quantities  based  on  negative  trials. 
The  first  row  in  section  (2)  shows  the  frequency  of  the  occurrence  of  the  correct  target 
identification  and  the  use  of  the  corresponding  category  when  the  ultrasound  method  is 
used.  The  second  row  is  the  sensitivity  or  conditional  probability  (P(T/t))  of  the  observer 
responding  target  T  at  the  correct  location  when  target  t  at  that  location  was  present.  The 
third  and  fourth  rows  are  the  same  quantities  as  the  first  two  but  the  method  used  is 
ultrasound-light.  The  first  row  in  section  (3)  shows  the  frequency  of  the  occurrence  of 
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false  positives  and  the  use  of  the  corresponding  category  when  the  ultrasound  method  is 
used.  The  second  row  is  1 -specificity  or  condition  probability  (P(T/n))  of  the  observer 
responding  target  T  when  no  target  was  present.  The  third  and  fourth  rows  are  the  same 
quantities  when  the  ultrasound-light  method  is  used. 

The  calculation  of  P(T/t)  using  both  methods  is  strictly  forward  when  the  observer 
identified  the  correct  target  or  did  not  detect  any  target  during  positive  trials.  But  there  were 
cases  in  positive  trials  in  which  the  observer  picked  up  artifacts  instead  of  the  correct  target. 
The  locations  of  artifacts  could  be  anywhere  when  ultrasound  is  used,  but  they  are 
constrained  along  the  object  line-of  sight  when  ultrasound-light  is  used.  For  these  cases, 
we  classified  them  into  category  1  in  positive  trials,  because  with  regard  to  the  correct  target 
the  observer  reported  no  target  present.  Since  these  cases  correspond  to  false  positives,  we 
added  them  to  the  21  negative  trials.  Thus  the  total  negative  trials  of  observer  1  in  first 
reading  is  29  and  second  reading  is  32,  and  observer  2  in  first  reading  is  33  and  second 
reading  is  34.  For  the  total  21  negative  trials,  optical  data  did  not  shown  any  phase 
transitions  and  amplitude  nulls.  Therefore,  these  trials  were  classified  as  category  1  under 
the  ultrasound-light  method. 

Average  empirical  ROC  curves  of  ultrasound  (solid)  and  ultrasound-light  (dashed) 
of  the  two  independent  readings  for  the  two  observers  are  shown  in  Fig.  10.  The  superior 
performance  of  ultrasound-light  is  shown  by  the  higher  ROC  curves.  The  fact  that 
ultrasound  ROC  curves  are  close  to  the  major  diagonal  indicates  that  the  detection 
information  is  very  poor  in  this  phantom  study  and  that  targets  can  only  be  recognized  at  a 
slightly  better  than  chance  level.  The  area  under  the  ROC  curve  provides  a  measure  of 
discriminatory  efficacy:  an  area  of  100%  denotes  an  ideal  system  and  50%  denotes  an 
ineffective  system.  The  calculated  areas  under  ROC  curves  of  ultrasound  and  ultrasound- 
light  of  observer  #1  are  53.9%  and  79.2%,  and  of  observer  #2  are  49.5%  and  85.8%.  The 


mean  values  are  51.7%  for  ultrasound  and  82.5%  for  ultrasound-light,  respectively. 
Following  the  statistical  test  procedure  given  in  [33],  the  difference  in  area  between  the  two 
methods  is  statistically  significant  (p<0.05,  see  Appendix  for  detail). 

In  the  human  observer  study,  each  observer  performed  two  independent  readings 
with  positive  trials  presented  in  slightly  different  orders.  In  the  first  reading,  each  positive 
trial  was  evaluated  in  the  order  of  ultrasound  and  ultrasound-light  instead  of  random.  This 
procedure  took  less  time  for  the  observers  because  they  were  already  familiar  with 
ultrasound  images  before  the  same  images  with  optical  line-of-sight  were  presented.  In  the 
cases  that  the  observers  detected  targets  along  the  line-of-sight  with  ultrasound,  this 
procedure  might  have  some  psychological  effect  on  the  category  they  chose  with 
ultrasound-light  because  they  might  still  remember  the  category  they  selected  with 
ultrasound  method.  To  test  this  effect,  we  presented  ultrasound  and  ultrasound-light 
images  as  well  as  negative  images  to  the  two  observers  in  random  order  at  the  second 
independent  reading.  The  ROC  data  obtained  from  the  two  readings  for  the  two  observers 
were  given  in  Tables  n  and  III  and  the  corresponding  ROC  curves  for  the  two  observers 
are  shown  in  Fig.l  1.  Procedure  #1  (first  reading)  is  the  ordered  test  with  ultrasound  and 
ultrasound-light  images  presented  in  sequence  while  procedure  #2  (second  reading)  is  the 
random  test  with  ultrasound  and  ultrasound-light  images  presented  in  random  order.  The 
results  with  both  testing  procedures  are  consistent  which  indicate  that  the  observers  are 
based  on  their  own  evaluation  criteria  and  are  not  influenced  by  the  ways  that  the  images 
are  presented.  The  areas  under  the  ROC  curves  of  procedures  #1  and  #2  for  observer  1  are 
57.3%  and  50.9%  for  ultrasound  alone  and  78.2%  and  80.2%  for  ultrasound-light,  while 
the  areas  for  observer  2  are  44.6%  and  50.0%  for  ultrasound  alone  and  85.8%  and  85.8% 
for  ultrasound-light.  The  fact  that  the  difference  between  ultrasound  ROC  curves  measured 
with  the  two  procedures  is  larger  than  the  difference  between  ultrasound-light  ROC  curves 


measured  with  the  two  procedures  for  both  observers  reflects  the  chance  performance 
characteristic  of  ultrasound  in  this  study. 

V.  DISCUSSION 

The  model  studies  described  in  the  paper  have  been  designed  to  assess  the  influence 
of  adding  object  line-of-sight  information  upon  detection  by  human  observers.  The  small 
objects  used  in  the  studies  are  optical  anomalies  with  acoustic  reflectivity  varying  from 
hyperechoic,  heterogeneous  to  slightly  hypoechoic.  These  anomalies  emulate  to  a  large 
extent  tumor  optical  properties  and  to  some  extent  tumor  acoustic  properties 
(inhomogeneous  echogenicity  and  slightly  hypoechoic  cases).  However,  these  anomalies 
as  well  as  the  scattering  phantom  in  which  they  are  embedded  are  obviously  inadequate  to 
reproduce  all  the  properties  of  breast  tissue.  Clinical  studies  are  planned  to  further  evaluate 
the  combined  imaging  system. 

The  ultrasound  scanner  used  for  this  study  is  not  a  state-of-the-art  breast  scanner, 
but  it  is  the  one  available  full-time  to  us  for  these  experiments  in  the  research  laboratory. 
However,  the  relative  improvements  in  detection  sensitivity  and  specificity  when  the  optical 
information  is  added  are  independent  from  the  scanner  quality  used  because  the  ROC 
comparisons  of  ultrasound  alone  and  ultrasound-light  are  based  on  a  pair  of  triple  images 
obtained  from  the  same  imagery.  With  a  better  scanner,  we  should  be  able  to  get  similar 
improvements  as  reported  in  this  paper  with  reduced  average  object  size. 

In  the  ultrasound-light  image  display,  the  optical  null  line  was  simply  coded  into  the 
corresponding  ultrasound  images  by  a  dashed  line.  This  display  method  is  simple  but  it 
does  not  use  optical  information  related  to  the  mass  of  the  object.  For  example,  a  large 
mass  with  higher  optical  absorption  and  scattering  will  produce  a  sharper  phase  transition 
and  deeper  amplitude  null  compared  with  a  small  mass.  In  our  future  studies,  we  need  to 
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code  optical  phase  gradient  information  with  a  gray  scale  har  into  ultrasound  images.  A 
higher  gray  level  (large  gradient)  may  suggest  that  observers  look  for  a  larger  object  and  a 
lower  gray  level  (small  gradient)  may  suggest  that  observers  search  for  a  smaller  object  in  a 
larger  area.  However,  the  approximate  relation  between  tissue  mass  sizes  and  optical 
properties  needs  to  be  investigated. 

To  estimate  the  sensitivity  of  OPAS,  we  measured  depths  of  amplitude  nulls  and 
spatial  intervals  of  phase  transitions  for  all  positive  trials.  The  depth  of  the  amplitude  is 
defined  as  the  difference  in  dB  between  the  peak  and  the  null  in  an  amplitude  profile.  Since 
the  two  peaks  on  each  side  of  a  null  were  not  exactly  the  same  (See  Fig.4a),  we  used 
smaller  peaks  in  all  cases.  The  spatial  interval  of  the  phase  transition  was  the  interval 
where  a  phase  curve  rose  from  one  stable  value  and  settled  down  at  another  stable  value 
(See  Fig.4b).  The  two  stable  values  were  obtained  by  averaging  5  to  10  points  on  each 
side  away  from  the  phase  transition.  If  we  denote  the  difference  between  the  two  stable 
values  as  D,  the  spatial  interval  of  the  phase  transition  is  quantitatively  defined  as  the 
difference  between  the  two  spatial  points  where  a  phase  curve  exceeds  10%D  and  reaches 
90%D.  The  measured  average  depth  of  amplitude  is  8.93  dB  (std.  4.58  dB)  and  the 
average  spatial  interval  of  phase  transition  is  4.78  mm  (std.  1.85  mm)  which  is  about  1  mm 
larger  than  the  average  size  of  testing  objects. 

The  system  noise  can  cause  pointing  error  by  shifting  the  null  line  of  OPAS. 
However,  as  long  as  the  null  line  is  confined  in  the  average  transition  region,  the 
ultrasound  detection  would  not  be  affected.  For  example,  if  the  null  line  were  shifted  to  the 
edge  of  the  phase-transition  interval  (2.39  mm),  the  pixel  shift  from  the  center  of  the  target 
in  the  image  were  only  6  points  at  the  typical  60  mm  depth  studied.  As  a  matter  of  fact,  the 
observers  had  the  tendency  to  search  and  circle  a  larger  area  alone  the  line-of-sight  than  the 
true  object  size  on  the  image 
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The  resolving  power  of  2  objects  by  OPAS  is  very  low  and  the  object  separation 
must  be  as  large  as  1  cm  [18].  When  the  object  separation  is  less  than  approximately  1  cm, 
OPAS  will  point  to  the  center  of  mass  of  the  objects.  For  this  reason,  the  combined  system 
is  optimal  for  the  localization  and  imaging  of  an  isolated  small  tumor  in  the  breast.  Several 
clinical  studies  have  been  devoted  to  the  determination  of  multicentric  occurrence  of  breast 
cancers  [25].  McDivitt  [26]  critically  reviewed  these  data  and  concluded  that  the  incidence 
of  invasive  foci  in  areas  of  multicentric  cancers  was  5-10%.  For  these  cases,  the  optical 
amplitude  and  phase  responses  may  show  multiple  peaks  and  the  locations  of  these  peaks 
may  also  be  used  as  references  to  determine  the  centers  of  the  anomalies. 

The  combined  system  has  been  implemented  by  mechanically  translating  an 
ultrasound  phased  array  with  OPAS.  This  choice  accommodates  the  current  set-up  of 
OPAS.  The  optical  null  line  can  be  electronically  steered  by  varying  the  amplitude 
modulation  of  the  two  sources  [18].  However,  scanning  parameters,  such  as  sensitivity  to 
off-line  object,  field  of  view  vs.  source  separation,  etc.,  need  to  be  studied  before  the 
OPAS  scan  mode  can  be  used.  This  is  a  topic  for  future  study. 

With  current  one-dimensional  implementation  of  the  combined  system,  false 
positives  happen  along  the  object  line-of-sight.  With  another  OPAS  deployed  from  a 
different  direction  than  the  combined  probe,  the  object  position  can  be  located  from  the 
intersection  of  two  lines  of  sight.  Therefore,  the  false  positive  rate  with  the  combined 
method  can  be  reduced  significantly. 

Transmission  mode  of  OPAS  is  used  in  our  current  system  set-up.  Reflection 
mode  is  also  possible  where  a  detector  is  placed  at  the  middle  of  two  laser  diode  sources 
[18].  This  deployment  requires  overlapping  geometry  of  an  ultrasound  array  and  the  light 
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detector,  and  can  be  done  by  coupling  the  detector  through  the  ultrasound  array  with  a 
small  light  guide.  Note  that  light  guides  can  be  quite  small,  i.e.,  on  the  order  of  100  micro 
optical  fibers.  The  penetration  depth  of  reflection  mode  is  about  half  that  of  transmission 
mode,  and  it  is  not  a  problem  for  most  breasts  under  soft  compression  when  a  patient  is  in 
a  supine  position. 

There  are  several  possible  applications  of  the  combined  method.  First,  the 
combined  method  may  allow  real  time  implementation  of  current  ultrasonic  wavefront 
compensation  algorithms.  The  female  breast  and  abdominal  wall  are  acoustically 
inhomogeneous  and  that  limits  the  ultrasound  resolution.  Ultrasonic  wavefront 
compensation  algorithms  adaptively  adjust  time  delays  [27-29]  and  amplitude  variations 
[30-31]  of  RF  signals  received  from  different  array  elements  to  compensate  tissue  induced 
distortion  and  therefore  to  improve  spatial  and  contrast  resolution.  But,  these  algorithms 
require  off-line  processing  because  of  extensive  computations  involved.  With  the 
combined  method,  the  compensation  can  be  implemented  in  real  time  when  the  need  for 
correction  is  confined  to  the  neighborhood  of  object  line-of-sight.  Second,  OP  AS 
combined  with  low-frequency  ultrasound  could  be  used  to  detect  and  monitor  a  small 
aneurysm  or  hemorrhage  in  the  brain.  Although  the  optical  method  alone  is  able  to  detect 
and  monitor  an  anomaly  in  the  brain  [32],  the  size  and  shape  information  provided  a 
combined  ultrasound-light  system  would  provide  valuable  information  to  guide  surgery  and 
treatment. 


VI.  SUMMARY 

A  novel  imaging  method  which  combines  ultrasound  and  near  infrared  diffuse  light 
has  been  introduced  in  this  paper.  The  combined  method  utilizes  object  line-of-sight 
information  provided  by  a  highly  sensitive  optical  positioner  to  confine  human  observers 
along  the  object  direction  to  search  for  possible  lesions  or  targets  in  ultrasound  images.  As 


a  result,  the  detection  sensitivity  and  specificity  of  a  small  target  of  a  few  mm  size  have 
been  improved  significantly  over  the  conventional  ultrasound.  The  average  areas  under 
ROC  curves  from  2  radiologists  are  much  higher  for  the  combined  method  (82.5%)  than 
for  ultrasound  alone  (51.7%).  Thus  the  combined  method  affords  a  significant 
improvement  in  detection  sensitivity  and  specificity  compared  with  conventional 
ultrasound. 


1  7 


APPENDIX 


We  have  adapted  the  approaeh  given  in  [33]  to  determine  the  statistical  significance 

of  the  difference  in  average  areas  measured  with  ultrasound  and  ultrasound-light  methods. 
In  the  approach,  the  statistic  power  is  calculated  via  given  by 

rnS-Z,W.  (4, 

where  =  1.645,  for  a  5%  one-sided  test  of  significance,  <5  =  02  ~  .  is  the  difference 

between  areas  measured  with  the  two  methods,  n  is  the  sample  size  used, 

y,  =  - 20f  and  V,  =  — ^  - W].  Z„=  0.84,  1 .28,  1 .645  for 

'  2-01  (l-t-0,)  '  '2-02  (1-1-02)  '  ^ 

80%,  90%  or  95%  statistic  power. 

Taking  the  averages  0i  =  51.7%  and  02  =82.5%  measured  with  ultrasound  and 
ultrasound-light  methods,  5=0.308,  y,  =  0.166  and  V2  =  0.087.  With  n=  30,  Z^  is 

1.469.  Therefore,  greater  than  90%  statistic  power  has  been  achieved,  i.e.,  we  have  a  high 
assurance  (>90%  probability)  that  the  test  of  significance  on  samples  would  yield  a 
statistically  significant  difference  (p<0.05). 
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TABLES 


Table  I.  Evaluation  form 


Please  check  for  possible  lesions  or  targets. 

Target#! 

1.  No  lesion  or  target  present 

2.  Probably  no  lesion  or  target  present 

3.  Lesion  or  target  may  be  present 

4.  Lesion  or  target  probably  present 

5.  Lesion  or  target  definitely  present 

Target  #2 

1 .  No  lesion  or  target  present 

2.  Probably  no  lesion  or  target  present 

3.  Lesion  or  target  may  be  present 

4.  Lesion  or  target  probably  present 

5.  Lesion  or  target  definitely  present 
Target  #3... 

Target  #4... 

Lesions  or  targets  may  be  hyperechogenic  or  hypoechogenic.  If  you 
identify  any  lesions  please  circle  and  number  them  on  images. 
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Table  II  ROC  data  of  observer  #  I  (First  Reading) 
A  total  of  30  positive  trails  and  29  negative  trails 


Section(l) 

Category 

1 

2 

3 

4 

5 

Section(2) 

US: 

Occurrence 

16 

3 

5 

3 

3 

P(T/t) 

1 

0.467 

0.367 

0.200 

0.100 

US-Light:  Occurrence 

8 

2 

6 

9 

5 

P(T/t) 

1 

0.733 

0.667 

0.467 

0.167 

Section  (3) 

US: 

Occurrence 

16 

9 

4 

0 

0 

P(T/n) 

1 

0.448 

0.138 

0 

0 

US-Light: 

Occurrence 

22 

1 

4 

2 

0 

P(T/n) 

1 

0.241 

0.207 

0.069 

0 

ROC  data  of  observer  #1  (Second  Reading) 

A  total  of  30  positive  trails  and  32  negative  trails 


Section!  1) 

Category 

1 

2 

3 

4 

5 

Section(2) 

US: 

Occurrence 

17 

2 

3 

4 

4 

P(T/t) 

1 

0.433 

0.366 

0.266 

0.133 

US-Light:  Occurrence 

7 

4 

7 

6 

6 

P(T/t) 

1 

0.767 

0.633 

0.400 

0.200 

Section  (3) 

US: 

Occurrence 

14 

7 

11 

0 

0 

P(T/n) 

1 

0.563 

0.344 

0 

0 

US-Light: 

Occurrence 

24 

3 

4 

1 

0 

P(T/n) 

1 

0.250 

0.156 

0.031 

0 
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Table  IE  ROC  data  of  observer  #2  (First  Reading) 
A  total  of  30  positive  trails  and  33  negative  trails 


Section(l) 

Category 

1 

2 

3 

4 

5 

Section(2) 

US: 

Occurrence 

15 

0 

4 

4 

7 

P(T/t) 

1 

0.5 

0.5 

0.367 

0.233 

US-Light:  Occurrence 

6 

2 

3 

7 

12 

P(T/t) 

1 

0.800 

0.733 

0.630 

0.400 

Section  (3) 

US: 

Occurrence 

4 

3 

19 

5 

2 

P(T/n) 

1 

0.879 

0.788 

0.212 

0.061 

US-Light: 

Occurrence 

29 

0 

1 

3 

0 

P(T/n) 

1 

0.121 

0.121 

0.091 

0 

Table  III  ROC  data  of  observer  #2  (Second  Reading) 


A  total  of  30  positive  trails  and  34  negative  trails 


Section(l) 

Category 

1 

2 

3 

4 

5 

Section(2) 

US: 

Occurrence 

11 

0 

4 

2 

13 

P(T/t) 

1 

0.633 

0.633 

0.5 

0.433 

US -Light:  Occurrence 

5 

0 

3 

8 

14 

P(T/t) 

1 

0.833 

0.833 

0.733 

0.467 

Section  (3) 

US: 

Occurrence 

5 

1 

11 

12 

5 

P(T/n) 

1 

0.853 

0.823 

0.5 

0.147 

US-Light: 

Occurrence 

29 

0 

0 

2 

3 

P(T/n) 

1 

0.147 

0.147 

0.147 

0.088 

FIGURES 


INTENSITY  MODULATED  LIGHT 


laser  source  (in  phase)  laser  source  (out  phase 

180  degree) 


o 


laser  detector 


laser  detector 


Fig.l.  Illustration  of  object  detection  by  the  optical  phased  array  positioner 


Fig.  2.  Optical  phased  array  geometry 


Y 


A 


f  (b) 

ultrasound 
probe 


Fig.3.  (a)  Coaxial  ultrasound  and  optical  probe  and  optical  detector  assembly,  (b)  front 
view  of  the  coaxial  probe. 
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response  and  (b)  phase  response.  The  amplitude  null  and  phase  transition  are  indications 
of  the  object  line-of-sight. 


50  100  150 


Fig. 5.  The  triple  images  taken  at  lateral  positions  16,  17,  18  mm.  Target  is  indicated  by 
arrow. 
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50  100  150 

Fig. 6.  The  same  triple  images  as  shown  in  Fig. 5  but  optical  line-of-sight  information  was 
coded  in  the  images  as  indicated  by  dashed  lines.  Target  is  indicated  by  arrow. 
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Amplitude  (volts) 


Fig.7.  Optical  response  when  a  4.8  mm  white  Teflon  rod  was  inserted  into  the  scattering 
phantom,  (a)  Amplitude  response  and  (b)  phase  response.  The  amplitude  null  and  phase 
transition  are  indications  of  the  object  line-of-sight. 
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Fig.  8.  The  triple  images  taken  at  lateral  positions  20,  21,  22  mm.  Target  and  artifact  are 
indicated  by  arrows.  Object  is  a  4.8  mm  white  Telfon  rod. 


3  1 


I  • 


•  % 


50  100  150  ■ 

Fi<^.9.  The  same  triple  images  as  shown  in  Fig. 8  but  optical  line-of-sight  infonnation  was 
coded  in  the  images  as  indicated  by  dashed  lines.  Target  and  artifact  are  indicated  by 


arrows. 
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sensitivity 


Fig.  10.  Average  ROC  curves  of  two  observers.  The  ROC  curves  of  ultrasound  method 
are  plotted  in  solid  lines  and  ROC  curves  of  ultrasound-light  are  plotted  in  dashed  lines. 


4  I 


ROC  curves  of  observer  1 


(a) 


Fig.ll.  ROC  curves  measured  with  two  different  test  procedures,  (a)  Obtained  from 
observer  1,  and  (b)  from  observer  2.  Procedure  #1  is  the  ordered  test  with  ultrasound  and 
ultrasound-light  images  presented  in  sequence  while  procedure  #2  is  the  random  test  with 
ultrasound  and  ultrasound-light  images  presented  in  random  order.  The  results  with  both 
testing  procedures  are  consistent. 
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